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Physical Limitations in Electron Ballistics* 


By J. R. PIERCE 


INTRODUCTION 


HE subject of this talk is “Physical Limitations in Electron Bal- 

listics”. It is pleasant to have a chance to talk about such physical 

limitations, because there is so little we can do about them. And, 
although these limitations are apt to be discouraging, a knowledge of them 
is very valuable, for it keeps us from spending time trying, like the in- 
ventors of perpetual motion machines, to do the impossible. 

As electron ballistics is particularly subject to physical limitations, there 
are so many that it is impossible to discuss all of them thoroughly at this 
time. Also, many of the limitations are of a rather complicated nature, 
and to deduce them from basic principles in a quantitative way requires 
much thought and patience. I think the best I can do is to try to mention 

* most of the chief limitations, as a warning to the uninitiated that rocks lie 
ahead in certain directions, but to concentrate attention on only a few of 
them. I have chosen this evening to devote particular attention to lim- 
itations that bear on the production and use of electron beams in which 
considerable current is required, such as those used in cathode ray tubes 
and high-frequency oscillators, and to mention only briefly as a sort of 
introduction problems pertaining more closely to low-current devices such as 
electron microscopes. 

THE WAVE NATURE OF THE ELECTRON 

One of the most important limitations in electron microscopy is the dual 
nature, wave and corpuscular, of the electron. Without making any 
attempt to justify or explain the combination of wave and particle con- 
cepts which is characteristic of modern physics, we may describe its con- 
sequence at once; very small objects don’t cast distinct shadows. This 
cannot be explained merely in terms of the physical size of the electron and 
the object. When an electron beam is reflected from a surface of regularly 

* A lecture given under the auspices of the Basic Science Group of the New York Sec- 


tion of the A.I.E.E., as a part of an Electron Ballistics Symposium, Columbia University, 
March 21, 1945. 
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spaced obstacles (the atoms in a crystal lattice, for instance) diffraction 
patterns are obtained, similar to those which may be obtained with waves 
of X-rays or light. It appears that electrons get around sufficiently smal] 
objects just as sound waves get around telephone poles, automobiles, and 
even houses, and if the objects are sufficiently small their effect on the 
electron flow will either be absent or will consist of a few ripples which are 
meaningless in disclosing the shape or size of the object. 

The electron wave-length, which varies inversely as the momentum of the 
electron, may be simply expressed in terms of the energy V in electron volts. 
A simple non-relativistic expression which is only 5% in error at 100,000 
volts (a high voltage for electron microscopes), is* 


N= V7150/V X 10° cm (1) 


Thus for 30,000-volt electrons the wave-length is 7 X 107! cm or about 
1.4 X 10°‘ times the diameter of a hair and 1.2 X 10~ times the length 
of a wave of yellow light. 

In terms of this wave-length \ and the half angle of the cone of rays 
accepted by the objective, a, we can express the distance d between point 
objects which can just be distinguished in an electron microscope. This 
distance is 


d = .61X/sina (2) 


For small values of a 
Qe == T/T (3) 


where f is the well known photographic f number, the ratio of the focal 
length to the lens diameter. We see that, just as with cameras, the smaller 
the f number the better. In electron microscopes a small f enables us to 
distinguish smaller objects. 


ABERRATIONS 


Just as in cameras, the limitation to the f number is imposed by lens 
aberrations. But in electron lenses the aberrations are much more severe. 
Why is this so? Because with electron lenses we have less freedom of design 
than with optical lenses. 

Consider an electric lens. The quantity analogous to the index of 
refraction for light is the square root of the potential with respect to the 
cathode. Now suppose that with a light lens we know the index of re- 
fraction at every point along the axis. Suppose, for instance, that the in- 
dex of refraction is 1 everywhere along the axis except for a space L long 


* The relativistic expression is 


N= (V150/V/V/1 + .98 X 10-*V) X 10°78 cm 
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where it is 2, asin Fig. 1. Our lens may be converging or diverging; strong 
or weak. In the analogous electric case, however, the potential throughout 
the lens space must satisfy Laplace’s equation, and this means that if it is 
specified along the axis it is known everywhere. We can easily see this by 
writing down Laplace’s equation for an axially symmetrical field. 


10 OV nf OV P 
= a rT - = { } 
r or or Oz? 4 
i LIGHT-CONDITIONS OFF 
add AXIS NOT FIXED BY 


CONDITIONS ON AXIS 


ELECTRIC FIELO-FIELO 
OFF AXIS SPECIFIED BY 
POTENTIAL ON AXIS 


T 
{ ° 
v-+/ #(3+ir COS e)ae 
0 
Fig. 1—Contrast between optical and electric focussing conditions 


The field near the axis may be expanded in powers of f 


OV a 
-=art+t::: (35) 
or 
Substituting this into (4), 
10 ) —d J 
— (ar) = 2a = ——_~ 
r Or Os" 
OV —1 Fe V 
—_— = (0) 
or 2 Of 


As a matter of fact, the potential V (2,7) remote from the axis can be 
expressed in terms of the potential V’,(s) on the axis as 


1 ® 
V = [ Vo(s + ir cos 0) dé (7) 
WT 0 


If we could introduce charges into our lens, Laplace’s equation would no 
longer hold and we would have more freedom of design. The methods 
proposed for the introduction of charges comprise the use of free charges 
(space charge) which are largely uncontrollable, and the use of curved grids, 
which do more damage than good. In other words, the cures are worse 
than the disease. 
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Similar limitations apply to magnetic lenses, and in the end we find that 
because of the simplest form of aberration, spherical aberration, best dei- 
inition is achieved in electron microscopes with f numbers of 100 or greater, 
while the f number of a light microscope objective corrected for spherical 
aberration and other defects as well may be around unity. Thus the elec- 
tron microscope is severely handicapped, and this handicap is overcome 
only because electron waves are much less than 1/100 the length of light 
waves. 
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Fig. 2—Approximate relation between beam size and angular spread. 


THERMAL VELOCITIES OF ELECTRONS 


In many electron-optical systems, and particularly in such devices as 
cathode ray tubes, it is desirable to focus an electron beam into a small 
area, SO as to produce a very small spot on a fluorescent screen, or to pass a 
considerable current through a small aperture. We might think at first 
that if our focusing system were good enough, that is, if it had very small 
aberrations, we could focus a current from a cathode of given area into as 
small a space as we desired. This, unfortunately, is not so. The obstacle 
is the thermal velocities of the electrons emitted by the cathode. 

A simple example will show the sort of thing we should expect to take 
place. Figure 2 shows a plane cathode and near to it a positive grid so 
fine as to cause no appreciable deflections of the electrons which pass through 
it. Farther on we have an aberrationless electron lens designed to focus 
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the electron stream at a spot a distance L beyond it. The electrons will 
leave the cathode with some slight sidewise velocity components; so, elec- 
tron paths will pass at several angles through a given point onthe lens. The 
lens will bend these paths approximately equally, and hence we can see that 
at the point where the beam is narrowest it will still have some appreciable 
diameter We. 

Now consider the beam at the lens. Suppose that through a given point 
all the paths lie within a cone of half angle 6. Then the width W, 
is approximately 


W. = 21h (8) 
We can also see that the paths at W, will lie within an angle approximately 

6. = Wi/2L (9) 
Hence we see that approximately 

AW, = We (10) 


In other words, we can have a small spot through which electrons pass 
over a wide angular range, or we can have a broad beam in which all paths 
are nearly parallel, but we can’t have a narrow spot and nearly parallel 
rays. 

We see that the actual width of spot will depend on the thermal veloc- 
ities, which are proportional to the square root of the cathode temperature, 
and on the forward velocity, which is proportional to the square root of the 
accelerating voltage. By using more involved arguments we discover 
that for any point in an electron stream, where the beam is wide, narrow, or 
intermediate, the current in an arbitrary direction chosen as the x direction 
can be expressed*'* 


dj = = forse ttre do, do, do, (11) 
v=m+y +; 

when v> — 2eV /m; (12) 

or dj = 0 (13) 

when v < VW2eV/m (14) 


Here j, is the cathode current density, V is voltage with respect to the 
cathode, T is the absolute temperature of the cathode in degrees Kelvin, 
and %z, vy, and 2, are the three velocity components; dj is the_ element 


* This expression neglects the effects of electron collisions, which may actually make 
the current density smaller. 
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of current density carried by electrons which have velocity components 
about 2, vy, 2%, lying in the little range of velocity dvz dvy doz. 

The reason for restriction (12) is that if an electron starts with zero 
thermal velocity from the cathode, it will attain the velocity given by the 
right side of (12) by falling through the potential drop V. As electrons 
cannot have velocities smaller than this, we have (13) and (14). 

By integrating (11) with appropriate limits we obtain a more specialized 
but very useful expression 


: ‘ 5 L1000V\.. 
JI <jn = fo (1 + - ; ) sine (15) 


For usual values of voltage, unity in the parentheses is negligible, and we 
can say that if all the electron paths approaching a given point in an electron 
beam lie within a cone of half angle 6, the current density 7 at that point 
cannot be greater than a limiting value j,, which is proportional to the 
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Fig. 3—Parameters important in determining spot size in a cathode ray tube. 


cathode current density, to the voltage, to sin?9, and inversely proportional! 
to the cathode temperature. 


Let us see what this means in some practical cases. 


Figure 3 shows a 


cathode ray tube. 


The electron stream has a width W at the final electron 


lens, and is focused on a screen a distance 1 beyond the lens. The half 
angle of the cone of rays reaching the screen cannot be greater than 


sin 0= 6 = W/2L 


Suppose the spot must have a diameter not greater than d. 
current be 7. Then from (15), 
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Thus if for a given spot size we want to increase the spot current, and if we 
are limited to a given cathode current density because of cathode life, 
we must make V larger, W larger or 1 smaller. 

Making W larger increases both lens and deflection aberrations. Making 
i smaller means that for a given linear deflection we must increase the 
angular deflection, and this too tends to defocus the spot. Because of these 
limitations, it is necessary to avail ourselves of the remaining variable and 
raise the operating voltage V. 

Another illustration, perhaps a little more subtle, of the effect of thermal 
velocities, lies in the analysis of the properties of a type of vacuum tube 
amplifier known as the ‘deflection tube”. In such a device, illustrated in 
Fig. 4, an electron stream from a cathode is accelerated and focused by a 
lens and deflected by a pair of deflecting electrodes so as to hit or miss an out- 
put electrode. Such a device may be used as an amplifier. 

Now it is obvious that as the output electrode on which the beam is 
focused is moved farther away from the deflecting plates, a given deflecting 
voltage will produce a greater linear deflection of the beam at the output. 


CATHODE . a. a 
A ? q OUTPUT ft 


euecrncn!” ~ELECTRODE 


LENS DEFLECTING 
PLATES 


Fig. 4—Amplifying tube making use of electron deflection 


As this at first sight seems desirable; it has been seriously suggested not 
only that this be done, but that an elaborate electron optical system be 
interposed between the deflecting plates and the output electrode to amplify 
the deflection. 

The merit of a deflection tube is roughly measured by the deflecting 
voltage required to move the beam from entirely missing the output elec- 
trode to entirely hitting the output electrode, and, of course, moving the 
output electrode farther away or putting lenses between the deflecting 
plates and the output electrode doesn’t reduce this voltage at all. As we 
improve the deflection sensitivity by these means, we simply increase the 
spot size at the same time. Focusing our attention on the beam between 
the deflecting plates, we appreciate at once that the electron paths through 
each point will be spread over some cone of half angle @, and that to change 
from a clean miss to a clean hit we must deflect the electrons through an 
angle of at least 20, regardless of what we do to the beam afterwards. 

Returning for a moment to equation (15), we see that it says the current 
density can be less than a certain limiting value depending on 4. Yet 
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expression (15) was obtained by integrating a supposedly exact expression, 
What does this inequality mean? 

The answer is that for the current to have the limiting value, electrons 
of all allowable velocities must approach each part of the spot from all angles 
lying within the cone of half angle 6. When the average current density in 
the spot is less than the limiting current density, the possibilities are 

(a) Electrons are approaching each point in the beam from all angles, but 
along some angles only electrons which left the cathode with greater than 
zero velocity can reach the spot. 

(b) Electrons leaving the cathode with all velocities can reach the spot, 
but at some portions of the spot electrons don’t come in at all angles within 
the cone angle @. 
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Fig. 5—Relation between nearness of approach to limiting current density and fraction 
of current utilized. 


Thus, we can have less than the limiting current either because electrons 
do not reach the spot with all allowable velocities or from all allowable 
angles. Of course both factors may operate. 

We can easily see how lens aberrations, which we know are present in al! 
electron-optical systems, can prevent our attaining the limiting current 
density. There is a more fundamental limitation, however. It can be 
shown that even with perfect focusing, we must sort out and throw away 
part of the current in order to approach the limiting current density, and 
we can even derive a theoretical curve for the case of perfect focusing re- 
lating the fraction of the limiting current density which is attained to the 
fraction of the cathode current which can reach the spot. Figure 5 shows 
such a curve which applies for voltages higher than, say, 10 volts. 

Usually, the failure to approach the limiting current density is chiefly 
caused by aberrations, and in ordinary cathode ray tubes the current 
density in the spot may be only a small fraction of the limiting value. A 
very close approach to the limiting current density has been achieved in a 
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special cathode ray tube designed by Dr. C. J. Davisson of the Bell Tele- 
phone Laboratories. 

When we become thoroughly convinced that these equations expressing 
the effects of thermal velocities very much cramp our style in designing 
electron-optical devices, as good engineers we wonder if there isn’t, after all, 
some way of getting around them. I don’t think there is. The suggestion 
illustrated in Fig. 6 is a typical example of such an attempt. We know 
that in a strong magnetic field electrons tend to follow the lines of force. 
Why not use a very strong magnetic field with lines of force approaching the 
axis at a gentle angle to drag the electron stream toward the axis? 

An electron off axis traveling parallel to the axis certainly will be dragged 
inward by such a field. The catch is that the field pulls the electron in 
because it makes the electron spiral around the axis. As the beam con- 
verges and the field becomes stronger, the pitch of each spiral decreases and 
the angular speed of each electron increases. Finally, if the field is strong 
enough, all the kinetic energy of the electron is converted from forward 


ELECTRON_«, ts MAGNE TIC! 


Z 4 Bol 


PATH >> x | LINES OF 
— FORCE 


Fig. 6—Reflection of an electron by a magnetic field with strongly converging lines of 
force. 





motion to revolution about the axis; the electron ceases to move into the 
field and bounces back out. It may be some small consolation to know 
that very high-current densities can be achieved by this means, but only 
because in their flat spiralling the electrons approach a spot at much wider 
angles with the axis than the small inclination of the lines of force. 


SPACE CHARGE LIMITATIONS 


In electron beam devices using reasonably large currents, the space 
charge of the electrons is a very serious source of trouble both in compli- 
cating design of the devices and in limiting their performance. 

Let us begin our consideration right at,the electron gun, the source of 
electron flow in many devices such as cathode ray tubes and certain high- 
frequency tubes. Electron guns are sometimes designed on the basis of 
radial space charge limited electron flow between a cathode in the form of a 
spherical cap of radius r, and a concentric spherical anode a distance d from 
the cathode. It can be shown that by use of suitable electrodes external to 
the beam, radial motion can be maintained between cathode and anode along 
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Straight lines normal to the cathode surface. A hole in the anode electrode - 
will allow the beam to emerge from the gun. Because of the change in 
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electron gun. They are plots of a factor called the perveance, which is 
defined as 

P = I/v™ (18) 
(that is, current divided by voltage to the 3/2 power) as a function of @, the 
half angle of the cone of flow, and d/r., the ratio of cathode-anode spacing 
to cathode radius. In getting an idea of the meaning of the curves, we may 
note that a perveance of 10° means a current of 1 milliampere at 100 volts. 
It is obvious from the curves that to get very high values of perveance, that 
is, high current at a given voltage, @ must be large and the cathode-anode 
spacing must be small. Making 6 large means that electrons approach the 
axis at steep angles; aberrations are bad and the beam tends to diverge 
rapidly beyond crossover. Moving the anode near to the cathode means 
that the hole which must be cut in the anode to allow the beam to pass 
through must be large, and cutting such a large hole in the anode defeats 
our aim of getting higher perveance; we can’t pull electrons away from the 
cathode with an electrode which isn’t there. Further, for ratios of spacing 
to cathode radius less than about .29, the lens action of the hole in the 
anode causes the emerging beam to diverge, which would make the gtin 
unsuitable for many applications. 

When we build guns for small currents at high voltages, such as cathode 
ray tube guns, space charge causes little trouble; when we try to obtain 
large currents at lower voltages, we find ourselves seriously embarrassed. 

Suppose we now turn our attention to the effect of space charge in beams 
when the beam travels a distance many times its own width. Consider, 
for instance, the case of a circular disk forming a space charge limited 
cathode. Suppose we place opposite this a fine grid, and shoot an electron 
stream out into a conducting box, as illustrated in Fig. 9a. We immediately 
realize that there will be a potential gradient away from the charge forming 
the beam. In this case, the gradient will be toward the nearest conductor; 
that is outwards, and the electron beam will diverge. 

How can we overcome such divergence? One way would be to arrange 
the boundary conditions in such a fashion that all the field would be di- 
rected along the beam instead of outwards; this might be done by sur- 
rounding the beam by a series of conducting rings and applying to them 
successively higher voltages as in 9b, the voltages which would occur in 
electron flow between infinite paralle] planes with the same current density. 
Another way in which the same effect may be achieved is through use of 
specially shaped electrodes outside of the beam, as shown in Fig. 9c."' In 
maintaining parallel flow by these means, the electric field due to the elec- 
trons acts along the beam, and increases continually in magnitude with 
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distance from the cathode. We can in fact calculate the potential at any [ 


distance along the beam by the well known Child’s law equation 


2.33 X 10° AV*?/x* 
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Fig. 9—Avoiding beam divergence by means of a longitudinal electric field. 


Here V is the anode voltage, « the cathode-anode spacing, J the current 
in amperes and A the cathode area, 
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Then 
V = 5,700 volts 


Thus to maintain parallel motion of the modest current of 10 milliamperes 
spread over an area of one square centimeter requires 5,700 volts. More- 
over, the requirement of distributing this voltage smoothly along the beam 
would make it very difficult to put the beam to any use. 

One means for mitigating the situation is to use an electron lens and 
direct the beam inward. Of course, the beam will eventually become par- 
allel and then diverge again, but by this means a fairly large current can be 
made to travel a considerable distance. Some calculations made by Thomp- 
son and Headrick” cover this type of motion, with an especial emphasis on 
the problem in cathode ray tubes, in which the currents are moderate. 

In order to confine large currents into beams, an axial magnetic field is 
sometimes used, as shown in Fig. 10 Here a cathode-grid combination 
shoots a beam of electrons into a long conducting tube. A long coil around 
the tube produces an axial magnetic field intended to confine the electron 
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Fig. 10 dui beam divergence by means of a longitudinal magnetic field 














paths in a roughly parallel beam. The radial electric field due to space 
charge will cause the beam to expand somewhat and to rotate about the 
axis. As the magnetic field is made stronger and stronger, the electrons 
will follow paths more and more nearly straight and parallel to the axis. 
For a given current and voltage, there is one sort of physical limitation in 
the strength of magnetic field we need to get a satisfactory beam. It is 
another effect that I wish to discuss. 

Suppose we have a very strong magnetic field, in which the electrons 
travel almost in straight lines. We know, of course, that the radial electric 
field is still present, and this means that the potential toward the center of 
the beam is depressed; this in turn means that the center electrons are 
slowed down. This slowing down of course increases the density of electrons 
in the center of the beam. The result is that if for some critical voltage or 
speed of injection we increase current beyond a certain value, the process 
runs away, the potential at the center of the beam drops to zero, and another 
type of electron flow with a “virtual cathode” of zero electron velocity at 
the center of the beam is established. Thus, although the magnetic field 
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has overcome the diverging effect of the space charge, we still have a space 
charge limitation of the beam current. C. J. Calbick has calculated the 
value of this limiting current.’ If the beam completely fills a conducting 
tube at a potential V’ with respect to the cathode, the limiting beam current 
is independent of the diameter of the beam and is 


I = 2993x 10°" (20) 


If the beam diameter is less than that of the conducting tube, the limiting 
current is lower. 

But perhaps we can completely overcome the effects of space charge. 
Suppose we put a very little gas in the discharge space. Then positive ions 
will be formed. Any tendency of the electronic space charge to lower the 
potential and slow up the electrons will trap positive ions in the potential! 
minimum and so raise the potential. Thus the gas enables us to get rid 
of the the slowing up effect of the space charge as well as its diverging 
effect. 

Before we congratulate ourselves unduly, it might be well to make sure 
about the stability of an electron beam in which the electronic space charge 
is neutralized by heavy positive ions. Langmuir and Tonks, in their work 
on plasma oscillations, introduced a concept, extended later by Hahn and 
Ramo, which enables us to investigate this problem. The concept is that 
of space charge waves. It is found that in a cloud of electrons whose net 
space charge is neutralized by heavy, relatively immobile positive ions, 
small disturbances of the electron charge density produce a linear restoring 
force; and this, together with the mass of the electrons, makes possible a 
type of space charge wave which may be compared roughly with sound 
waves, although much of the detailed behavior of space charge waves is 
quite different from that of sound waves. We may express a disturbance in 
an electron beam in terms of these space charge waves and then examine the 
subsequent history of the disturbance as a function of time. This has been 
done" and the perhaps surprising result is that even when the electronic 
space charge is neutralized by heavy positive ions, the flow tends to collapse 
if the current is raised above a limiting value 


I = 6x UV (21) 


It is true that this current is 6.5 times the limiting current in the absence of 
ions, but it is a limit nevertheless. 

If this limit in the presence of ions seems unnatural, perhaps we should 
recall a mechanical analogy. Consider a vertical long column subjected to 
aload /. If we subject it to a sidewise force aF proportional to F, as shown 


in Fig. Ila, the behavior on increasing F will be a gradual deformation 
(analogous to the space charge lowering of potential in the absence of ions) 
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ending in collapse. However, even if, as in 11b, there is no sidewise loading 
and no bending during loading, we know from Euler’s formula that beyond a 
certain loading the column will still collapse. This behavior is analogous to 
that of an electron beam in which the electronic space charge is neutralized 
by positive ions and there is no depression of potential in the beam. 

This space charge limitation either in the presence or absence of ions 
allows the passage of quite a large current through a tube, as the table 


below will show: 


Voltage Current, amperes, no ions Current, amperes, ions 
1000 927 6.01 
100 .029 .190 
10 .009 .060 


We might therefore feel that the space charge is disposed of in a practical 
sense, and so it is in many cases.* 


F F 
AF 
T=29.3 x1i0 Vv”? I=190 xio Vv”? 
Qa b 


Fig. 11—-Comparison of limiting stable beam currents with and without positive ions. 
Power DisstpaATiON LIMITATIONS 


Having talked about various limitations imposed by wave effects, aber- 
rations, thermal velocities and space charge on the electron flow in the 
beam itself, I want to close by discussing briefly a topic which seems hardly 
included in electron ballistics but yet is vital to any application in that 
field. I refer to the problems associated with power dissipation when 
electrons strike something and stop. This is a good deal like the problem 
imposed by suddenly coming down to earth while studying the sensations 
of a free fall. It is inevitable and may be fatal unless satisfactory provision 
is made for the dissipation of kinetic energy. 

What I want chiefly to bring out are the consequences of scaling a given 
electronic device down in size. If we change the size of each part of an 

* Tt appears that in many gas discharges, including those in which plasma oscillations 
are observed, the current is too high to allow persistence of the homogeneous flow upon 
which the plasma oscillation equations are based. 
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electron device in the ratio R, if we keep all voltages the same, and if we 
change all magnetic fields in the ratio 1/R, electron current will remain the 
same (provided the cathode is still capable of giving space charge limited 
emission). Electron paths will remain exactly similar, though smaller; 
the power into the electron beam will remain the same, but what will happen 
to the power dissipation capabilities of the device and what will happen to 
the temperature? 

In a device cooled by radiation alone and with cool surroundings, the 
radiating area varies as R®, and since the radiation per unit area varies as 
T*, the temperature will vary as R*. 

In considering a case of cooling by conduction alone, think of a rod 
carrying a certain amount of power away. If all the dimensions of a rod 
are changed by a factor R, the length will be changed by a factor R, the cross 
sectional area will change by a factor R’, and if the thermal conductivity 
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Fig. 12—Variation of magnetic field and temperature in scaling an electronic device. 





remains constant the temperature will vary as R-!. This is a faster rate of 
variation than in the case of cooling by radiation, and hence as the system 
is scaled to a smaller and smaller size, cooling by conduction will become 
negligible and radiation cooling only will remain effective and will determine 
the temperature. 

Figure 12 gives an idea of the variation of various quantities discussed. 

We want to make electronic devices smaller for a number of reasons; 
perhaps chiefly to reduce transit time and so to secure operation at higher 
frequencies. In doing this, we encounter the fundamental limitation of 
reduced power dissipation capabilities and increased temperature. What 
is the trouble? We have scaled everything. Or have we? The answer 
is, we have not. The electrons, atoms, and quanta are still the same size. 
Had we been able to scale these, we should have increased the heat con- 
ductivity and the radiating power of our device, and all would have been 
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well. As it is, if we make a tube for given power smaller and smaller, using 
the most refractory materials available we eventually reach a size of tube 
which will, despite our best efforts, melt, thaw, and resolve itself into a dew. 


CONCLUSION 


Perhaps after these somewhat gloomy words concerning physical lim- 
itations in electron ballistics, you may wonder how it is at all possible to 
surmount the difficulties mentioned. It certainly is not easy; all electronic 
devices represent compromises of one sort or another between fundamental 
physical limitations of electron flow on the one hand and structural com- 
plications on the other. In working with vacuum tubes one is perhaps 
troubled more by physical limitations, difficulties of construction, inade- 
quacy of materials and the lack of quantitative agreement between compli- 
cated phenomena and relatively simple theories than in any other part of 
the electric art. It is for this reason that a friend of mine twisted an old 
aphorism into a new one and said, ‘“‘Nature abhors a vacuum tube”’. 
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Electron Ballistics in High-Frequency Fields* 


By A. L. SAMUEL 


HIS, the final lecture of a series on Electron Ballistics, is not a summary 

of the material which has been previously presented but rather it is an 
attempt to show how the ballistic approach can be extended to the analysis 
of high-frequency devices. Much that might otherwise be said about ultra- 
high frequencies cannot be said because of secrecy requirements. However, 
there is considerable material which can be presented, within the limits of 
the necessary security regulations, which may be of interest to those who are 
not already well acquainted with the subject. I will, perforce, not be able 
to say anything specific about actual devices utilizing the principles to be 
discussed. 

Many of the ultra-high-frequency devices which have come into use 
during the last few years have employed electron beams of one sort or 
another. These devices can be analysed in any one of a number of ways. 
For example, we can write the equation of space-charge flow. This ap- 
proach considers the electric charge as a continuous fluid subject to Poisson’s 
equation. The small-signal theory of Peterson and Llewellyn is an example 
of this type of analysis. Or if we wish we can consider the various types of 
wave motion which can exist in a space-charge region. The space-charge- 
wave analysis of Hahn and Ramo as applied to velocity-variation tubes is an 
example of this. In addition there is an electron-ballistic approach to the 
problem and it is with this method that we will be concerned in the present 
lecture. 

Before we become involved in the details of the analysis, we should perhaps 
spend a few moments considering the relationship between these various 
methods. If we have an interaction taking place between electric fields 
and moving charges, we know at once from Newton’s second law that the 
forces acting on the electrons must of necessity be equal and opposite to 
= : those acting on the fields. It is therefore a matter of small concern whether 
we consider the forces acting on the electrons and the effects of these forces 
on the electron motion or whether we consider the alteration in fields which 
the electron motion produces. We can, if we wish, compute the energy 
transfer to an electric field by the motion of an electric charge or we can 
compute the change in energy of the electron which accompanies this trans- 








* Originally presented on April 11, 1945 as the concluding lecture of a symposium on 
Electron Ballistics sponsored by the Basic Science Group of the American Institute of 
Electrical Engineers. 
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fer. Iwas tempted to say “which results from this transfer’ but this implies 
a cause and an effect, a notion which has no place in the present discussion. 
The dual aspect of any energy-transfer problem must always be kept in 
mind. Much needless discussion frequently arises between proponents of 
one point of view and those preferring the other when the only difference 
is one of language and both groups are really saying the same thing. The 
electron-ballistic approach yields a simple physical picture; it is capable of 
being applied to widely differing situations, but it is not well suited for a 
determination of the reactive contributions of an electron stream. 


Basic CONCEPTS 


There are several concepts which we will find useful in our analysis. 
These concepts are extremely simple, so simple in fact that one is tempted to 
assume that they are well known. However, these concepts are so basic 
to the subject, and their results so far reaching that we must pause to 
consider them. 

The first is the concept of total current, as distinguished from its com- 
ponents. One way of writing Kirchhoff’s second law is 


Div. J = 0 (1) 


This simply says that the total current entering or leaving any differential! 
region in space is zero. This expression must of course be generalized by 
including displacement currents as proposed by Maxwell if applied to 
alternating currents. The current J/ is the total current density as here 
defined. An important consequence of equation (1), actually only an 
alternate way of stating it, is that the total current always exists in closed 
paths. Let us take a simple case of a two-element thermionic vacuum tube 
connected to a battery. Visualize the situation existing if but a single 
electron leaves the cathode and travels to the plate. The electron takes a 
finite time to cross from the cathode to the plate. During this time a current 
exists, the magnitude being given by the relationship 


I = a 


and according to our premise this current is the same in every part of the 
circuit. The current begins at the instant that the electron leaves the 
cathode and it ceases when the electron arrives at the plate. In the appar- 
ently empty region ahead of the electron there must exist a displacement 
component, numerically equal to the conduction, or perhaps we should say 
convection component accounted for by the moving electron. An ammeter, 
were there one sufficiently sensitive and fast, connected in the external leads 
would read a current during this same interval of time. 

I have chosen to talk about but a single electron to emphasize the electron- 
















324 BELL SYSTEM TECHNICAL JOURNAL 


ballistic aspect; however, the concept is much broader than this since it js 

not at all dependent upon a corpuscular concept of the electron. Asa result 

of this property of the total current, the current to any electrode within a 
| vacuum tube does not necessarily bear any relationship to the number of 
electrons which enter or leave it. Obviously then, currents can exist in the 
grid circuit of a three-element tube even though none of the electrons are 
actually intercepted by the grid. This current may have any phase rela- 
tionship to an impressed voltage on the grid so that the grid may draw power 
from the external circuit, or it may deliver power to the external circuit, 
all without actually intercepting any electronic current. The grid-current 
component resulting from the electronic flow between cathode and plate 
may equally well bear a quadrature relationship to the impressed voltage, 
in which case it will either increase or decrease the apparent interelectrode 
capacitance. If these effects seem queer it is because one is still confusing 
the electronic component with the total current. 

A second basic concept once stated becomes self-evident. This is to the 
effect that the only one thing which we can do to an electron is to change its 
velocity, that is, if we are to confine ourselves to the classical concept of 
an electron. We can change its longitudinal velocity, that is, alter its speed 
but not its direction other than possibly to reverse it, or we can introduce a 
transverse component to its velocity, that is, alter its direction as well as its 
speed. Thought of in this light all electronic devices in which a control is 
exercised over an electron stream are velocity-modulated devices. It might 
be argued that one could equally well say that all we can do is to change the 
electron’s acceleration (derivative of velocity) or its position (integral of velocity). 
The singling out of velocity is in a sense arbitrary. It does, however, have 
some very interesting ramifications. 

I might digress for a moment to elaborate on this idea. Since some of 
the newer devices have been labeled velocity-modulation tubes, there is a 
perfectly understandable tendency on the part of the uninitiated to assume 
that these tubes differ from earlier known devices, such as, for example, the 
space-charge-control tubes, the Barkhausen tube or the magnetron in the 
fact that they employ velocity modulation. The real difference lies else- 
where as we shall see in a few moments. At the same time that these newer 
devices were introduced, there was introduced a new way of looking at 
something which is very old in the art. This newer viewpoint, to my way of 
thinking, constitutes a far greater fundamental contribution than do the 
specific devices which have received so much attention. The pioneers in this 
new approach: Heil and Heil, Bruche and Recknagel, the Varian Brothers, 
Hahn and Metcalf, to mention a few, and the many other workers who lost 
in the race to publish their independent contributions in this field—all of 
these people deserve the greatest of praise for their stimulating contributions 
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to our thinking. My only point in all this discussion is to emphasize that 
the basic method of acting on the electron stream has not really been changed 
at all. The entire matter is summarized in the original statement that the 
only thing which we can do to an electron is to change its velocity. 

Before going on to the next aspect of the problem there is a closely related 
concept which should be mentioned. This concept is that a change in the 
component of the velocity of an electron along one space coordinate does not 
introduce components of velocity in directions orthogonal to the first. For 
example, if an electron beam is deflected by a transverse electric field, there 
will be no accompanying change in the longitudinal velocity. The difficulty 
in the way of doing this in a practical case has nothing to do with the concept 
but only with the problem of producing unidirectional fields. Analyses of 
deflecting field problems which ignore the longitudinal components of the 
fringing fields are apt to be wrong. The problem of high-frequency deflect- 
ing fields has been treated in great detail in the literature and frequently 
with more acrimony than accuracy. 

One further note should be added at this point. In an earlier lecture it 
was pointed out that the magnetic effects of an electromagnetic field are in 
general very much smaller than the electric effects. We will not stop to 
prove that this is still true at the frequencies which now interest us but will 
accept it without further discussion. 

For our next concept we leave electron flow for a moment and consider the 
fields within a resonant cavity. You may very properly object that this 
has nothing to do with electron ballistics, and indeed it does not. However, 
we will find it necessary to discuss problems involving cavity resonators, and 
a failure to understand some of the properties of these circuit elements can 
cause a great deal of trouble. There are two conflicting approaches to this 
problem which I will attempt to reconcile. 

The physicist when first presented with the problem of a resonant cavity 
is inclined to say: This is a boundary value problem. The solution consists in 
writing Maxwell’s equations subject to the conditions that the tangential com- 
ponent of E must be zero along the conducting walls. While a scalar and a mag- 
netic vector potential can be defined, the field is not related to the former in the 
simple manner used in electrostatic problems. 

The engineer, on the other hand, is inclined to say: This looks like an 
extension of the usual resonant circuit. A capacitance exists between the top 
and bottom walls of the cavity; charging currents will flow through the single 
turn toroidal inductance formed by the side walls. I would like to know 
what voltage difference exists between the top and bottom walls, and what 
currents exists in the side walls. 

Now, actually, I am maligning both the physicist and the engineer by my 
statements; nevertheless, there are these two approaches. Which is cor- 
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rect? Well, they both are. It is not correct to speak of an electrostatic 
potential within a resonant cavity; nevertheless, we may and do talk about 
the voltage between the top and bottom of a resonant cavity. What do 
we mean? Simply the maximum instantaneous line integral of the electric 
field taken along some specified path. In any practical device utilizing 
electron beams we are naturally interested in the path taken by the elec- 
trons. The fact that the line integral is different for different paths is of no 
great concern. We are interested in but one of these paths. We shall 
therefore have occasion to talk about voltages in cavities but we must always 
remember what is meant, and we must never for one instant forget that this 
voltage is not unique but that it depends upon some assumed path. 

The second peculiarity of this voltage must also be emphasized. The line 
integral must be taken at a specified instant in time. In effect one takes a 
photograph of the field at some instant in time and then at one’s leisure 
performs the integration. 

Now, of course, an electron when projected through such a cavity will 
perform yet another type of integration. The change in squared velocity 
of the electron as expressed in volts will be given by the line integral of the 
field encountered by the electron; that is, integrated not instantaneously 
but with the electron velocity. This is not a simple process, because the 
electron velocity is continuously being changed by the field interaction and 
therefore the velocity with which the integration is performed depends 
upon the integrated value of the field up to the point in question. This 
has nothing to do with the concept of voltage in a resonant cavity. The 
cavity voltage can, however, be considered as the maximum change in 
squared velocity expressed in volts which an electron could receive if its 
entrance velocity was very large so that the transit time was small compared 
with the period of the cavity field. 

The four basic concepts which I have chosen to recall to your mind are, 
by way of summary: (1) the total current is the same in all parts of a circuit, 
that is div. J = 0; (2) the only way we can act on an electron is to change its 
velocity; (3) the changes in the velocity component of an electron along 
any one rectangular coordinate have no effect on the velocity components 
along any other coordinate; and (4) for convenience, a voltage can be defined 
in a resonant circuit as the line integral of the electric field taken along some 
prescribed path. 


TRANSIT ANGLE 


Since we are to deal with the interaction of electrons and high-frequency 
fields, we frequently find it convenient to measure electron velocity not 
directly but in terms of the equivalent potential difference through which an 
electron must fall to obtain the velocity in question, and the unit of measure 
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will be a volt. Instead of measuring the time required for an electron to 
traverse any given distance in seconds, it is also convenient to use, as a 
unit of time, one radian of angle at the operating frequency. We frequently 
refer to the transit angle of an electron rather than the transit time, although 
both terms are used. In fact, we may on occasion measure distances in 
terms of transit angle, and this usage is extended to measure dimensions 
transverse to the direction of travel of the electron beam. When used in 
this fashion, we mean that the dimension in question is such that were an 
electron to be projected in this direction with a velocity equal to that of the 
electrons in the main beam, the high-frequency field would change through 
the stated number of radians during the transit time. 
THE Five FUNCTIONS IN AN ELECTRONIC DEVICE 

With this preliminary discussion out of the way we can now answer the 
question which has probably been troubling quite a few of you. If the only 
thing we can do to an electron is to change its velocity, then in what basic 
way does the velocity-modulation tube differ from the conventional negative 
grid tube or from the magnetron? 

Well, this is an involved story. If we are to make any use at all of an 
electron beam we must in general perform five distinct operations or func- 
tions. First we must produce the beam. Then we must impress a signal 
of some sort onto the beam. From what I have just said this can be done 
only by varying the velocities of the electrons contained in the beam. The 
third operation consists in converting this variation into a usable form. 
It is in this way that the diverse forms of electronic devices differ to the 
greatest degree. We will go into this matter in more detail shortly. The 
fourth operation consists in abstracting energy from the beam, and the final 
operation consists in collecting the spent electrons. While these operations 
are distinct from an analytical point of view, in many actual devices they 
are performed more or less simultaneously and more than one operation 
may be performed by certain portions of the tube structure. In fact, in 
some devices, for example in the space-charge-control tube, the confusion 
is so great as to make the separation seem rather forced. This very confu- 
sion may partly explain why vacuum-tube engineers who were steeped in 
the art were so slow to realize the advantages of this new way of looking 
at things which I will call the velocity-modulation concept. 

By way of mental exercise in this new way of thinking let us see how 
we can analyze a simple space-charge-control triode. Well, first of all we 
have to identify the electron gun which produces the beam. The electrons 
most certainly come from the cathode, but where is the first accelerating 
electrode? Actually there isn’t any unless we think of the combined d-c 
field resulting from the d-c potentials on the grid and plate as assisted by 
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the initial emission velocities as performing this function. The next func- 
tion, that of varying the electron velocities, is performed by the grid which 
varies the potential gradient in the vicinity of the cathode and hence the 
velocity of the electrons as they approach a potential minimum or virtual 
cathode which is formed a short distance in front of the cathode by the 
action of space charge. This virtual cathode performs the third function, 
that of conversion, by sorting out the electrons and allowing only those elec- 
trons with emission velocities greater than some specific value to pass. 
This, then, is one of the conversion mechanisms which we will call virtual- 
cathode sorting. In this example the virtual cathode occurs very close to 
the real cathode but this is not always the case. The fourth function, 
that of utilization, is performed by allowing the sorted electrons to traverse 
an electromagnetic field between the virtual cathode and the plate. This 
operation is completed by the time the electrons have reached the plate. 
Of course in the triode the plate then performs the final operation, that of 
collecting the spent electrons and dissipating the remaining energy as heat. 
It should be clearly realized, however, that this last function need not neces- 
sarily be performed by the same electrode which provides the output field. 
Indeed the so-called inductive-output tube proposed by Haeff is a space- 
charge-control tube in which these two operations are separated. 


CONVERSION MECHANISMS 


But now to get back to a cataloguing of the different kinds of conversion 
mechanisms. ‘The first general type involves sorting. The first kind which 
we have mentioned is by virtual-cathode sorting. A second kind of sorting 
might involve deflecting the electron beam in proportion to its longitudinal 
velocity instead of reflecting or transmitting it. Various deflection tubes 
have been proposed from time to time using this mechanism. We shall 
be forced to neglect this phase of the problem this evening because of time 
limitations but those of you who are interested will find the literature filled 
with detailed discussions. Still a third type of sorting, sometimes called 
anode sorting, is used in certain Barkhausen tubes when the plate is oper- 
ated at or near the cathode potential so that fast electrons are collected while 
slow electrons are reflected and caused to retraverse a high-frequency field. 
There are still other types of sorting mechanisms but I will not burden 
you with these. 

A second general type of conversion mechanism I will call bunching, to 
distinguish sorting in which electrons are separated according to their 
velocities from bunching in which electrons of differing velocities are brought 
together. Now it just happens that many of the older devices used sorting, 
while many of the newer devices use bunching but this is not universally 
the case. For example, the magnetron as used at high frequencies and the 








SERRA N SS Ra Ras Ntestetsn a 































Sarna we at cae sek ea 






ier Dns 2 


















ETNA SANS PRRs Neat et 


; 
i 
f 
: 
q 
% 


ELECTRON BALLISTICS IN HIGH-FREQUENCY FIELDS 329 


cyclotron both employ a combination of sorting and bunching. A peculiar 
property of the motion of an electron in a magnetic field lies in the existence 
of the so called Larmor frequency. You will recall that the angular velocity 
of an electron in a magnetic field depends only upon the field-strength and 
not at all upon the electron’s linear velocity. This time in seconds is 
given by 


1 = 0357 _X 10 . 
H ; 


or in radians 


10600 
AH 


= 2 


Electrons of widely differing velocity can thus revolve together in spoke- 
like bunches with the faster electrons going around larger circles than the 
slow ones, but just enough larger to keep them together. This, then, is 
one kind of bunching, which for simplicity we shall call magnetic bunching. 
It is used in the magnetron and in the cyclotron. We will have more to say 
on this subject a little later. 

A second kind of bunching was used in some of the early Barkhausen 
tubes where the plate electrode was operated at a fairly high negative poten- 
tial so that none of the electrons were able to reach it. Under such condi- 
tions a uniformly spaced stream of electrons with varying velocities is re- 
flected as a bunched stream, the slower electrons being reflected almost at 
once and the faster electrons penetrating the retarding field for a greater 
distance and hence taking longer to return. This same type of bunching is 
used in a newer form of oscillator, commonly referred to as a reflex tube 
which was suggested by Hahn and Metcalf in 1939, and by others at about 
the same time. The reflex tube differs from the Barkhausen tube, not in 
the basic mechanisms so much as in the fact that the conversion mechanism 
occurs in a different region in the tube from the region devoted to velocity 
modulation and to energy abstraction. A second kind of bunching is then 
reflex bunching. 

A third type of bunching was used in the diode oscillators of Muller and 
of Llewellyn. The mathematical research done by W. E. Benham may be 
mentioned as of interest in this connection. In these tubes a uniform stream 
of electrons becomes bunched simply through the fact that faster moving 
électrons overtake slower ones which precede them. In these earlier forms 
of tubes we again have the case where this conversion is performed simul- 
taneously with one or more of the other processes so that it is very difficult 
to separate them. However, in 1935 Heil and Heil;proposed a tube in 
which the conversion region’ was separated from the other regions of the 
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tube. This tube, the velocity-modulation tubes of Hahn and Metcalf, and 
the klystron tubes of the Varian Brothers, are alike in their use of transit- 
time bunching in a relatively-field-free drift tube. Since this separation of 
functions renders these devices much easier to analyze and since the struc- 
tures are quite interesting in any case we will spend most of our time con- 
sidering them and will, I fear, rather neglect some of the other types of tubes. 

We will, of course, keep our analysis as general as possible so that the 
results may be applied to a variety of different devices. 


Input GAP ANALYSIS 


Let us begin by a small-signal consideration of a uniform electron stream 
entering a region in which there is a longitudinal field defined as some func- 
tion of the distance and of time. This can be the entire Llewellyn diode 
or it can be the input region of a klystron. We ask ourselves with what 
velocity will the electrons leave this region and what will be the net exchange 
of energy between the electrons and the field. At any point within the field 
a typical electron will experience an acceleration given by 


1 
} = 
" 


where 7 is proportional to the maximum amplitude of the h.f. field, but con- 
tains a numerical constant so that ¥V is expressed in centimeters per second 
per second. Now in the usual case f(t) will be a simple sine function but 
f(y) may assume a variety of forms. Again, by way of simplifying our 
work we will assume that it is also a sine function. Let us consider how 
we can go about solving this apparently simple equation. Unfortunately 
this expression can not be solved directly because the value of ¢ at any plane 
(that is, the time of arrival of an electron at this plane) depends upon the 
interchange of energy between the electron and the field. Here we are 
forced back to the time-honored mathematical device of assuming a solu- 
tion in the form of a series and then evaluating these coefficients. There isa 
large number of ways in which this can be done, and consequently a large 
number of different solutions which look very different but which all give 
comparable answers. Usually when such solutions are published, the arith- 
metical work is omitted leaving one with the feeling that there is something 
involved that is not within the ken of ordinary mortals. The fact is that 
the work is usually extremely tedious but actually very simple. It will be 
instructive to follow through one form of such an analysis in just enough 
detail to see the amount of work involved. 

Since we are interested in the energy which is proportional to y? we will 
write at once 


(W)ya = K = Kot Kit Ke + Az +... 


E+ nfly) f(t) (1) 
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where the K’s are a function of the transit time, of the field distribution and 
of the entrance phase, and we will proceed to evaluate these coefficients. 
The average energy per unit of change as expressed in volts is then simply 


at the end of the field while the gain is: 


((/2)K — K,) = Set 4 


where the bar means that we are averaging over all values of the entrance 
phase. 

It is of interest to evaluate the value of velocity \? which individual elec- 
trons receive as a function of the entrance phase. For small signals it is 
usually sufficient to evaluate ¥’ maximized with respect to the starting 
phase, then 


ntK, , 1 tKe 


2 


“ 


oe 


V won _ (g 2K — Ko) max _ | 


a --- |max 


We can further define the ratio of Vinax to the largest value it can have as a 
coefficient 8, sometimes called the modulation coefficient. 

But now to evaluate the A’s. There are many ways of doing this as I 
have intimated. We will proceed by writing 


y = yot) + nyr(t) + mPye(t) + nays(t) + 


where the y’s are coefficients depending upon the transit time ¢ which in 
itself is a function of the applied field thus 
t= tot nti t+ te + n°t3 + 
We can then expand each function of time into a series remembering that 
f'(x)d ‘i I" (x) d 
1! 


f(x + d) = f(x) + 


or for our particular case 


‘io(to) [nt ” ty 7 

yo(t) = yo(to) + —— In + 0 Ju ht +**) 

Vo(to){nty + n te a n ts a e. ./ 
2! 

Now we can expand y,(¢), yo(f) etc. in exactly the same way. Finally we 

get a collection of terms which can be grouped in like powers of n thus 


y = Yo(to) + n [terms in 4, ¥, f1, fo, etc.] + 7? | ae 
The coefficient of the 7 is in fact Vo(to) 4: + yilto). We will not bother to 
write the rest. This expression can then be differentiated to get y and then 
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squared. However, we still have some undetermined coefficients the /, , 
tgetc. terms. These we can evaluate by noting that we wish these values at 
y = a, where a is a fixed distance in the actual device. 
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At this distance 


the ¢ coefficients in the expression for y must have such values that the value 


of y does not change with the value of 7. 


This can only be true if the 


individual expressions multiplying each power of n are each equal to zero. 


Equating these expressions to zero one can evaluate all of the ?’s. 


ple the first term yields 


or 


Jo(to)ti + yullo) = 0 


_ illo) 


stele)” 


For exam- 


Introducing these values, differentiating and squaring, one finally gets an 
expression for (*)y = 0 as a power series in y, the coefficients all being of a 


form easily evaluated for any specified field distribution. 


Since we have by 


definition called these coefficients Ko , Ki, etc. these values are then 


r | 
Ko = Yo 
K, = 2(suii — yor) 
> asf - a? es oN 
Ko = (41 _ 211 +- 240 Ve) inte 20 V2 + ¢ — 
Vo 


This then constitutes the formal solution of the problem. 


We must 


now particularize our problem to some specific field distribution and evaluate 


the y coefficients. 


Suppose, for example, that there is a uniform d.c. field 


(E of equation 1) and an alternating field which varies as some cosine func- 


tion of distance. 


and 


Then the latter is 


f(y) = cos (2 + :) 


y=, E+ 1008 (wt + 9) cos 


( 


b 


ns + :) 


we must eliminate the y which appears in this expression and replace y by 


its equivalent 


and expanding 


Y= Vt mt vet ... 
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cos (2 + :) = cos (= + :) 4. 


Fs T Vo 2 
-sin{ —— + ¢)[nyi + aye + ---] +> 
b b 
1! 
and as before equating like powers of n with ¥) defined as 
JraA+vhit+ Th + 


we finally arrive at 


Ba 
W=-E 
g 
V1 = cos (wt + ¢) cos (= + :) 
. 
V2 => <7) a/b cos (wt = ¢) sin (= + ‘). 


Now we need only integrate these expressions to obtain the values of the j’s 
and the y’s needed to evaluate the K’s. 

If we average y* over all values of the starting phase we can write the 
energy contributed by the field to the electron’s velocity. When this is 
done one finds that the odd powers of n are identically zero leaving only the 
even powers to be considered and for small signal analysis purposes we need 
only consider K,. ‘The energy per electron expressed in volts is 


V = 2.49 X 10 °E*r*4(6) 


where f(@) = w?K2 , and the power is obtained by multiplying this expression 
by the beam current in amperes. 

The end results can be expressed as curves of f(@) against @ as shown in 
Fig. 1. Three examples are shown: the uniform field case and two different 
harmonic distributions as indicated by the smaller plot in the lower left- 
hand corner. You will note that there exist regions of positive f(@) where 
the net transfer of energy is from the field to the electron and regions in 
which the transfer is in the other direction; the former portions are of con- 
siderable interest in connection with the input gaps in velocity modulation 
tubes, and for that matter in the cathode grid region of the negative grid 
tube although this is more complicated than is here indicated, as this trans- 
fer of energy constitutes a loss to the field which loads the input circuit. 
The latter portions may be utilized as was done in the Muller and Llewellyn 
diodes to obtain sustained oscillations. 

If, as I have indicated, we maximize 7” as a function of the starting phase 
we can evaluate the modulation coefficient. The value for the uniform field 
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case, as shown in Fig. 2, is simply, 8 = 
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sin 6/2 
2” 


write the loss expression for this case as 


Now let us consider the conversion region in a typical velocity-variation 
tube. Figure 3 is a drawing of several such devices with the conversion 


16 


f(0) = 2 (1 — cos @) — @sin @. 


Drift SPACE ANALYSIS 


For future reference we wil] 





P=2.49x10 €71#(6) 

















VALUE OF f (6) 





Fig. 1—The energy transfer between an initially uniform electron stream and a longitudi- 
nal electromagnetic field as a function of transit angle. 


regions indicated. 





We will assume for the moment that the electrons enter 
this region with a small variation in velocity and at a perfectly uniform rate. 
Since the total number of electrons entering the region must be equal to the 











TRANSIT ANGLE (6) 


number of electrons leaving the region we may write 


or 


1, dt, = lod 
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—_ dto 
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@ IN RADIANS 
Fig. 2—The (velocity) modulation coefficient between an initially uniform electron stream 
and a uniform electromagnetic field as a function of transit angle. 
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Fig. 3—Typical velocity variation devices employing transit-time bunching. 
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However, a relationship exists between ¢; and ty , 


f 


9) 


& 


t, = lo + 
Where 


/ ; 
v=% Vi1l+asinah, 
{ 


th = + - : . 
tw V1+ asinatl; 


Now if a < 1 
t ; 
h=h& + - (1 ~ $ sin at ++) 
and 
dty = £ aw 


— = ] cos wt 
dt, — 


‘ : law 
1) = 10 1 -+ cos aly 
Y%-2 


and finally 


so that finally 


5 . 6 
4 = a(t + > cos ot) 


This says that the velocity variation impressed on the beam at the en- 


trance to the drift space or conversion region has resulted in a current varia- 
tion at the output. For those of you who think in vacuum tube parameters 
it is of interest to differentiate this expression with respect to the a-c voltage 
and obtain the transconductance 

diy 
dV ac. 


i| = ig (1 + vet 


diy an Bio 
dVacl 2V° 
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This result is obtained by neglecting all of the higher order terms and is 
therefore only a small signal theory of a very restricted sort. 


Now let us consider what we have done. Well, we have followed a small 


interval of time through the drift tube. At the input this time df had a 
current 7%) associated with it; at the output the size of this unit of time is 
different—it is now dt; and the current associated with it is 7;. The physi- 
cal picture corresponding to this phenomenon is that of a uniform distribu- 
tion of electric charge becoming bunched with time as it traverses the drift 
space. 

The next step in the analysis is to carry our approximation a step further 
and consider higher-order terms. Expanding the expression for 7; and using 
our nomenclature the desired expression is 


; . 6 6 6 
ly =m) 1-+2 n(% Jos wt + Js (2% Jos 2wt + (1 = Jos nat) |. 


This equation is not exact since it neglects space charge effects but it does 
indicate the presence of harmonics in the beam current and it reveals cer- 
tain non-linear effects which can also be illustrated by the so-called phase- 
focusing diagrams of Bruche and Rechnagel. 


PuHase-Focustnc DIAGRAMS 


Bruche and Recknagel pointed out that an analogy exists between the 
focusing in space of a parallel light beam and the focusing in phase of the 
electrons in a uniform electron beam. In fact a small-signal theory can be 
developed entirely in terms of optical equations. We will not go into 
this aspect in detail but we will use their diagram (Fig. 4) to illustrate the 
bunching effect graphically. A uniform beam of electrons is represented by 
a series of parallel lines in distance and time coordinates, focus being indi- 
cated by a crossing of these lines after they have been deflected by the veloc- 
ity modulation. 

This general type of diagram has been popularized in this country by the 
Varians, and their associates under the name Applegate diagram, the only 
difference being an interchange of axis. Figure 5, taken from a recent paper 
by Dr. A. E. Harrison, illustrates this version of the Bruche and Recknagel 
diagram. 

Now if instead of judging the current density by the density of the lines 
on the diagram, we make a plot of the current density as a function of time 
for different fixed distances from the input gap, the pictures are somewhat as 
shown on Fig. 6. Figure 7 represents a plot presented by Kompfner and 
combines in one illustration the type of presentation used by Tombs. 
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PHASE FOCUSING IN A REFLEX TUBE 
It might be well to pause for a moment in our discussion of transit time 
bunching to consider how the phase focusing diagrams can be applied to a 
reflex tube. The elements of a modern reflex tube are shown in Fig. 8 
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Fig. 4—The phase-focusing diagram of Bruche and Recknagel showing the analogy to 
optical focusing. 

which was taken from a recent I.R.E. paper by Dr. J. R. Pierce. Electrons 

from the cathode pass through an input gap defined by two grids where they 

are modulated in velocity. In traveling in the retarding field produced by 

the repeller those electrons which passed the gap when the field was becom- 

ing progressively less accelerated, become bunched; the faster electrons 
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Fig. 5—Applegate’s version of the phase-focusing diagram (Harrison). 
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Fig. 6—The conduction-current distribution at different distances along the beam as 
predicted by thz phase-focusing diagram. 
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Fig. 8—The elements of a modern reflex tube (Pierce). 
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penetrating the field to a greater extent and waiting, as it were, for the 
slower electrons which follow to catch up. The electrons which pass across 
the gap while the field is becoming progressively more accelerating are 
spread out. If the retarding field is uniform it can be likened to the earth’s 
gravitational field and the phase-focusing paths on our time-distance plot are 
parabolas. Figure 9, taken from Pierce’s paper, illustrates this while Fig. 
10 is such a plot taken from the paper by Harrison. One interesting and, 
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Fig. 10-——The phase-focusing diagram for a reflex oscillator (Harrison 








in a way, unfortunate difference between reflection bunching and direct 
transit-time bunching is the fact that for reflection bunching the slow elec- 
trons catch up with the fast ones while the reverse is true for the other type. 
This means that if both types of bunching are present as shown in Fig. 11, 
) (also taken from Harrison’s paper) one will tend to undo the effect of the 
| other. 

Another way of combining effects of separate bunching actions is to build 
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a cascade transit-time-bunching amplifier in which a series of three gaps is 
used together with two drift spaces. The first gap velocity modulates the 
beam; this modulation is converted into a current modulation in the first 
drift space. The beam then excites the second cavity, which again velocity 
modulates the beam in quadrature with the original modulation. This 
action of course occurs in the output gap of a two-gap tube but it is not 
there used. Here this second and larger velocity modulation is converted 
to current modulation in the second drift space. The output is finally 
taken off the beam by the third gap. A phase-focusing diagram of this 
sort (again taken from Harrison’s paper) is shown in Fig. 12. 


SPACE-CHARGE-WAVE ANALYSIS 


This phase-focusing approach is rather intriguing as one feels that one 
has a physical picture of what is going on. The picture is, however, very 
inexact except under certain highly specialized cases, as it completely ignores 
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Fig. 11—Diagram showing reflex bunching combined with field-free transit-time bunching 
(Harrison). 


space-charge effects. These space-charge effects are of two sorts: a d-c 
effect, if you will, and an r-f effect; that is, the presence of the electrons of 
the beam will alter the average velocity of the electrons at different parts of 
the beam, and will tend to undo the bunching action. Because of this 
second effect, the electrons are effectively prevented from passing each other 
as the graphical solution suggests. Instead, as the density of the electrons 
in the bunch becomes greater, the mutual repulsion forces tend to prevent a 
further concentration of charge. The electron bunch then tends to disperse. 
The action could be likened to the propagation of a sound wave in a moving 
column of air. While there are several approximate ways to handle this 
problem, Hahn was the first to propose a really satisfactory theory. Inci- 
dentally it should be noted that the Benham, Muller, Llewellyn and Peter- 
son type of theory is capable of treating this aspect of the problem in a 
rigorous way and including all space-charge effects, but unfortunately these 
theories are limited in that they have been applied only to the parallel- 
plane case, and of course they are only small-signal theories. 
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APs Is — Hahn’s analysis starts by treating an infinitely long electron beam, using 
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: first b motions are small compared to the d-c but it does not ignore the r-f effects 
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Fig. 12—Diagram for a cascade amplifier (Harrison). 
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tric wave guide willdo. He assumes an axial magnetic field and a stream of 
positive ions having the same velocity axially and the same charge density. 
These ions are assumed to have infinite mass. The solution is much too 
complicated and involved to present here even in abstract. It involves the 
complete solution of Maxwell’s equations subjected to the stated assump- 
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tions as restricted by the assumed boundary conditions at the edge of the 
beam. 

It is found that two waves are possible, one traveling slightly faster than 
the electron beam and the second traveling slower. A point where the 
velocity components are in phase will correspond to the input to the beam, 
while points where the current components are in phase correspond to the 
desired positions for the output. The propagation constants for these two 
waves in a simplified special case where the magnetic field strength js 
infinite are given by Hahn, as well as expressions for the optimum drift 
tube length. He goes on to consider the case where the magnetic field is 
zero and finds that for this case the density of the charge does not vary much 
but instead the beam swells in and out so that instead of being lumps of 
charge with spaces between, the lumps appear in the outer boundary. Hahn 
has extended his general method of analysis to consider the modulation 
coefficient of gaps through which the beam must pass. His results are a 
great deal more general than those we have presented. 

Ramo has reformulated Hahn’s theory by means of retarded potentials 
for the most important case. This results in some simplification of the 
theory. He computes the more important design constants for a velocity 
modulated tube, such as the optimum drift tube length and the amount and 
phase of the transconductance. Those of you who are particularly inter- 
ested are referred to the original paper. An interesting aspect brought 
out rather forcibly by Ramo’s analysis is the existence of higher-order waves 
on the beam, always occurring in pairs, one faster and the other slower than 
the beam velocity. 


THE MAGNETRON 


In what time remains I want to say just a very few words about the mag- 
netron. This is a very complicated subject and one which cannot be ade- 
quately dealt with in an entire evening, and certainly not in the time 
remaining. 

As you all know, the magnetron was invented and named by Dr. A. W. 
Hull. Habann, Zacek, Okabe and others pioneered in the use of the mag- 
netron as an ultra-high-frequency oscillator. As envisioned today a 
magnetron is a two-element device, usually cylindrical with a centrally 
located cathode and a surrounding anode. The anode may be continuous 
or it may be split into a number of segments as suggested by Okabe, and 


these segments joined together either externally or internally by resonant 
circuits. 


The basic ballistic problems of the magnetron, and hence the only prob- 
lems which directly concern us at this time are (1) that of determining the 
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electron paths within the magnetron and having determined these paths (2) 
that of getting an understanding of the mechanism whereby electrons in 
traversing these paths are able to deliver energy to the connected high- 
frequency circuits. One might think that the first problem would be a 
relatively easy job. As a matter of fact the literature is surfeited with 
papers purporting to give the answer. Unfortunately almost all of the 
published work ignores the effect of space charge. A few moments’ thought 
will suggest that space charge may be a controlling factor because of the 
long electron paths which are sure to result in crossed electric and magnetic 
fields, and indeed more detailed computations bear this out. Nevertheless 
the neglect of space charge greatly simplifies the problem. There are those 
who believe that the no-space-charge theories have no bearing on the way 
actual magnetrons work and that any correspondence between the predic- 
tions of such theories and the actual behavior of magnetrons is simply the 
result of an unfortunate coincidence. In fact Brillouin points out that 
the simplified form in which the Larmor theorem is applied by many, is in 
itself an approximation which was perfectly valid as originally applied by 
Larmor to the electronic orbits within the atom but which does not apply 
to conditions as they exist in the magnetron. 

A number of recent workers have attempted to.include the effects of 
space charge but have unfortunately largely restricted themselves to small 
signal theories while the magnetron is seldom operated under small signal 
conditions, at least not intentionally. Most theories are further restricted 
to a consideration either of the coaxial case where the cathode radius is 
small compared to the anode radius or of the plane case. Most practical 
structures are intermediate between these extremes. 

As an example of the difficulties involved, Fig. 13, reproduced from a 
paper by Kilgore, shows the electron paths as computed neglecting space 
charge and also shows experimental proof that these paths actually exist. 
This illustration has been frequently reproduced and widely accepted. 
The experimental picture was obtained in the presence of gas, to make 
the electron beam path visible, and unfortunately the ionization which 
makes the beam visible also tends to neutralize space charge effects. The 
experimental arrangement departs still further from reality in that the 
electron emission from the cathode was restricted to a limited region so 
that the space charge forces were still further reduced. Now it is probably 
true that some magnetrons operate with electron paths as shown; still it is 
not true that all magnetrons operate in this way. 

Contrasting with this picture which was until recently commonly ac- 
cepted, Brillouin, Blewett and Ramo, and others have shown that stable 
distributions are possible in which a space charge of almost uniform density 
rotates with a uniform angular velocity about the axis. Brillouin goes so 
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far as to label the curves due to Kilgore as wrong, and pictures the possible 
electron trajectories as shown in Fig. 14. 

One of the earliest papers to consider this newer picture of the electron 
paths in the magnetron was published by Posthumus in 1935. This was 
definitely a ballistic approach and hence suitable for discussing tonight. 


+100 ELECTRON 





Typical electron paths in a two-segment magnetron showing how electrons arrive 
at the plate-half of lower potential (Kilgore). 


Posthumus limits his discussion to but one type of oscillation which can be 
obtained in the split-anode magnetron. Those of you who are familiar 
with the early literature on the magnetron will recall that two distinct types 
of oscillations were frequently described. One type usually called “elec- 
tronic” was found to occur under conditions when the magnetic field was just 
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high enough to cut off the anode current under static conditions. This field 
has the value computed by Hull: 


Hull’s first computation, by the way, was made neglecting space charge, 
but, strangely enough, the result is not changed by space charge. These 
electronic oscillations were assumed to be related in frequency to the time 
of transit of an electron from the cathode to the anode, and at cutoff this is 
inversely proportional to the field strength, as expressed by the empirical 
relationship 


AH = 13,100: 





Electronic trajectories for different magnetic fields 
A—small magnetic field L >>} 
B—moderate magnetic field L ~b 
C—strong magnetic field L<b 
D—critical magnetic field L=0 


Fig. 14—Electronic trajectories for different magnetic fields varying from weak fields to 
the critical field shown to the right (Brillouin). 


In general, it was found that best operation occurred when the magnetic 
field was not quite perpendicular to the electric field. The efficiency and 
outputs as reported for this type of oscillator were always low, in spite of 
the large amount of effort devoted to it by an equally large number of work- 
ers. <A second type of oscillation, usually referred to as negative resistance 
oscillations, has also been the subject of considerable study and some practi- 
cal use has been made of it at relatively low frequencies. 

Contrasting with this, Posthumus described a third kind of oscillation 
which he called rotating field oscillations. As in the electronic oscillations 
the preferred frequency is determined by the magnetic field-strength and the 
anode potential, the frequency being inversely proportional to the magnetic 
field-strength. Contrasting with the electronic oscillations, the rotating 
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field oscillations occur with the magnetic field-strength very much above the 
critical cutoff value and the efficiency on occasion reached as much as 70! 
While a careful reading of the literature will reveal that some of the earlier 
experimenters were occasionally dealing with these oscillations, Posthumus’ 
observations represent a new departure in magnetron theory and practice 
and one which we might do well to investigate. 

Posthumus’ approach consisted in studying the electron paths in a mag- 
netron in detail in order to find the conditions under which electrons may 
reach the plate with considerably less energy than that corresponding to the 
plate potential. He assumed a magnetron having & pairs of plates and 
based his calculations on the supposition of a rotating electric field with & 
pairs of poles. In reality there exists a simple alternating field but this 
can be resolved into two rotating fields rotating in opposite directions. 
Power engineers will recognize this as identical with the procedure used in 
analyzing single-phase rotating machinery. Posthumus neglected the field 
opposite to the static angular velocity and considered only one component. 
This is an approximation but a fairly plausible one which can be partially 
justified. 

In the absence of oscillations there is a radial electric field independent 
of the angular position and inversely proportional to radius (for the coaxial! 
cylindrica! case). When oscillations are present there is an additional radia! 
field which varies as some periodic function of the angle and with a period 
2r, and a tangential component of the same general type. For simplicity 
these functions are taken to be simple harmonic functions and can therefore 
be split into two circular rotating fields. 

Posthumus writes the two simultaneous differential equations determining 
the path of an electron, neglecting space charge, and inquires if a solution 
is possible for an elecron path which travels at approximately the same 
angular velocity as the rotating field but lags it by an angle a. An equally 
satisfactory way of looking at this is to say that we transform our coordinates 
from a fixed system to one rotating with the field and inquire if a solution 
is possible where a the angular motion is always small. He finds that such a 
solution is indeed possible and that for the electron motion to be stable the 

value of a must be such that the electrons are somewhat behind the line for 
which the field has its maximum retarding value. The electrons are thus 
in a position to lose energy to the field and to spiral out toward the anode. 

Posthumus defined the value of the electron’s radia! velocity squared at 
the anode as P and the total velocity squared at the anode as VY. Nor- 
malized plots of these two parameters are shown in Fig. 15 as a function of 
frequency. The upper plot shows the radial velocity. Obviously for elec- 
trons to reach the plate at all they must have a positive velocity at the plate. 
Electrons can therefore reach the plate with any given field value, say Z = 2, 
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the abscissa axis. The line for P = 0 appears on the lower curve as the 
dotted line s. Here the ordinate is the total velocity squared, normalized 
with respect to the value without oscillation. Efficiencies can therefore be 
put on the plot directly as shown by the right-hand scale in per cent. The 
line s is therefore a plot of the maximum possible efficiency. This refers 
to what we might call the electronic efficiency since no account is taken of 
circuit losses. Now in any physical device there are some circuit losses and 
hence a lower value of electronic efficiency for which sustained oscillations 
are not possible. 


The dotted line p is Posthumus’ experimental value for 
this lower limit. 


Between the lines p and s, then, oscillations are possible 
at frequencies given by the abscissae and with field values shown on the 
solid lines. Actual data for an experimental tube are shown on the plot, 
oscillations occurring at the wavelengths indicated and over the ranges in 
field shown by the lines terminating in arrows. 

One additional line ¢ is shown on the plot connecting points on the different 


Z lines for which the efficiency isa maximum. The optimum design would 


be one based on the intersection of this line with the p line. Still other facts 
will appear from a detailed study of these results but we shall not be able 
to devote any more time to this interesting subject. 


CONCLUSION 


In concluding a talk of this sort and particularly in concluding a series of 
talks, it is usually appropriate to look ahead to the future and predict the 
trend of affairs, or perhaps to point out certain fruitful fields of research. 
I find this a singularly difficult thing to do. However, it is not revealing 
any military secrets to say that much of the progress of the last few years 
has been in the direction of making things work and not toward getting a 
clearer understanding of the underlying theory. If, for example, an il- 
luminating approach could be devised which would make the problems 
associated with transverse fields, both electric and magnetic, appear as 
simple and straightforward as do longitudinal-electric-field problems, as a 
result of the velocity-modulation concept, then I believe even more striking 
advances could be made in the ultra-high-frequency field than those which 
the war years have brought forth. 
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Dynamics of Package Cushioning 


By RAYMOND D. MINDLIN 


INTRODUCTION 


ECHANICAL damage is a common occurrence in the transportation 

of packaged articles. The causes of failures are generally inadequate 

protective cushioning, lack of ruggedness of the outer packing container, 

or occasional abnormal weakness of the packaged article. The first of these 
difficulties is the subject of this paper. 

One of the major influences in reducing the incidence of mechanical failures 
of packaged articles in recent years has been the use of the drop test. The 
drop test is performed simply by raising the package to a specified height and 
dropping it to the floor. The package and its contents are then examined 
fordamage. This is a go-no-go test and requires a large number of samples 
before a reliable estimate of quality can be made. An adequate number of 
tests is prohibitive when the article packaged is costly. In such cases it is 
important, and in any case it is useful, to supplement the drop test data with 
measurements and calculations. It is also possible to evolve rational pro- 
cedures for designing packages, as described in the present paper, so that 
a particular product will survive a drop test at any specified height, with a 
known factor of safety and with a minimum amount of space assigned for 
cushioning. The drop test then becomes only a check instead of playing an 
integral role in a cut and try design procedure. 

Assuming that the outer container is adequate, the survival of a packaged 
article in a drop test still depends upon a large number of factors descriptive 
of the mechanical properties of both the cushioning medium and the pack- 
aged item. However, the more important properties can be grouped so 
that they may be replaced by knowledge of only the following factors: 

(1) The magnitude of the maximum acceleration that the cushioning 

permits the packaged item to reach. 

(2) The form of the acceleration-time relation. 

(3) The strengths, natural frequencies of vibration and damping of the 

structural elements of the packaged article. 

Part I of this paper is concerned primarily with methods for predicting 
maximum acceleration of the packaged article with emphasis on non-linear 
cushioning. Part II deals primarily with the prediction of the form of the 
acceleration-time relation. Part III deals with the effect of acceleration on 
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the packaged article and gives methods for determining whether or not th 
strength of the packaged article will be exceeded. The strength deter 
minations themselves are not dealt with here; but the information in Part I/| 
is essential in interpreting and applying the data obtained in strength meas- 
urements. In Part IV some consideration is given to the influence of dis- 
tributed mass and elasticity. 

It cannot be emphasized too strongly that the determination of the mech 
anical properties of the packaged article, not dealt with in this paper, is an 
essential preliminary to a rational design procedure for packaging. The 
whole purpose in designing package cushioning is to limit the forces which 
may act on the packaged item. If one does not know to what values to 
limit the forces, a rational design procedure cannot be applied. 

It is interesting to observe that the methods described here for analyzing 
and designing package cushioning are directly applicable to the design of 
shock mounts intended to protect equipment from the effects of a sudden 
change in velocity. All of the principles, formulas and design curves given 
here may be used in the shock mount problem with the simple substitution 
of V*/2g for h, where h is the height of drop in the packaging problem, g is 
the acceleration of gravity and V is the velocity change in the shock mount 
problem. 

This paper is essentially a report on a study undertaken at the Bell Tele- 
phone Laboratories, Inc., in the Electronic Apparatus Development 
Department. The results have been applied to the packaging of large 
vacuum tubes and all of the examples used to illustrate the analysis and 
design procedures in the paper are taken from vacuum tube applications. 

Miss H. A. Lefkowitz, Member of the Technical Staff, Bell Telephone 
Laboratories, assisted in the mathematical studies. The oscillograms, used 
as illustrations, were prepared under the supervision of Mr. F. W. Stubner, 
Member of the Technical Staff, Bell Telephone Laboratories. Figure 
3.8.2 was taken from a thesis submitted by Mr. C. Ulucay in partial fulfill- 
ment of the requirements for the degree of Master of Science in the Depart- 
ment of Civil Engineering at Columbia University. The calculations for 
Figs. 3.5.1 to 3.5.6 and Fig. 3.2.2 for 8; > 0 were performed on the Westing- 
house Mechanical Transients Analyzer under the supervision of Dr. G. 
D. McCann, Transmission Engineer, Westinghouse Electric and Manu- 
facturing Company. 


ASSUMPTIONS 


The procedures to be described for the analysis and design of package 
cushioning are based on applications of a few simple laws of mechanics to 
an idealized mechanical system representing the package and its contents. 
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DYNAMICS OF PACKAGE CUSHIONING 


Essentially, a package consists of 

1. Elements of the packaged article which are susceptible to mechanical! 

damage. 

2a. The packaged article as a whole. 

2b. A cushioning medium (excelsior, cardboard spring pads, metal springs, 

etc.) 

3. An outer container (cardboard carton, wood packing case, etc.) 
The four major components are illustrated schematically in Fig. 0.2.1. 
The system is further idealized by “Jumping the parameters’; for example, 
the outer container is considered as a single mass, the cushioning is con- 
sidered as a massless spring with friction losses. The result of this idealiza- 
tion is to lose some of the fine detail of the real distributed system such as 
wave propagation through the cushioning and higher modes of vibration in 

















Fig. 0.2.1—Schematic representation of a package. 
1. Element of packaged article 
2a. Packaged article as a whole 
2b. Cushioning 
3. Outer container 


the package structure and in the packaged article. Some consideration of 
these details is given in Part IV. 

The major components 
of the system are as follows: 

1. A structural element of the packaged item is represented by a mass 
(m,) supported by a linear massless spring with or without velocity 
damping. The mass m, is assumed to be small in comparison with 
the mass of the whole packaged item. 

. The whole packaged item is represented by a mass mp. 

. The cushioning is represented by a spring which may have a linear 
or non-linear load-displacement characteristic and which dissipates 
energy through velocity damping or dry friction. Permanent de- 
formation of the cushioning is not considered, that is, in a repetition 
of the drop test it is assumed that the package has the same properties 
as before the first test. A properly designed package will have essen- 
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tially this characteristic. The mass of the cushioning is assumed to 


be small in comparison with me, except in Section 4.2. 

3. The outer container is represented by the mass m;. The impact oi 
m; on the floor is assumed to be inelastic and during contact the rela 
tive displacement between m3 and the initial position of the floor 
is assumed to be small in comparison with the relative displacement 
between ms. and m3. In other words, no spring action is assigned to 
the outer container and the floor is considered rigid. 
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Fig. 0.2.2—Idealized mechanical system representing a package in a drop test. 
PART I 
MAXIMUM ACCELERATION AND DISPLACEMENT 


1.1 INTRODUCTION 


Most of Part I is concerned with the prediction of the maximum accelera- 
tion that the cushioning permits the packaged article (m») to attain. In 
many instances this will be all the information necessary for judging the 
suitability of a cushioning system. It will be all that is necessary if the 
shape and scale of the acceleration-time function satisfy certain criteria 
which are treated in detailin Parts II1and IV. If these criteria are satisfied, 
the effect of the drop on the packaged article is found by multiplying the 
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DYNAMICS OF PACKAGE CUSHIONING 


dead load stresses (obtained in the usual manner) by the ratio of the maxi- 
mum acceleration to the acceleration of gravity. If the criteria for the use 
of maximum acceleration alone are not satisfied, then Parts II and III will 


supply a numerical factor (the Amplification Factor) by which the maximum 
acceleration should be multiplied, and the remainder of the procedure is the 
same as before. 

The determination of the maximum acceleration is founded on a knowledge 
of the load-displacement characteristics of the cushioning. When the cush- 
ioning system is simple enough, the load-displacement relation may be found 
or designed by purely analytical procedures. The tension spring package, 
discussed in Sections 1.7 and 1.8, is an example where such a treatment is 
possible. In many instances, as with distributed cushioning, the load- 
displacement relation is more easily found by test. 

A load-displacement test is made by applying successively increasing 
forces, with weights or in a load testing machine, to the packaged item 
completely assembled in its package, and measuring the corresponding 
displacements. The force is applied usually by means of a rod inserted 
in a hole cut through the outer container and the cushioning to the packaged 
item. It is convenient to use a low loading rate in the test, and, in doing so, 
the effect of resisting forces that depend on velocity is lost. These forces 
are often of little importance but, in certain designs, it is necessary to con- 
sider them. This is done for velocity damping in Sections 2.5, 2.6, 3.2 and 
3.5. 

Most of Part I is concerned with cushioning having non-linear load-dis- 
placement characteristics. Linear cushioning is rarely encountered, but 
it will be treated first because of its simplicity and because it will be con- 
venient later to express the maximum acceleration in non-linear cases in 
terms of the maximum acceleration in a hypothetical linear case. 


1.2. DERIVATION OF EQUATIONS OF MOTION 


To introduce the method of analysis that will be used in Part I, the sim- 
plest possible system is considered first. The m, system is omitted entirely, 
the mass of the outer container (ms) is neglected, and the cushioning is 
assumed to have no damping or friction. There remain only the mass 
my (the mass of the packaged item alone) and the supporting spring, as 
shown in Fig. 1.2.1. If the spring is linear its displacement is proportional 
to the applied load throughout the range of use (see Fig. 1.4.1). The spring 
rate (k») of a linear spring is a constant usually expressed in terms of pounds 
per inch. The force (P) transmitted through a linear spring is therefore 
given by 

P = koxa, (1.2.1) 
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where x2 is the displacement of m, measured downward from its position a 
first contact of the spring with the floor (see Fig. 1.2.1). For a non-linear 
spring P will be some other function of x: 


P = P(x»). (1.2.2 


To write the equation of motion for the mass m2, we consider the forces 
acting on it at any instant. These are (see Fig. 1.2.2(b)) the spring force 
P and the weight m2g, where g is the acceleration of gravity. When x, 
is positive (i.e., a downward displacement of m2 from its position at first 
contact of the spring with the floor) the spring exerts an upward force P 

















Fig. 1.2.1—Elementary system. 




















pt ida 


t 


Pp 








(a) (b) 
Fig. 1.2.2—Free body diagram for elementary system. 
(a) Spring not in contact with floor. 
(b) Spring in contact with floor. 


on the mass, opposing the weight. The total downward force on mz is 


thus m2.g — P. By the second law of motion, the product of the mass and 
its acceleration at any instant is equal to the applied force: 


mt. = meg — P, (1.2.3) 


where the symbol :, representing the acceleration of m2, stands for the 
second derivative of displacement with respect to time (d’x2/dt). Equation 
(1.2.3) is the law governing the motion of m, as long as the spring is in con- 
tact with the floor. When the spring is not in contact with the floor, it can 
exert no force on the mass so that, in writing the equation of motion that 
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governs before or after contact, the free-body diagram of Fig. 1.2.2(a) should 
be used. Then 


Xo = g. (1.2.4) 


Equation (1.2.4) holds (neglecting air resistance) from the instant the 
package starts to fall until the instant it strikes the floor and from it we 
can find the package velocity at the instant of first contact. Integrating 
(1.2.4) with respect to time, we find 


i = git A, (1.2.5) 


where #2 is the velocity (dx2/dt) and A is a constant of integration whose 
value is found from the initial condition that when ¢ = O (the instant of 
release) 42 = 0. Thus A = Oand 


2 = gt. (1.2.6) 
Integrating again, 
xo= Sg? + B. (1.2.7) 


The value of the integration constant B is found from the initial condition 
that x. = —/h (the height of drop) when¢ = 0. Hence B = —hand 


x2 = 3 gl —h. (1.2.8) 


At the instant of contact, x2 = 0 and, from (1.2.8), the time at first contact 
is given by ¢3 = 2h/g. Substituting this value of ¢ in (1.2.5) we find, for the 
velocity at first contact, 

[42] rom0 =v 2¢gh. (1.2.9) 


We now have the initial conditions for finding the values of the integration 
constants in the solution of equation (1.2.3), which we proceed to obtain. 


dt \ dt 


First multiply both sides of (1.2.3) by dx2/dt and write #2 = d (2) 


x: dXs x: dx 
.—— — = = wt 
My = (2 I n (1.2.10) 


3M - (2) +P 
¢ 


Multiplying by dé and integrating once: 


Lis &3 fe | ; P dx = | mog dx, + C, 
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where C is a constant of integration whose value is determined by the initia 
conditions that #3 = 2gh and x. = Oat the instant of contact. Hence 


C = mogh +f P dx2. 
Substituting the above value of C in (1.2.11), we have 
Ling t2 + [ P dx. = mog(h + x2). (1.2.12 
“0 


It may be observed that (1.2.12) is an energy equation in which the 
terms have the following meanings: 


SMox> is the instantaneous kinetic energy of mz, 
[ P dx, is the energy stored in the spring at any instant. It is also 
Jo 


equal to the area under the load-displacement curve up to 
the displacement 2, 
m2g(h + x2) is the potential energy of the mass at its initial height h + x, 
above the instantaneous position x2. 
Hence (1.2.12) expresses the law of conservation of energy. 
Ordinarily / is very much larger than x2 so that we may write, with good 
accuracy, 


Ling dt. + | P dx. = mogh. (1.2.13) 
0 


To the same approximation, equation (1.2.3) becomes 
Ms Xo -+- P =0. (1.2.14) 


Equation (1.2.14) and its first integral, equation (1.2.13), are convenient 
forms for calculating events at any time during contact. Their use will be 
illustrated in Part II. For calculating only maximum displacement and 
acceleration, the equations become simpler. Let 

W. = weight of the packaged article (= mg), 

d», = maximum displacement of the packaged article, 

G,, = absolute value of maximum acceleration of the packaged article 

in terms of number of times gravity (Gn = | ¥2/g |max), 

P,, = maximum force exerted on packaged article by cushioning. 

We shall limit our study to the practical regions where P > 0 when 
xy > 0. Then it may be seen from (1.2.13) that x2 is a maximum when 
2 is zero, hence 


dm 
[ P dx. = Woh, (1.2.15) 
0 
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and, from (1.2.14), 


Pm 


WW.’ (1.2.16) 


Gn, = 
where P,, is the maximum value of P. If P(x2) is a monotonic function, 
P,, may be obtained from (1.2.2) by substituting d,, for xe: 


Pm = P(dm). (1.2.17) 


In the unusual case where P(x) is not monotonic, the maximum value of P 
in the interval 0 < x2 < d, must be chosen instead of equation (1.2.17). 

The general procedure is to calculate dy from (1.2.15), P» from (1.2.17) 
and then G,, from (1.2.16). If P can be expressed analytically in terms of 2x» 
and if the integral in (1.2.15) can be evaluated in terms of elementary func- 
tions, simple formulas can be found for d,, and G,,. If this is not possible, 
then the integration can be performed graphically or numerically. Both 
of these procedures will be illustrated. In either case the maximum accel- 
eration and displacement are obtained in terms of the weight of the pack- 
aged item, the height of drop and parameters descriptive of the load-dis- 
placement characteristics of the cushioning. 


1.3. Linear ELASTICITY 


For cushioning with a linear load-displacement relation, equation (1.2.1) 
applies. Substituting this value of P in (1.2.15), and performing the in- 


tegration, we find 
2hW's 
ay. cose Salta (1.3.1) 
/ ks 


From (1.3.1) and (1.2.17), 
Py =. 2hW oko, 
and, from (1.3.2) and (1.2.16), 
hk» 
g«~/—. (1.3.3) 
Vw: 
Notice that equation (1.3.3) holds only if there is space available for a 
displacement d,, and if the cushioning is linear and capable of transmitting 
a force P,,. Also, from (1.3.3) and (1.3.1), 
2h 


Gm 


dm (1.3.4) 


5G,  2hW2 
by = W2Gm _ 2AWs 


2h! dm 
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Example: Find the properties of the linear cushioning required so that 
the maximum acceleration will be 50g in a o ft. drop of a 20 lb. article. 
From (1.3.4), 


necessary travel, d, = - rT = 1.44 inches. 
From (1.3.5), 
; 20 x (50)° ; 
spring rate, k, = aS = 694 lbs/in. 
a 30 


From (1.2.16) 
Maximum force P,, = 20 K 50 = 1000 lbs. 
1.4 CUSHIONING WITH NON-LINEAR ELASTICITY 


In practice it is rarely that a packaging system has linear spring charac- 
teristics. Departure from linearity may be due to 


Pp 





= 
Fig. 1.4.1—Linear elasticity. Class A. 


1. Non-linear geometry, such as in the tension spring package described 
in Section 1.7. 

2. Non-linear characteristics of distributed cushioning materials such as 
excelsior and rubber. 

3. Abrupt change of stiffness such as occurs if the packaged item can 
strike the wall of the container. 

For the purpose of developing design formulas it is desirable to have 
analytical functions to represent load-displacement characteristics. It is 
not feasible to have only one family of functions with adjustable parameters 
to fit all possible shapes of load-displacement curves. Therefore, all the 
practical shapes have been divided into six general classes, most of which 
are associated with simple functions having one or two adjustable param- 
eters. The six classes are as follows: 

Class A—Linear Elasticity. This has already been treated. Its load- 
displacement function is 


P= RoXo. (1.4.1) 
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Class B—Cubic Elasticity. 
tom in the anticipated range of use, but the slope of the load-displacement 


This includes cushioning which does not bot- 





function generally increases with increasing displacement as in the curved 
f full line of Fig. 1.4.2. A suitable load-displacement function is 







P = koxe+ rx. (1.4.2) 






| ko is the initial spring rate of the cushioning, as shown by the slope of the 
| dashed straight line in Fig. 1.4.2, and r determines the rate of increase of the 
» spring rate. The same function can be used if the slope of the curve de- 
* creases gradually with increasing load as shown by the curved dashed line 
in Fig. 1.4.2. In this case the parameter r is negative. 










X> 4. % 
Fig. 1.4.3 
Fig. 1.4.2—Cubic elasticity. Class B. 
Fig. 1.4.3—Tangent elasticity. Class C. 


















Class C—Tangent Elasticity. Cushioning that bottoms, but not very 


abruptly, can be represented by the load-displacement function 










2kyod, % 

. < ody 1X2 

P= tan 
Tv 


é . 
2d, 


(1.4.3) 












Referring to Fig. 1.4.3, ko is the initial spring rate and d, is the maximum 
available displacement. ‘The figure shows how the stiffness of the cushion- 
ing (i.e., the slope of the curve) increases as the displacement approaches 
- the maximum available (d,) at hard bottoming. The shape of the curve 
is typical of load-displacement curves for a great variety of packages with 
distributed cushioning. 

Figure 1.4.7 illustrates the wide variety of shapes of non-linear cushioning 
' characteristics that can be obtained with the single function given by equa- 
} — tion (1.4.3) simply by varying the parameter fo; and a similar set is given by 
each value of d,. Although these families of curves do not include all pos- 
sible shapes, one of them can usually be found to fit a practical shape for 
cushioning of this class over the anticipated range of use. 
Class D—Bi-linear Elasticity. This is characterized by a load-displace- 
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ment curve consisting of two straight line segments. 


The load displacen en 
function is (see Fig. 1.4.4) 


P = Ro x2 0 — Xo < ds ) 


(1.4.4 


ll 


P Rp x2 —_ (Ry —_ ky)d, Xo = d,| , 


It is useful especially in situations where very abrupt bottoming is possil)le 
Class E—Hyperbolic Tangent Elasticity. When the mechanism of the 


cushioning is such as to limit the maximum force that can be transmitted 
over a considerable displacement range, the load-displacement function 
ko x2 
P = P tanh P (1.4.5 
0 


is useful. Po is the asymptotic value of the force and &p is the initial spring 
rate (see Fig. 1.4.5). 





Fig. 1.4.4 Fig. 1.4.5 


Fig. 1.4.4—Bi-linear elasticity. Class D. 
Fig. 1.4.5—Hyperbolic tangent elasticity. Class E. 


Class F—Anomalous Elasticity. In occasional instances the load-dis- 
placement curve of the cushioning cannot be matched accurately enough 
by any of the five preceding functions. In such cases a numerical integra- 
tion procedure can be used, as described in Section 1.15. 


1.5 CUSHIONING WitH CuBiIc Exasticity (CLAss B) 


Substituting (1.4.2) in (1.2.15) and performing the integration, we have: 


» 2 ds 
on 40m = Woh. (1.5.1) 
Now, let 
a , 
dy = me (1.5.2) 
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that is, do is the displacement that would take place if the elasticity were 
linear (See equation (1.3.1)) with a constant spring rate ky equal to the initial 
spring rate of the cubic elasticity. Also let 


' 412 hr 
ie 
Then, from (1.5.1), (1.5.2) and (1.5.3) 


dn /2 1 (1.5.4) 
z V BS +vV1i+ B)-. a 


Equation (1.5.4) is plotted in Fig. 1.5.1 which shows graphically how the 


| maximum displacement d,, compares with the ‘equivalent linear displace- 


ment dy’’ as the parameter B is varied. Note that B depends on the weight 
of the packaged item, the height of drop and the shape of the load displace- 
ment curve (as determined by & and r). 


Pp 








Xo 


Fig. 1.4.6—Anomalous elasticity. Class F. 


Similarly we can compare the maximum acceleration G,, with the maxi- 
mum (Go) that would obtain if the load displacement curve were linear with 
spring rate ky. The latter acceleration is given by 


G- Ve (1.5.5) 


and the former is obtained by finding P,, from (1.2.17) and then, from 
(1.2.16), 
Gm 2 ; ‘ 
== B (1+ B(-1+ V1 +4 B). (1.5.6) 


Equation (1.5.6) is plotted in Fig. 1.5.2. 


1.6 PROCEDURE FOR FINDING MAXIMUM ACCELERATION AND DISPLACEMENT 
FOR CUSHIONING WITH CUBIC ELASTICITY 

If the load-displacement curve of a cushioning system has the general 

appearance of Fig. 1.4.2 (where the slope increases or decreases gradually 
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with displacement) the following procedure may be used for estimating 
effectiveness of the cushioning. 


150(0) 





140(0) 


120(0) fp 


110(0) 


100(0) 


dy 


Fig. 1.4.7—Family of oad displacement curves for cushioning with tangent elasticity. 


a. Select the point on the load-displacement curve for which the load 
is equal to the weight of the packaged item multiplied by the allowable 
Gm. Call this load Ps and the corresponding displacement d2. 
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Fig. 1.5.1—Maximum displacement for cushioning with cubic elasticity. See 


equation (1.5.4). 











Fig. 1.5.2—Maximum acceleration for cushioning with cubic elasticity. See 
equation (1.5.6). 
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b. Select another point (d;, P;) about half way toward the origin fro: 
(dz, Po). See Fig. 1.4.2. 
c. Calculate 


dy 


: d3 — di’ 


r (1.6.2 


d. Using the known weight, HW’2, of the packaged item, the specitied height 
of drop /, and ko and r from (1.6.1) and (1.6.2), calculate B, dy and Go from 
(1.5.2), (1.5.3) and (1.5.5). Then calculate the maximum acceleration G,, 
and maximum displacement d,, from (1.5.6) and (1.5.4) or find their values 
from Figs. 1.5.1 and 1.5.2. 

Example: A large vacuum tube, weighing 22.5 Ibs, was packed in a 7” x 
7¢” x 15” carton which was supported on corrugated cardboard spring pads 
ina 103”x113”x18}” carton. The latter was, in turn, packed in 28 pounds 
of excelsior in a 25” x 25” x 30” carton. The tube is rated at 50g and _ the 
package is intended for a drop of three feet. 

A rod was inserted through a hole cut through the three cartons to the 
tube. Load was applied to the rod and the displacement of the tube was 
measured. The data obtained were 


P 


x2 
load in lb displacement in inches) 


100 
200 
300 
400 
500 
600 
700 
800 
900 
1000 
1100 
1200 
1300 
1400 


a 


Cleon ed ede 


Te ne ee 


The data are plotted in Fig. 1.6.1. The resulting curve is suitable for 
classification as either Class B or Class C cushioning. Considering it, for 
the present, as Class B, we take P, = 22.5 X 50 = 1225, and from the curve, 
d. = 1.9 inches. Also, from the curve, take d; = 1 inch and P; = 365 
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lbs. Substituting these values in (1.6.1) and (1.6.2), we find ky = 255 and 
y = 108. Then 


4s hr $x 22.5 & 36 X& 108 


= 28.6. 


-O. 


Ib 

-“- Vw, V 22.: 

Entering Fig. 1.5.2 with B = 5.4 we find G,,/Go = 1.9. Hence 
Gw = 73.6 X19 = S55 


This is close enough to the 50g rating of the tube to call the cushioning 
safe insofar as maximum acceleration is concerned. 


1400+ 

1200 

1000 * 
800} 
600} 
400; 


200} 





0 02 04 06 08 10 12 14 16 18 20 
x 


2 
Fig. 1.6.1—Experimental load-displacement curve for a corrugated cardboard spring pad 
and excelsior cushion. 

The maximum displacement, obtained by entering Fig. 1.5.1 with B = 
5.4 and finding d,,/dyo = 0.75. Then d,, = 0.75 X 2h/Go = 1.95 inches. 
Hence, the package is much larger than necessary since approximately 8 
inches of cushioning thickness is supplied to accommodate 2 inches of 
displacement. 

1.7 THe TENSION-SPRING PACKAGE (CLAss B) 

The tension spring package is useful when the allowable G,, is so small 
and height of drop so great that a large displacement (say d, > several 
inches) is required. The decision as to whether or not a tension spring 
package is indicated may be made on the basis of a preliminary estimate 
of displacement based on the linear case. Suppose the height of drop is 
to be 60 inches and the allowable acceleration for the packaged item is 
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20g. Then, from Equation 1.3.4, the approximate displacement that wl! 
be required is 
2 X 60 ; " 
dn = * = 6 inches. (1.7.1 


The actual maximum displacement in a tension spring package will prove 
to be somewhat more than 6 inches, but the preliminary calculation shows 
the displacement to be large enough to warrant the use of this type of 


cushioning. 
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Fig. 1.7.1—Schematic diagram of a tension spring package. 


A schematic diagram of a typical tension spring package is shown in Fig. 


1.7.1 and a photograph of one design is given in Fig. 1.7.2. The packaged T 
item is suspended on eight identical helical tension springs which diverge by : 
to the outer frame. The analysis and design procedures described in this vert 
and the following section apply equally well if the springs converge from the I 
packaged item to the outer frame. With a slight modification, indicated x 


in the next section, the procedure also applies if four of the springs (say, x 
BJ, DL, EM, OG in Fig. 1.7.1) are omitted. 

In all cases, however, we shall consider only systems having reflected 
symmetry about each of three mutually perpendicular planes through the 
center of gravity of the packaged article. 
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Se 


Fig. 1.7.2—A tension spring package. 


The load-displacement characteristics of the spring system may be found 
by statical considerations. We shall examine, first, the displacement in the 
vertical direction in Fig. 1.7.1, using the following notations: 

P = force applied to the suspended object, 

x2 = displacement of suspended object, 

¥%o = perpendicular distance (JR, Fig. 1.7.1) from inner spring support 

point (J, Fig. 1.7.1) to nearest plane, perpendicular to displacement 
direction and containing four outer spring support points (A, B, C, 


D, Fig. 1.7.1); 
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| 


f; = distance (JA) between spring support points when suspended art cle 
is in equilibrium position, 

f = projection of ¢; on plane ABCD, 

f = ¢€ minus length (between hooks) of unstretched spring, 

k = spring rate of each spring. 


ll 


Consider, first, the action of one pair of springs, say EM and GO of lig 
1.7.1, independent of the remainder of the suspension. Since EM and GO 
lie in parallel vertical planes and the points M and O remain in the initia 
planes of their respective springs during a vertical displacement, the tw; 
springs may be considered to lie in the same plane, and to be translated hori- 
zontally in this plane so that their outer ends are separated by a distance 2/ 
Hence Fig. 1.7.3 may be used to represent the independent action of this 
pair of springs and it is required to find the force Y’ needed to transform Fig. 
1.7.3(a) to Fig. 1.7.3(b). Initially there are two springs, each of length ¢ — ; 





L Bt Z 
E G 
a b 


Fig. 1.7.3-——-Diagram used in discussion of tension spring package. 


and spring constant &, with no initial tension in them. One end of one spring 
is fixed at point & and one end of the other spring is fixed at a point G 
distant 2ffrom £. The springs are then stretched so that the two initially 
free ends are located at a point V equidistant from E and G and distant », 
from line EG. The axis of each spring makes an angle a with EG, where 
, 
; Xe os 
sina = (1.7.2 


VO + x" 
In this state the axial force F in each spring is 
FH=khVF 4a —C4+4) (1.7.3 


and the force 0’, required to equilibrate the two forces F is 2F sina. Con- 
e . e . * U . 
sidering the force (’ as a function of the displacement x. , we write 


, 
2x» 


Q'(%) = : - ie ha me ee} fl (1.7.4) 
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, | | 
O'(s’) = 2k 2’ 
: V1 + 2” 


/ 


r 

Consider, next, the configuration shown in Fig. 1.7.4(a), where one end of 
each of four springs is fixed at a corner of a rectangle of length 2 and width 
2x). Each spring is again of length f— f. The four free ends of the springs 
are drawn together at a common point V at the center of the rectangle (see 
Fig. 1.7.4(b)). The system is in equilibrium in this position. A force ( 
is then applied at Y in the plane of the rectangle and normal to the side 2¢. 


Fig. 1.7.4—Action of springs in a tension spring package. 


The common point X is displaced a distance x2 to V’ (see Fig. 1.7.4(c)). 
Writing s = x2/f, a = xo/¢, we observe that 


Q(z) = O'(s + a) + O(c — a), (1.7.6) 


or, from equation (1.7.5), 


( 3 + a = © ) oe 
2kl< 25—(1—)) _+ ee Chass) 
( V1i-+ (s+ a>) Vi-F (¢ — aF)) 


The standard tension spring package has two sets of four springs so that 
the force P required to displace the common point X a distance x2 is 
P(z) = 2Q(z). (1.7.8) 
If x2 is small in comparison with ¢ (i.e., 2 is small in comparison with 
unity), equation (1.7.8) may be written approximately as 
(1 3 * 
P(s). = 4k€] 22 — — ; (1.7.9) 
Even when x» becomes almost as large as f, equation (1.7.9) has been found, 
experimentally, to be remarkably accurate. 
Writing 


(1.7.10) 











we 
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and 
2\3/2 
— sue ‘ ~ if, (1.7.11 


3 
Pp = K(s + -) (1.7.12 


It is seen, by comparison with (1.4.2) that this is Class B cushioning 
(cubic elasticity). A is the initial spring rate and c determines the rate of 
increase of stiffness with displacement. With the notation kp , r of Section 
1.5, we see that 


equation (1.7.9) becomes 


tro ob 


k=K (1.7.13) 
K aod 
r= Ieee (1.7.14) 


Hence equations (1.5.6) and (1.5.4) may again be used to calculate maximum 
acceleration and displacement. #8 has the same meaning as before (Eq. 
1.5.3). 

To predict the performance, in the vertical direction (Fig. 1.7.1), of an 
existing tension spring package the same procedure as outlined in Section 
1.6 may be used, except that it is not necessary to have a load-displacement 
curve for calculating ky andr. Instead, these parameters may be calculated 
directly from equations (1.7.10), (1.7.11), (1.7.13) and (1.7.14). The 
remainder of the procedure is the same as in Section 1.6(d), 

To predict the performance perpendicular to another face, say AEHD 
of Fig. 1.7.1, it is only necessary, in the calculation of ko and 7, to substitute 
xo for xo, (’ for ¢ (see Fig. 1.7.1) and, in place of bd: 


t , 
’ = 1 — pill — bd). (1.7.15) 
The initial spring rate A for any direction of acceleration may be calcu- 
lated from the initial spring rates K, , Ke, A; in the three directions normal 
to the faces of the frame by using the relation 
1 s? # u 


eo Rtmt x 


(1.7.16) 


where s, ¢, « are the direction cosines of the acceleration direction with 
respect to the normals to the faces of the frame. It is seen, from (1.7.16), 
that the spring rate is given by the radius to the surface of an ellipsoid whose 
principal semi-axes are AK, , Ke, Kg. 
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The displacement direction does not necessarily coincide with the acceler- 
ation direction. The angle 6 between them is given by 


s Cw 
cos 9 = K|] — (4:7:37) 
cos E + K + =| hae 


: where K is defined by equation (1.7.16). 


The spring characteristics may be made the same in all directions and the 


) displacement direction may be made to coincide with the acceleration 
' direction by setting 


5 


3 / Jf £ 
Mm = % = %X% = V/2 
and 
fa fF =" 
(see Fig. 1.7.1). This makes 6 = 6’ = 6”, ¢ = 0.828 and k/K = 0.274 
in the calculations of the next section. 
Re OCEDURE FO SIGNING TENSION SPRING PACKAGES 
1.8 PROCEDURE FOR DEsI I S is 


The design of a tension spring package, as contrasted with the analysis of 
one, must proceed without initial knowledge of values for the parameters ky 
and r, since these cannot be known until the springs are designed. There- 
fore equations (1.5.4) and (1.5.6) cannot be used directly. For design pur- 
poses they are transformed to the following set of formulas: 


/(1 + a’)? 


Vege = (1.8.1) 

L=-—7 =! /( N+ J/0 48) +2 (1.8.2) 

in 2 tt Ve Se TR 82 

din /~ ~ oa 

ee AOS EVES D) (1.8.3) 

uM N W2Gn 2 (1 . B ’ 1 + 14 B 1 84 

——— ot ae i ae + B) 8.4) 
2K OB = 


€ = hn, «De Gu + xX » =) 
( b= 1+ 4/14 7: 1.8.5) 


These formulas have been converted to design curves which are given in 
Figs. 1.8.1 to 1.8.5. The curves are for use in connection with the following 
routine procedure which has been found useful in designing the springs for 
tension spring packages. Reference should be made to Table I. 
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1. Enter, on Line 1, Table I, the weight (IV) in pounds, of the sus- 


w bh 


pended item. This includes the weight of the cradle or other holding 
arrangement and one-third the estimated weight of the springs. 


b=0.3 
1.6 pe er beo.2 


iy b=0.1 


 b=0.0 
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Fig. 1.8.1—Tension spring package design curve. Equation (1.8.1). 


Enter, on Line 2, the height of drop (/) in inches. 


. Enter, on Line 3, the maximum allowable acceleration (G,,) in units 


of ‘number of times gravity.”” This should be determined before- 
hand from tests on the item to be packaged. 
Enter, on Line 4, the dimension x» (inches). 
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€ Sus- 5. Enter, on Line 5, the dimension ¢ (inches). For a package to have 
olding the same spring rate in all directions, f = x9 2 is a necessary con- 
:. dition. 


6. Enter, on Line 6, the value chosen for 6. As 6 becomes greater than 
zero, the stiffness of the whole suspension increases for a given stiff- 
ness of individual springs. The reverse happens for 6 less than zero. 

25 


2.4 bt 


23 


1.6 
1.5 


1.4 


1.2 


1.1 





Fig. 1.8.2—Tension spring package design curve. Equation (1.8.2). 


7. Calculate x/f. 

8. Enter Fig. 1.8.1 with x9/f and find Vc. 

9. Calculate L = h/(VWc (G,,). 

10. Enter Fig. 1.8.2 with Z and find NV. 

11. Calculate A = (W2G,,")/(2hN). This is the initial spring rate of the 
suspension in the direction of x9 . 

12. Calculate f = 3.13 (K W.)'. This is the natural frequency of vibra- 
tion (cycles per second) of the suspension for small amplitudes in the 
xo direction. This should not be close to the natural frequency of 
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vibration of any element in the packaged item, which should be 
determined by test beforehand (see, also, Section 4.2). In any case 
it is advisable to provide damping for the suspension. 


3.4——_,—JJ 


| | | 
3.2+—+ 4 } } + + + + + 
| | 

















Fig. 1.8.3—Tension spring package design curve. Equation (1.7.10). 


13. Enter Fig. 1.8.3 with x»/¢and find k/K. If a four-spring package is 


14. 


desired, instead of an eight-spring package, (see Section 1.7) the 
value of k/K found on Fig. 1.8.3 should be multiplied by two before 
entering it on Line 13 in Table I. This is the only change required 
in the procedure. 

Calculate k = «(£) lhis is the spring rate of each of the springs in 


pounds per inch. 
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15. Calculate B = (2W»h)/(Kce€?). 

16. Enter Fig. 1.8.4(a) or (b) with B and find d,,/(/c ¢). 
17. Calculate dn/f = VWe-dm/Vcl. 

18. Calculate d, = €-d,,/¢. 

19. Calculate (dn/f) + (x0/f). 
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Fig. 1.8.4(a)—Tension spring package design curve. Equations (1.8.3) and (1.8.4 


20. Enter Fig. 1.8.5 with (d,/f) + (x0/f) and find (e/f) — 6. e is the 
stretch of each spring (in inches) when the displacement is d,, inches. 

21. Calculate F,, = k-(e/f)-¢. This is the maximum load (in pounds) 
on each spring. 

22, 23, 24, 25, 26. These are the coil diameter, wire diameter, number of 
turns, fiber stress and length of coils. These quantities are calculated 
from the ordinary formulas, charts or slide rules for helical springs, 
using the values of & and F,, from Lines 14 and 21. 

27. The length inside hooks is entered on Line 27 to group all of the spring 

specifications. 
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Fig. 1.8.4(b)—Tension spring package design curve. Equations (1.8.3) and (1.8.4 
I £ I £ £ j 


1.3 RSE Tnpastan 20. 
6 : aa 3 
11 Z 


09s 








ie) O2 04 06 08 20 12 14 286 12.8 20 


Fig. 1.8.5—Tension spring package design curve. Equation (1.8.5). 
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As an example of the calculations, Table I contains the entries for the 


design of springs for a 21 pound article (including } the estimated spring 


weight) which is to be packaged so as not to exceed 35g in a five-foot drop. 


TABLE I 
Com putation Form for Tension Spring Packages 

1. We (Ibs.). 21 

2. h (ins.). 60 

3. Gu. .s.+ 35 

4. xo (ins.) ; 5 

5. @ (ins.). 7.07 
0. , 0 

7. Cale. xo/f : 0.707 
8. Find Vc from Fig. 1.8.1 0.91 
9, Cale. h/ Ve (Gn = L 0.269 
10. Find V from Fig. 1.8.2 1.265 
11. Cale. W2Gy?2/2hN = K (\bs/in.) 169.0 
12. Cale. f = 3.13 (K/W2)} (cyc./sec.) 8.9 
13. Find &/K from Fig. 1.8.3 0.274 
14. Calc. k = K - k/K (lbs/in.) 46.5 
15. Calc. B = 2W oh/Kef?. 0.368 
16. Find dn/ Wc ¢ from Fig. 1.8.4 0.575 
17. Cale. dm/t = Ve + dn/V cl 0.518 
18. Calc. dm = € + d»/f (ins.) 3.68 
19. Cale. dm/f + xo/€ 1.220 
20. Find e/¢ from Fig. 1.8.5 and line 6 0.580 
21. Cale. Fm = k + €/€ + € (Ibs) 191.0 
22. Coil diameter (ins.) 1.40 
23. Wire diameter (ins. ) 0.207 
24. Number of turns 19 
25. Fiber Stress (Ibs./sq. in.) - 10-8 80 
26. Length of Coils (ins.) 3.93 
27. Length inside hooks (ins.) 7.07 


1.9 CUSHIONING WITH TANGENT ELAsticity (CLAss C) 


This is one of the most frequently encountered classes of cushioning since 
it includes a very common type of bottoming (Figs. 1.4.3 and 1.4.7). The 
load-displacement function (equation (1.4.3)) takes into account the fact 
that the cushioning can be compressed only to a definite amount d, . 

To find formulas for maximum acceleration and displacement, we pro- 
ceed as follows. Substitute equation (1.4.3) in (1.2.15) and perform the 
integration, obtaining 
mete log cos - = —Woeh, (1.9.1) 
ie 2d» 
which may be written as 


rd», rw? Woh 


tan 2d; = V exp( kod ) — i. (1.9.2) 
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Equation (1.9.2) can then be substituted into (1.4.3) to obtain the maximum 
force P», in accordance with (1.2.17): 


dhods . / fe Wah\ 
Pa = —_ , rt ry ee ° oF 
= exp (SU) 1 (1.9.3 





So 
do 
Fig. 1.9.1—Curve for finding maximum acceleration for cushioning with tangent elasticity. 
See equation (1.9.4). 


The maximum acceleration is then obtained from (1.2.16) and may be 


written in the form 
G, 2d, 7, dy 5 
oe xp{=—} —1, 1.9.4) 
Go tily y/o (5 I ( : 


where dp and Gp are defined just as in (1.5.2) and (1.5.5). Go is the maxi- 
mum acceleration that would obtain if the cushioning did not bottom, that 
is, if the spring rate remained constant at its initial value ko. do is the 
maximum displacement that would be reached under the same linear con- 
ditions. Hence G,,/Gp is a multiplying factor to be applied to a hypothetical 
linear cushioning to take into account the effect of bottoming. The multi- 
plying factor depends only on the ratio (d,/do) of the amount of space 
actually available to the amount of space that would be used under linear 
conditions. 
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he ratio G,,/Go is plotted against the ratio d,/dp in Fig. 1.9.1. It may 
be seen that the multiplying factor increases very rapidly as the displace- 
ment ratio (d,/do) falls below unity. For example, if the cushioning, with 


tangent elasticity, reaches hard bottoming (d,) when only 80°% of the 
required displacement (do) is attained, the acceleration is multiplied by 
3.5; if only 60% of the required displacement is available, the acceleration 
is multiplied by 11.5. 

Example: To illustrate with a numerical example, consider the case already 
discussed in Section 1.3, where we found that a spring rate of 694 |bs/in 
and a displacement of 1.44 inches were required to limit a 20-pound article 
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Fig. 1.9.2—Curve for finding maximum displacement for cushioning with tangent elasticity. 
See equation (1.9.5). 
to an acceleration of 50g in a 3-foot drop with linear cushioning. Let us 
suppose that only 1.15 inches are available, instead of 1.44 inches, and that 
the cushioning has tangent elasticity starting with a spring rate of 694 
lbs./in. Entering the curve of Fig. 1.9.1 at d,/do = 1.15/1.44 we find 
G»,/Go = 3.5. Hence the maximum acceleration will be 175g instead of the 
required 50g. This illustrates the wide variations in acceleration that may 
occur as a result of minor variations of dimensions in high G packages. 
It is not necessarily true that the 175g test is 3.5 times as severe as the 
50g test for all elements of the supported structure, since the severity de- 
pends also on the shape and scale of the acceleration-time relation. The 
factor may be more or less than 3.5 but it will be very close to this value for 
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all high-frequency elements of the structure. This subject is treated jp 
detail in Parts II and III. 

The maximum displacement d,, , in the case of tangent elasticity, may be 
calculated from equation (1.9.2) or, in terms of d,/do , from 


a ee w (dV , 
i. 2 exp 8 (2) | ‘iia 


The ratio d,,/d» is plotted against d,/do in Fig. 1.9.2. 

The use of Fig. 1.9.2 can be illustrated with the example already calcu- 
lated, in which d,/dy = 1.15/1.44 = 0.8. Entering the abscissa of Fig, 
1.9.2 with d,/do = 0.8 we tind d,,/d, = 0.915. Hence the maximum dis- 
placement will be 0.915 & 1.15 = 1.05 inches. 


1.10 Opriwum SHAPE OF LOAD-DISPLACEMENT CURVE FOR TANGENT 
ELASTICITY 

It is possible to choose the best shape for the load-displacement curve 
of the cushioning from those represented in Fig. 1.4.7. This will be, of 
course, not the best of all possible curves, but only the best among “tangent 
elasticity” curves. The best shape is detined as the one that yields the 
smallest maximum acceleration (G,,) for a given weight (IW), height of drop 
(4) and available space (d,). This leaves the initial spring rate (ko) as the 
only remaining variable. ‘To find its optimum value (say ka), set equal to 
zero the derivative of G,, (equation (1.9.4)) with respect to ky , remembering 
that Go and do are functions of ky. The result is 


(_eWeh\ (eal) | _ ae 
~ hn PL ae * on 


; 3.1W oh 
ga 

Substituting (1.10.2) in (1.9.4) we find the minimum value (G;,,) of maxi- 
mum acceleration to be 


from which 


> 


(1.10.2) 


4 3.9h 


ian = 
ag dy, 


(1.10.3) 


To illustrate the application of equations (1.10.2) and (1.10.3), consider 
again the case of the 20-pound article dropped from a height of three feet. 
We found that a linear spring, with a spring constant of 694 Ibs/in, would 
limit the maximum acceleration to 50g if 1.44 inches of displacement were 
available. If only 1.15 inches of displacement are available, and the initial 
spring rate is kept at 694 lbs/in, we found the maximum acceleration to be 
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175g if the cushioning bottoms with tangent elasticity. Now, according to 
equation (1.10.2), the best initial spring rate for cushioning with tangent 
elasticity would be 

31 MIO X SB 


ko > \e 
(1.15)? 


= 1690 lbs/in. 


In this case, equation (1.10.3) gives, for the maximum acceleration, 


G’. ia 3.9 X 36 = 1220. 

1.15 rm 
Hence, confronted with a space limitation less than that required for a 50g 
linear spring, it is better to use an initial spring rate higher than that for the 
50g linear spring in order to strike an economical balance between displace- 
ment and bottoming. The best balance, among cushionings having tangent 

elasticity, is obtained by using equation (1.10.2). 

If no factor of safety is considered, it would be still better not to use a 
bottoming type of cushion at all. From equations (1.3.5) and (1.3.3) 
it can be seen that a linear spring with a constant of 1090 Ibs/in will give 
only 63g with a displacement of 1.15 inches. Such a spring, though, would 
bottom very sharply at a drop slightly higher than 3 ft. and would give 
an acceleration much greater than cushioning with tangent elasticity which 
bottoms more gradually. This may be important if there are high-fre- 


quency, brittle elements in the packaged article (see Part III). 


1.11. PROCEDURE FOR FINDING MAXIMUM ACCELERATION AND 
DISPLACEMENT FOR CUSHIONING WITH TANGENT ELASTICITY 
(CLAss C) 


To illustrate the use of the equations and curves for Class C cushioning, 
the same example used for Class B will be used, as it was observed that the 
experimental load-displacement curve in that example (Fig. 1.6.1) is a 
border line one which can be treated as either B or C. 

By laying a straight edge along the first part of the curve (Fig. 1.6.1), 
the average initial spring rate is found to be 305 Ibs/in. This value is taken 
as ko in the present case. 

The next step is to find a value of d, such that a graph of P/d, vs X2/d> 
will fall slightly above the curve ky = 30(0) in Fig. 1.4.7; d, must be greater 
than 2 inches, since that displacement was obtained in the experiment. 
As a trial take d, = 2.25 inches and test it at one point, say the experi- 
mental point P = 300 lbs., x. = § in. Then P/d, = 133 and x2/d, = 
0.39. The point (0.39, 133) falls below the curve ko = 30(0) in Fig. 1.4.7. 
Next try d, = 2.5 inches. In this case, for the experimental point P = 














| 
| 
/ 
} 
| 





386 BELL SYSTEM TECHNICAL JOURNAL 


300, v2 = §, we find P/d, = 120, x2/d, = .35. This falls slightly above th: 
ky = 30(0) curve as required. The whole experimental curve is the Th 
plotted to the coordinates P/2.5 vs x2/2.5 and is found to fit as closely a 
necessary. Hence the parameters are adopted as ky = 305, d, = 2.5. 
We can now calculate the maximum acceleration that the tube wi 
receive in, say, a three-foot drop test. First calculate, from equations (1 
and (2), 


/y # ? > ? — 
ko 305 
/2hk ‘2 5 
Go = / 2hko = ys x 36 x 3 : = 31.3. 
Vw.” 4 22.5 
Then d,/dy = 1.08. Entering Fig. 4 with this value we find G,,/Gp = 1.82. 


Hence the maximum acceleration is: 
Gm = 31.3 K 1.82 = 57g. 

Finally, entering Fig. 5 with d,/do = 1.08 we find d,,/d, = 0.8. Hence the 
maximum displacement is d,, = 0.8 * 2.5 = 2.0 inches. This indicates 
that the load-displacement test was carried far enough to cover the range 
up to a three-foot drop. 

It may be observed that the results obtained, by treating the same data 
as Class B or Class C cushioning, agree within a few per cent. This is ri 
because, in the example chosen, both B and C curves can be made to fit the 
experimental load-displacement curve. 


1.12 CONSEQUENCES OF ABRUPT BoTTromMING (CLAss D) : 
It is useful to examine cushioning systems that can bottom more abruptly 
than Class C cushioning. Abrupt bottoming is possible, for example, in a 
tension spring package lacking a snubbing device. An estimate of the 
increase in acceleration can be made by studying the case of bilinear elasticity . 
(Fig. 1.4.4). Here we have a spring rate ky up toa displacement d,, follow- ‘ 


ing which the cushioning has a different spring rate ky. ko represents the 
average spring rate before bottoming and &, can represent the much greater 
stiffness of the wall of the container. 


If do > ds, that is, if 
Wo 
,/™ allt A (1.12.1) 
ho ) 





the suspended article will bottom and the maximum displacement and ( 
acceleration are obtained by using both of the equations (1.4.4) in evaluating ( 
the integral in (1.2.15). Thus, 


j lin 
| Ry xX dxXo + | [kino — (Ry — Ro)d,| dxo = Weh. (¥:32:2) 
0 


“ds 
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The remainder of the procedure for finding G,, is the same as before. 
The value of d,, found from (1.12.2) is substituted for x. in the second of 


O|P 


© O01 02 03 04 05 06 07 08 09 10 
ds 
co 


Fig. 1.12.1—Curves for finding maximum acceleration as a result of abrupt bottoming. 
See equation (1.12.3). 


(1.4.4) and the value of P,, , thus obtained, is used, in (1.2.16), with the 
result: 


/k d ky 
Tm = Gp : : ——}), (1.12.3) 
Gu = Go + E(t) L124 


where Go is the acceleration that would be reached if a displacement d, 
were available: 
/ 
iy afte. (1.12.4) 
V WwW, 

The ratio G,,/Go is plotted against d,/dp in Fig. 1.12.1 for several values 
of k,/ko. Since, in practice, k, might be thousands of times as great as hy . 
it may be seen that the increase in maximum acceleration can be very large 
even when d, is only slightly less than dy). It is apparent that a snubbing 
device is desirable in a tension spring suspension. This is especially true 
when considering high-frequency elements of the packaged article. It will 
be shown, in Part III, that low-frequency elements are not affected as much 
as might be expected from consideration of maximum acceleration alone. 
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1.13) CUSHIONING WitH HypEeRBoLIC TANGENT ELasticity (CLass E 

In the preceding sections, there have been considered four types of elas- 
ticity (linear, cubic, tangent and bilinear) that fit the load-displacement 
characteristics of the more common cushioning materials and devices 
There now remains the problem of finding more nearly ideal shapes of 
elasticity. By ‘tmore nearly ideal” is meant a shape which will result in a 
smaller maximum displacement for a given maximum acceleration. This is 
important in the packaging of very delicate articles if shipping space is 
limited. 

It may be observed (equation (1.2.15)) that the total area under the 
load-displacement curve is equal to the maximum energy of the system. 
The maximum ordinate of the enclosed area is proportional to the maxi- 
mum acceleration. Hence, if we wish to (1) limit the maximum acceleration 
(2) accomodate a given kinetic energy and (3) have as small a displace- 
ment as possible, the best shape for the load displacement function is P = 
constant, where the constant is the product of the supported mass and the 
maximum allowable acceleration. 

It is not practical to obtain this ideal shape exactly, for there will always 
be a finite initial spring rate and a rounding off of the load-displacement 
curve to the limiting maximum load. A function which represents this 
practical condition (and also includes the ideal case) is the hyperbolic 
tangent function mentioned in Section 1.4: 


ko Xo 


P = Po tanh P, 


(1.13.1) 


The formulas for maximum acceleration and displacement are found in 
the same way as for the other classes of cushioning with the results: 


Po = Wohkp 
dm = cosh exp { 53 (1.13.2) 
ko 0 
or 
ae do Po —] : WiG: 
dn = WG, cosh exp ( 2P2 ) (1.13.3) 
and 
Po Ro dm 
im = => tanh - 13.4 
C iW. anh P, (1.13.4) 
or 


P| en(—) ass 


Gn = 


wile 


W 


P 
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where, as before 


dy = Go = 
ky” 
Equations (1.13.3) and (1.13.5) are plotted, in Figs. 1.13.1 and 1.13.2, 
against the dimensionless parameter Po/W Go. The latter is the ratio of 
the maximum force, that the hyperbolic tangent cushioning will transmit, 





1.5 


W2Go 
Fig. 1.13.1—Maximum displacement for cushioning with hyperbolic tangent elasticity. 
See equation (1.13.3). 
to the force that linear cushioning would transmit under the conditions 
specified. 

To find the value of ko which yields the minimum value of acceleration 
for a given maximum displacement, differentiate (1.13.4) with respect to ko 
and set the result equal to zero: 

2 Rod, 
sech 


= (), (1.13.6) 
Po 


This is satisfied by ko — ~, which represents the rectangular load dis- 
placement curve and confirms the conclusion reached from energy 
considerations. 
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Taking the limit of (1.13.4) as kj — *, we find the optimum acceleration 
to be 


./ 0 = 
Gn = We (1.13.7 
The corresponding maximum displacement is found, from (1.13.2) to be 
Pie Weh " he (1.13.8) 
- P, Gi. Ris 





Fig. 1.13.2—Maximum acceleration for cushioning with hyperbolic tangent elasticity. 
See equation (1.13.5). 


1.14 Mixrum Space REQUIREMENTS FOR VARIOUS CLASSES OF 
CUSHIONING 
It is interesting to compare the minimum amount of space for displace- 
ment that can be attained with the various kinds of cushioning that have 
been discussed. 


, h 

Hyperbolic Tangent Elasticity dn = G’ 
, ae 2h 
Linear Elasticity dm = G,, 
hr . ——_— , 3.9h 
Tangent Elasticity dn = —) 
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Cubic elasticity will give a d,, somewhat more or less than 24/G,, depending 
upon whether the parameter r is positive or negative. 

It is seen that a factor of almost four can be gained, in the linear dimensions 
of the cushioning space required, by replacing the tangent type of cushioning 
with the hyperbolic tangent type. 

There are several ways of obtaining a load-displacement curve with a 
shape similar to the hyperbolic tangent curve. One of the most interesting 
is suggested by the fact that the load-displacement curve of a strut has 
approximately this shape. Hence a bristle brush has the proper 
characteristics. 


TABLE II 


_ 0 0 0 
0.500 .500 120 30 6. 
0.100 600 150 9<s 43. 8 
0.100 700 205 i 61. ay. 
0.100 . 800 290 24.8 86. 15 
0.100 900 410 3S. 121. 22 
0.100 000 585 8 170. 31 
0.050 050 730 32. 203. 39 
0.050 100 950 2. 245.3 51. 
0.050 .150 1370 58 303.5 74. 
0.025 175 1680 38. 341. 91 
0.025 . 200 2240 390. 21, 121. 
0.0125 .2125 | 2620 30. 421.. 22 141.: 
0.0125 .225 3200 36.4 457.7 24. 173. 
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1.15 NuMERICAL MetHop FOR ANALYZING CLASs F CUSHIONING 

The numerical method to be described is one that has been adapted from 
a graphical one used by the Committee on Packing and Handling of Radio 
Valves of the British Radio Board. The method has advantages of sim- 
plicity in concept and ease of application, especially when the load-displace- 
ment curve of the cushioning does not resemble closely one of the Classes A 
to E described above. It has the disadvantage that it does not yield, 
directly, numerical factors by which the spring rate or depth of cushioning 
should be changed in the event that the analysis reveals inadequate or more 
than adequate protection. 

The method is based on the fact that the area under the load-displace- 
ment curve of the cushioning represents the energy stored in, or absorbed 
by, the cushion. The total amount of energy that must be transferred is 
equal to the product of the weight (W,) of the suspended item and the 
height (4) of drop. By finding the abscissa (22) and its ordinate (P) which 
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include an area W’yh, the maximum displacement is immediately x2 and the 
maximum acceleration is the quotient P/W, , in accordance with equations 
(1.2.15) and (1.2.16). 

As actually applied in the present instance, the British method was 
modified slightly to make the procedure a routine numerical one. The 
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Fig. 1.15.1—Experimental load-displacement curve for Table II. 


computing form is given in detail in Table II, in which the data are taken 

from an experimental load-displacement curve (Fig. 1.15.1) for the end spring 

pads of a vacuum tube package. The load-displacement data are listed in 

Columns 3 and 4 of Table II. The meaning of each column jn the table is as 

follows. 

Column 1. n(= 1, 2, 3---) is the number that identifies the displacement 
(and corresponding load) up to which the area under the load- 
displacement curve is to be calculated. 
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Column 4. 
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. A(x), is the increment of displacement between (1x2),; and 


(t2)n- A(ve)n = (Xe)n — (X2),-1, See Fig. 1.15.2. Note that, 
as the curve becomes steeper, A(xv2), is taken smaller for better 
accuracy. 


(v2), is the displacement associated with the rn" point (see 
1.15.2). 


P,, is the load that produces displacement (.r2),, . 


} 


men 





OA, | dA,| | | oA, 
(XQ), XD), (XD), 2), 
‘ beet (XIng 
A(X), 4(XQ), AIX) SA0G), 





Xp = DISPLACEMENT 


Fig. 1.15.2—Graphical illustration of numerical method of calculating area under load 


Column 5. 


Column 6. 


Column 7. 


Column 8. 


displacement curve. See Table IT. 


AA, = $A(x2)n(Pr-1 + Pn) is the area of the trapezoid with 
altitude A(x,), and bases P,; and P,,. It is approximately the 
energy absorbed by the cushioning in displacing from (x2), ; 
to (X2)n. 

A,, is the sum of all the trapezoidal areas from x2 = 0 to x2 = 
(x2), . It is approximately the total energy the cushioning can 
absorb in displacing an amount (x2), beginning at zero dis- 
placement. Note that Ao is always equal to zero. 

h, = A,/W2 is the height of fall that will cause the cushioning 
to displace an amount (x2), . In Table II, HW’, = 18.5 pounds. 

G, = P,/W, is the maximum acceleration experienced by the 
suspended mass when dropped from a height /,, . 
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Fig. 1.15.3—Maximum acceleration vs. height of drop for an 18.5 pound article supported 


HEIGHT OF DROP (INCHES) 


on cushioning with the load-displacement curve of Fig. 1.15.1. See Table IT. 


From the last two columns of the table a curve of height of drop vs. the 


corresponding acceleration may be plotted as in Fig. 1.15.3. 


PART II 


ACCELERATION-TIME RELATIONS 


In Part I we were concerned primarily with the maximum acceleration of 
In this part we shall study the details of the variation 


the packaged item. 


of acceleration with time in order to have this information available for our 


2.1. INTRODUCTION 
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study, in Part III, of its influence on the response of elements of the packaged 
item. 

The first case to be considered will be the simple single mass and linear 
spring example described in Sections 1.2 and 1.3. Following this the 
phenomenon of rebound of the package will be considered. The influence 
of velocity damping and dry friction will be studied; and, finally, the effects 
of non-linearity of the cushion elasticity on the acceleration-time relation 
will be investigated. 


2.2 ACCELERATION-TIME RELATION FOR LINEAR ELASTICITY 


Returning to the elementary example studied in Sections 1.2 and 1.3, we 
first write the equation of motion for the mass m:, on its linear spring of 
spring rate ke (see Fig. 1.2.1.). Equation (1.2.3) becomes 


Mot, + kext. = mog. 
Using the initial conditions 
[Xe]imo0 = 0, 
[te]no = V/2gh, 
the solution of (2.2.1) is 


sin (wot — a) + g 
Qs 


Vet Dai gh 


w} 


x2. = 


WT 4 Mah Get ~ 0) + . 


k3 


Qn 


We . 
ky dm, 


—j 
= tan’ — 


Bis, oad 
we V 2gh 


tan 


we is the circular frequency, f2 is the frequency and 7» is the period of vibra- 
tion of the mass mz on its spring; d,, has the same definition as in Section 1.3 
(equation (1.3.1)). 

Now, W2/ke is the static displacement of the mass m2 on its spring. This 
is usually very small in comparison with the maximum displacement (d,,) 
during impact. Hence W2/k, will be neglected, and (2.2.5) becomes 


x4 = 7. i sin wot. (2.2.8) 
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Differentiating (2.2.8) twice with respect to /, we find, for the acceleration 
Be = —wadm SIN wel = —w2 2gh sin wot. (2.2.9) 


Hence the absolute magnitude of the maximum acceleration is 


Xe max > ben 2h Ry 
G, = ims = i = = (2.2.10) 
g g Ws 


as before. 


Gm 
¥ 


¥ — 


Fig. 2.2.1 Fig. 2.2.2 











Fig. 2.2.1—Half-sine-wave pulse acceleration. See equation (2.2.9). 
Fig. 2.2.2—Oscillogram of a half-sine-wave pulse obtained with a piezo-crystal 
accelerometer. 

Equation (2.2.9) shows that the acceleration varies sinusoidally with time. 
It rises from its initial zero value to its maximum in a time 7/2w» , at which 
time the displacement also reaches its maximum value. The acceleration 
returns to zero again at time 7/w.. At this time the displacement is also 
zero. This is the end of the range of applicability of equation (2.2.9); for 
when ¢ becomes slightly greater than m/w. , a tension in the spring is required. 
Since no mechanism, such as a large mass ms, (Fig. 0.2.2), has been supplied, 
to allow a tension in the spring to develop, the system will rebound from the 
floor at the end of the half period r/w.. The acceleration is thus a half- 
sinusoidal pulse of duration m/w. = 72/2 and amplitude G,,g as illustrated 
in Fig. 2.2.1. An oscillogram of such a pulse obtained with a piezo-crystal 
accelerometer is shown in Fig. 2.2.2. 


2.3 PACKAGE REBOUND. 


The presence of the mass of an outer container will affect the acceleration 
after the first half cycle of displacement. The outer container is represented 
by the mass m; in the general idealized system illustrated in Fig. 0.2.2 and in 
the simpler system (Fig. 2.3.1) that we shall consider now. 
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Two pairs of equations are necessary to describe the action of the system; 


one pair applies during the time of contact of m; with the floor and the 
second pair applies if rebound occurs. 

The mass m3 is assumed to be inelastic (see Section 0.2) so that, during the 
interval of its contact with the floor, the equation of motion for mz will be 
the same as before (2.2.1). In addition, there will be an equation of equili- 
brium for the mass ms : 


R = koxo + mag ; 8) 


where R is the upward force exerted by the floor on ms . 














MA 
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3 
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X3 








Fig. 2.3.1 


Two-mass system representing packaged 4x ticle, linear cushioning and 
outer container. 


Equations (2.3.1) and (2.2.1) will hold as long as Ris positive. 


To find out 
when R > 0, solve (2.2.1) for k2x2 and substitute in (2.3.1): 


R = W, + Ws — MoXo ° €2:3 2) 


That is, a necessary condition for rebound is that the mass of the cushioned 
article, multiplied by its maximum acceleration, exceeds the total weight 
of the package. The condition for rebound may be written 


Ga . Bs n ss (2.3.3) 


This is a necessary, but not a sufficient, condition for rebound because there 
will be energy losses as a result of damping and permanent deformation. 
G» will generally have to be considerably greater than the right hand side 
of (2.3.3) for rebound to occur. 

If rebound does not occur, equation (2.2.9) 


continues to apply, except for 
damping which will be considered in Section 2.5. 





| 
| 
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2.4 Motion AFTER REBOUND 


If rebound occurs, the equations of motion for the two masses, m2 and m3 . 
are 


Mok, + ko(x2 — X3) = mog, (2.4.1) 
M3i3 — ko(x2 — x3) = meyg. (2.4.2) 
Multiplying (2.4.1) by m3 and (2.4.2) by mz and subtracting, we find 
my + key = 0, (2.4.3) 
where 
y = %2— %, (2.4.4) 
m = ae (2.4.5) 





Fig. 2.4.1—Oscillogram illustrating the half-sine pulse followed by the higher frequency, 
lower amplitude vibration of the packaged article in a rebounding package. 


Equation (2.4.3) is the equation governing the vibration of the two-mass 
system asa simple oscillator. The circular frequency of the vibration is 


o = y/" (2.4.6) 
m 


and it may be noticed that this frequency is always greater than we (equation 
(2.2.6)). This fact is important in estimating the effect of vibrations on 
elements of the packaged item (Section 3.5). 

w is also the frequency of vibration of the packaged article during the 
interval of free fall. This vibration (usually of small amplitude) is initiated 
by the sudden release of the dead load displacement of the packaged article. 

As an intermediate step in obtaining the acceleration after rebound we 
shall find the magnitude of the relative displacement (y) of the two masses. 
To do this it is necessary to solve equation (2.4.3) with the appropriate 
boundary conditions. Calling ¢, the time at which ms; leaves the floor, we 
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must find y and y at ¢ = ¢,. Since ms; is motionless at ¢ = ¢, , the relative 
displacement and velocity at that time are identical with x. and #2 respec- 
tively. The former is simply the stretch of the spring necessary to just 
pull the mass ms off the floor, i.e., 
_Ws 
To find the velocity at ¢ = ¢, , substitute (2.4.7) in (2.2.4) and also substi- 
tute ¢, for? in the latter. This gives an equation fordetermining/,. Then, 
returning to (2.2.4), differentiate it once to obtain %2 and substitute for ¢ the 
value ¢, just found. The result is 


[ylene, = ([xe]enr, = (2.4.7) 


an gW3(2We + W3) 
Weke F 


The solution of (2.4.3) with initial conditions (2.4.7) and (2.4.8) is 


y= od 4/2" ae sin (wt — ¢), (2.4.9) 
Ww ko 


1 alylene, 

[Vlemt, 

We are now in a position to find the acceleration of the packaged item 
after rebound. Substitute y of (2.4.9) for x2 — x, in (2.4.1) to obtain 


(2.4.8) 


[Xo]eme, = [vent =— V/ 2gh 


where ¢ = wf, — tan 


a 
- 
9 


f=g+™ 


w 


2gh — 188 sin (wt — ¢). @ (2.4.10) 


To obtain a simple formula for the ratio of the maximum accelerations 
after and before rebound, let us assume that both are much greater than 
gravitational acceleration. Then if 


G, = maximum number of g’s after rebound, (maximum of (2.4.10)) 
Gn = : maximum number of g’s before rebound, 


we find, from (2.4.10), neglecting the term g outside the radical, 


ee /! me (2.4.11) 
Gn W2 othe W3 2hke 
Hence, the maximum acceleration after rebound is always less than the 
maximum acceleration before rebound. Therefore, conditions after re- 
bound need only be examined when the frequency after rebound (see equa- 
tion (2.4.6)) is near the natural frequency of vibration of a critical element 
of the packaged item (see Section 3.5). 
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The complete acceleration history of a rebounding package with un- 
damped linear cushioning is thus a half sine wave pulse of amplitude G,, = 
V 2hke We and duration +/w. followed by an oscillating acceleration of 
amplitude given by (2.4.11) and frequency given by (2.4.6). Such a wave 
shape is shown in Fig. 2.4.1. 


2.5 INFLUENCE OF DAMPING ON ACCELERATION 

The presence of damping in cushioning is always desirable to prevent the 
building up of large amplitudes as a result of periodic disturbances. How- 
ever, damping also has an effect on the maximum acceleration that is at- 
tained in a drop test. From the latter point of view there is an optimum 
amount of damping and an amount that should not be exceeded if the maxi- 
mum undamped acceleration is not to be exceeded. 

We shall consider the case of a linear cushion with damping proportional 
to velocity. The system is represented in Fig. 2.5.1. With the addition 


ms 
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Fig. 2.5.1—Idealization of linear cushioning with velocity damping. 


of the damping term the equation of motion of my» , during contact of the 
package with the floor, is 


MoXo + CoXe + Rox = QO, (2.5.1) 


in which c. is the damping coefficient of the cushioning. Equation (2.5.1) 
is more conveniently expressed as 


a + Bowers + were = 0, (2.5.2) 
where 
ko S 
wo = /*, (2.5.3) 
Mo 
a= — 2.5.4) 
o = * 5. 
“ —— Q2mow ( 


we is the undamped circular frequency of vibration of mz, on its spring and 2 
is the fraction of critical damping. 62: = 0 means no damping and ~2 = 1 
means just enough damping so that there will be no oscillation if the pack- 
aged article is displaced and released. 
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The acceleration solution of (2.5.2), with the initial conditions of the drop 
test (see (2.2.2) and (2.2.3)) is 


‘ We 2eh Zawe 
n= Vv +8 e §?*2* cos (wolv/1 — B82 + y) (2.5.5) 
Vi— B : 


where 














Fig. 2.5.2—Acceleration-time curves for linear cushioning with various amounts of 
damping (no rebound). See equation (2.5.5 
The acceleration is thus a damped sinusoid with an abruptly reached 
initial value whose magnitude depends upon the amount of damping. For 
small damping, the initial acceleration is small and then the acceleration 
increases, but never reaches the value that would be reached without any 
damping. For high damping (8: > 0.5) the initial value is greater than 
without any damping and falls off thereafter. Figure 2.5.2 shows the shapes 
of the acceleration time curves for several values of 82. All of the curves 
are for no rebound. It may be seen, from equation (2.5.5) and Fig. 2.5.2 
that the addition of damping changes the shape of the acceleration-time 
relation in three ways. First, a damped sinusoid replaces the pure sinusoid ; 
second, the frequency is reduced; and, third, the initial phase is changed. 
It is useful to consider in detail the effect of damping on maximum accel- 
eration. Let 
Gm = Maximum number of g’s with damping 


‘Qhke ; ee 
= W. > maximum number of g’s without damping. 
2 


Go 
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Fig. 2.5.3—Influence of velocity damping on maximum acceleration. See equation (2.5.5). d 
Then, from (2.5.5), at 4 = 0 e 
Gm i 
~ = 2p, (2.5.7) 4 
Go ; 
and, after ¢ = 0, | 
| G lex ; : 
a eg Pavttn (2.5.8) 4 
| Go ‘ 
) where ¢,, , the time at which the maximum occurs, is given by ; 
» _ (l — 462)-V71 — B? a 
tan wlm V1 — 63 = A Vv ~ Ba (2.5.9) | Fig 
- B2(3 — 463) 
The largest value of G,,/Go from (2.5.7) and (2.5.8) is plotted against 82 ; 
in Fig. 2.5.3. It is shown there that, as the damping is increased from zero, f gt 
the maximum acceleration first decreases to a minimum of 80% of Gp and } cc 
then increases to Go at 50% of critical damping. In this interval the maxi- 4 
mum acceleration occurs after ¢ = 0. For damping greater than 8. = oa 
0.5 the maximum acceleration occurs at the instant of contact and increases Ww 
in direct proportion to Be . is 
, is 
2.6 INFLUENCE OF DAMPING ON REBOUND , 
In considering rebound without damping, it was found that rebound does ‘ 


not occur unless the product of the maximum acceleration and the sus- 
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pended mass exceeds the total weight of the package. It was not necessary 
to distinguish between maximum acceleration on the first downstroke and 
first upstroke, since these are the same when there is no damping. With 
damping, however, the maximum acceleration on the first downstroke is 
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Fig. 2.6.1—Influence of velocity damping on maximum upstroke: acceleration. See 
equation (2.5.5). 


greater than that on the first upstroke (Fig. 2.5.2) and it is the latter that 
controls rebound. Hence damping inhibits rebound. 

For example, with 50% of critical damping (82 = 0.5), equations (2.5.8) 
and (2.5.9) and Fig. 2.5.2 show that for the first downstroke G,,/Go = 1 
while for the first upstroke G,,/Go = 0.164. Hence the tendency to rebound 
is reduced by a factor of six when damping to the extent of 50% of critical 
is added to an undamped package. 

The ratio of the maximum acceleration on the first upstroke to the maxi- 
mum undamped acceleration is plotted in Fig. 2.6.1 for various values of 8: . 
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2.7 INFLUENCE oF Dry FRictTION ON ACCELERATION AND DISPLACEMEN) 

By “dry friction” is meant friction that is independent of velocity excep) 
for sign. During contact of the package with the floor the motion of m 
might be opposed by a constant friction force F. Such a force is developed, 


for example, in a package with corrugated spring pad cushioning by rubbing 


against the side and end pads in a top or bottom drop. A typical idealized 


ic 
te 


a4 Xo 





Tr 
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a 


Fig. 2.7.1—Load vs. displacement for cushioning with dry friction. 





!oad-displacement curve is shown in Fig. 2.7.1. For the first downstroke 
of mz, the equation of motion of mz is 


MoXo + Rov. = —Ff, (2.7.4) 
With initial conditions 


[ve]--0 = 0, [¥2}e9 = V 2eh (2.7.2) 


the solution of (2.7.1) is 


XX. = / oe (; ) sin (wot +a) — ; (2.7.3) 
2 2 
2W 2 
dy = 4/7, 


ee 
ko dy We Go 4 


where 
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2hke 


Go = V We ° 
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d) and Go are the maximum displacement and acceleration that would obtain 
if no friction were present. From (2.7.2) the maximum displacement with 


friction is 


dn = 4 i + (; ) = ; : (2.7.4) 


Hence, the presence of friction decreases the maximum displacement since 
dm < do . 
From (2.7.3) the acceleration is 


: = — y/ «i 4 (7) sin (wet + a), (2.7.5) 


so that the maximum acceleration is 


Gn = 4/ G2 + (F,) (2.7.6) 


which is greater than the maximum acceleration without friction. 

It would appear, at first glance, that cushioning with friction always 
gives a greater acceleration than the corresponding cushioning without 
friction. However, the reverse is actually true provided we allow the same 
displacement in both cases. This may be done, as may be seen from (2.7.4), 
by decreasing the spring rate in the cushioning with friction to 

2F 


k = ay ° 
r = kp d, ( 


i) 
~~ 
~J 


The maximum acceleration in the cushioning with friction is then, from 

(2.7.6), 

F F —G 
— <= —, 


ame (2.7.8) 
W. , 


Gr = G — 


That is, for the same maximum displacement, the maximum acceleration 
is reduced by the addition of dry friction. 


2.8 ACCELERATION-TIME RELATION FOR CUBIC ELASTICITY 
As an example of the effect of nonlinearity of the cushioning on the shape 
of the acceleration-time function, the case of cubic elasticity (Class B) 
will be analyzed. The system to be considered is illustrated in Fig. 1.2.1, 
and the load-displacement relation for the cushioning is given by 


P = Roxe + rx. (2.8.1) 








406 BELL SYSTEM TECHNICAL JOURNAL 


Substituting (2.8.1) in (1.2.13) and performing the indicated integration, 


we find 


(2.8.2 


Remembering that #2 = dx2/dt, we solve (2.8.2) for dt: 


dxe 


i= 
| 2h = Stat ~ 7 x! 
g m2” 2me 





(2.8.3 





Fig. 2.8.1—Duration of acceleration pulse for cushioning with cubic elasticity. See 


equation (2.8.10). 


Then, with the initial condition x. = 0 when ¢ = 0, the integral of (2.8.3) is 


dx2 


tor 


Ms 


dx 
of F-f 
y/ gh — ho x3 


To integrate (2.8.4), let 


Xe 
V k*(:3 “3 ay din ¥ + <” 


where 


es se VIi+B 


(2.8.4) 
a 
2m," 
(2.8.5) 
(2.8.6) 
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Fig. 2.8.2—Acceleration-time curves for cushioning with cubic elasticity. See 
equation (2.8.14). 


and B and d,, are as given in Part I: 
4Wehr 
a 


ie on : 
dn = iy 4/5 (-1 + V1+B)- 
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Then (2.8.4) becomes 


1 otf 5 a m 
ie os oe BD)’ (2.8.7 

in which the integral is the elliptic integral of the first kind (see Hancock 

“Elliptic Integrals,’ John Wiley and Sons, New York, 1917). In (2.8.7), 


we = wo (1 + B)’, (2.8.8) 


where wo = ~/ ko, mz is the radian frequency that would obtain if the cushion- 
ing were linear with spring rate ko. The motion for the linear case has a 
half period, or pulse duration to = 4/wy. The half-period (72) of the motion 


with cubic elasticity is twice the time required for x2 to increase from (0) 
to dm 


From 2.8.5 
[Z]2.—0 = 0, 
(2.8.9) 
Zissite - 1. 
Hence, from (2.8.7), the half-period is 
=f — = (2.8.10 
a ab Vi ee mone 


where K is the complete elliptic integral of the first kind. The duration of 
the acceleration pulse is therefore 2K/w.. We can define a radian frequency 
of the acceleration by 
1/4 
w= : = 4 = pao ~: St (2.8.11) 

The ratio wo/we (i.€., 72/79) is plotted in Fig. 2.8.1 which illustrates how the 
pulse duration decreases as the parameter B increases. Hence, for a given 
cushioning with cubic elasticity, the pulse duration decreases as the height 
of drop increases. This is in contrast to the linear case in which the dura- 
tion is independent of the height of drop. 

To find the acceleration #2 , we return to (2.8.7) and write it in the form 
of an elliptic function: 


SNnwt = Z. (2.8.12) 
Substituting the expression for Z given in (2.8.5) and solving for x2 , we find 
Xo = dmcn(wt — K). (2.8.13) 


Finally, differentiating (2.8.13) twice with respect to ¢, we find the accelera- 
tion to be 


#2 = widm[2k2sn2(wt — K) —1]cn(w-t — K). (2.8.14) 
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The ratio — #:/Gog is plotted in Fig. 2.8.2 against a radian coordinate 
(wot) for several values of B. It may be seen that, as B increases, the maxi- 
mum acceleration increases, the duration of the pulse decreases (see Fig. 
2.8.1) and the acceleration-time curve becomes bell shaped. For reference, 
the sinusoid for the linear case (B = 0) is plotted in the figure. 

Figure 2.8.2 is plotted for perfect rebound. If rebound does not occur, 
the curves continue, mirrored in the time axis, so as to form a periodic 
vibration of period 272. 


2.9 ACCELERATION-TIME RELATION FOR TANGENT ELASTICITY 


In this section the effect of tangent elasticity on the shape of the accelera- 
tion-time relation will be studied. The shape of the load displacement 
curve is given by 


Qk 1 TX 
P = ~* tan od, (1.4.3) 
T 2d> 


The system considered is again that shown in Fig. 1.2.1. Referring to the 
energy equation (1.2.13): 


9 
) 


as +[0 Pdx. = mz gh, (1.2.13) 


we substitute the above value of P and perform the indicated integration to 
obtain, for the velocity, 


8 ¢ Xo 
x» = 4/2! + ee dy i, og COs a (2.9.1) 


Then, as in Section 2.8, 
dx: 2 


= ['¢ dx. -[ 
4/ h + She < log cos 


and the half-period (72) of the motion is again twice the time required 
for x2 to increase from 0 tod,,. Hence 


(2.9.2) 


TXe 
2d, 


dx 


1m 
n= 2] —= 
4/ h += Sac * log yg 7 


> ds 


(2.9.3) 


where, from Section 1.9, 


dm = ams cos ' exp 7 
us 
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The radian frequency of the acceleration is 
Qo = 


and this is to be compared with the frequency 


¥ ho 


To me 


that would obtain if the cushioning were linear with spring rate ko. The 
ratio wo/we (i.e., 72/T 9) was obtained by numerical integration of (2.9.3) 
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Fig. 2.9.1—Duration of acceleration pulse for cushioning with tangent elasticity. See 


equation (2.9.3). 


and is plotted in Fig. 2.9.1 against the ratio d,/dy. The figure shows that 
for d»/dy < 1, the pulse duration varies almost linearly with d,/do. As the 
bottoming distance becomes larger than that required for linear cushioning 
with spring rate ko, the pulse duration approaches asymptotically the 
duration +/wo for the linear case. 

As dy/do decreases, the pulse duration becomes shorter, but the maximum 
acceleration increases, in accordance with equation (1.9.4) and Fig. 1.9.1. 
The shapes of the acceleration-time curves for several values of d,/do are 
illustrated in Fig. 2.9.2. They are more sharply peaked than the corre- 
sponding curves for cubic elasticity (Fig. 2.8.2) as might be expected from 
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the fact that the load-displacement curve for tangent elasticity rises more 
rapidly than that for cubic elasticity; that is, the bottoming is harder. 
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Fig. 2.9.2—Acceleration-time curves for cushioning with tangent elasticity. 


The curves of Fig. 2.9.2 were obtained by numerical integration of equa- 
tion (2.9.2), to obtain x2 as a function of ¢, following which these values 
were substituted in the equation 

2ko dy 1X2 
sf. + —— tan— = 0 
Tv 2d 
to obtain #2. It may be observed that the maximum values of the curves 
are the values dictated by equation (1.9.4). 
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In performing the numerical integrations of equations (2.9.2) and (2.9.3), 
it is found that the integrand becomes infinite when x. = d, since at this 
point the velocity is zero. In order to avoid this difficulty, it was assumed 
that, for a small distance in the neighborhood of d,, , the acceleration is 
constant with magnitude G,,g as obtained from equation (1.9.4). The 
procedure is described in further detail in Section 2.12. 

Figure 2.9.2 gives the acceleration-time curve for perfect rebound. [i 
rebound does not occur, the acceleration is a periodic vibration, each suc- 
cessive half period having the shape shown, with alternating sign. 


2.10 ACCELERATION-TIME RELATION FOR ABRUPT BOTTOMING 


By abrupt bottoming, we mean bilinear cushioning (Class D) as treated 
in Section 1.12. The load-displacement relation is (see equation (1.4.4) 
and Fig. 1.4.4) 


P = ko Xe 0 s X2 


IA 


d, ) 
ae, (1.4.4) 
P = kyxe — (Ry — Rods x2 5 d,| 


Considering, again, the system illustrated in Fig. 1.2.1, the equation of 
motion of m, , before bottoming, is 


My ko + ko x. = O OS n.d, (2.10.1) 
with initial conditions 
[welemo = 0, — [a] pno = VW 2gh. (2.10.2) 
The solution of (2.10.1) is then 


/ 
f= Vv 2gh SIN wo t, 0 
Wo 


IA 
2 
lA 


< d., (2.10.3) 


wy = /*. (2.10.4) 
mo 


The time (¢,) at which a» reaches d, is found from (2.10.3): 


Sa ie od 
mel Vin” sin 4,’ 


where 


(2.10.5) 


where 


! See Timoshenko, ‘Vibration problems in Engineering,” D. Van Nostrand Co.. New 
York, Second Edition (1937) page 123. 
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i.e., do is the displacement that would have been reached if the spring rate 
remained constant. 
The velocity of m, at time /, is 


[te]er, = V 2gh cos wols = /28h(1 _ 3). (2.10.6) 


If V/2gh > wod,, the displacement will exceed d, and the equation of 
motion becomes 


MoXo + Rpx2 — (ky — ko)ds = 0, vs (2.10.7) 


IV 


The solution of (2.10.7), with initial conditions 


[xe]r—e, =d, 


. . 2 (2.10.8) 
is 


xX. = as (eres _ i’) sin (wef + a — awyls) 


P (2.10.9) 
+ (1 — *) ds x, 2d, 
where 
— ky 
tan a = d? 
ky {2 — 1 f 
d; : (2.10.10) 


ky 
Op = o 
Ms 


By differentiating (2.10.3) and (2.10.9) twice with respect to ¢, the 
accelerations for the two regions are found to be 


IIA 


Xo = — Gog sin wol, 0sS2x2Sd,, (2.10.11) 


stat A € b\ A = 
2 = Gog f/e+e(s _ ) sin (asf +a wy ts), (2.10.12) 


ae 
. 2hko 
Go = VE ° (2.10.13) 


Typical shapes of the acceleration pulse represented by equations (2.10.11) 
and (2.10.12) are shown in Fig. 2.10.1. The curves are drawn for d,/dy = 


Vil 


where 
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0.5 and for several values of ky/ko. The peak values of the curves are th: 


same as given by equation (1.12.3). The curve marked k,/ko = 1 is the 


sinusoid of the linear case with duration 
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Fig. 2.10.1—Acceleration-time curves for cushioning with bi-linear elasticity. 
d./dy = 0.5. See equations (2.10.11) and (2.10.12). 
As before, if the package does not rebound, the acceleration shown is mir- 
rored in the time axis after each half cycle, to forma vibration of period 275. 
It is useful to know the duration of the complete pulse (aa’ in Fig. 2.10.1) 
and also the duration of bottoming (6d’ in Fig. 2.10.1). Calling the former 
tr. and the latter rz , we have, from equations (2.10.11) and (2.10.12) 
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ds 
do 
Fig. 2.10.2—Pulse durations for cushioning with bi-linear elasticity. See equations 
(2.10.15) and (2.10.16). 


TT. 7 d, ko 2 : l ky 
<a & /® “—— d? (2.10.15) 
a(S —1 
oe ko 2 —1 Ro 
thes Vi es / d3 (2.10.16) 
ky es 1 


V %\e 
These two equations are plotted in Fig. 2.10.2 for several values of ky/ Ry . 
2.11 ACCELERATION-TIME RELATION FOR HYPERBOLIC TANGENT 
ELASTICITY 
The relation between acceleration and time for hyperbolic tangent 
elasticity is found by the same procedure that was used for tangent elasticity 
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in Section 2.9. The system considered is that shown in Fig. 1.2.1 and th: 
load displacement curve of the cushioning is given by 


Ro wo 


P = Po tanh - .- 


Substituting the above expression for P in the energy equation (1.2.13), we 


find the velocity to be 


_. a Ro xe 
to = eh — ; 2.11.1) 
TV y/ gh a log cosh P, ( 
Then, as before, 
72 dx 
t= | — (2.11.2) 
0 v2 


and the half period (72) of the motion is twice the time required for x2 to 


increase from 0 to d,, , or 
dm d 
«3 os = (2.11.3) 


where, from Section 1.13, 


> ro? 
dy = - * cosh | exp (4). (1.13.3) 
2uU0 ~ 0 
The radian frequency of the acceleration is detined as 


us 
— = - 
T2 


and this is to be compared with the frequency 


o* aa 

yee V/! ms 

that would obtain if the cushioning had a constant spring rate equal to the 
initial spring rate (ko) of the hyperbolic tangent cushioning. The ratio 
wo/we (or 72/79) is plotted, in Fig. 2.11.1, against the dimensionless param- 
eter Po/W2Go (see Section 1.13). It may be observed that the pulse 
duration becomes very long when Po/W 2Gp is small, i.e., when the horizontal 
portion of the load displacement curve (Fig. 1.4.5) comes into play. The 
influence on the shape of the acceleration-time curve is illustrated in Fig. 
2.11.2. The curve marked Po/W2Go — © is the sinusoid for the linear 
case. For small values of Po/W 2G the curve approaches a square wave. 
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Fig. 2.11.1—Duration of acceleration pulse for cushioning with hyperbolic 
tangent elasticity. 


2.12 NUMERICAL PROCEDURE FOR FINDING ACCELERATION- TIME RELATION 
FOR CLAss F CUSHIONING 


When the load-displacement curve does not resemble one of Classes A to 
E, the acceleration-time relation may be found by numerical integration. 
Combining the energy equation, 





= + I Pp dx» = mogh, (1.2.13) 
2 0 
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with the equation relating time and velocity, 


72 dx 
as I _ (2.12.1 


Xe 
we find 


/wW. {7 dx» 


die / “4 | re . )) 
4 = fi - | P dxz eee 
0 


As an example, consider the problem of a 15-pound article supported on 
cushioning with the load-displacement curve shown in Fig. 2.12.1. The 
package is to be dropped from a height of 3 feet. The computations are 


—- ————\ 7 


} 
} 





4 
Fig. 2.11.2—Acceleration-time curves for cushioning with hyperbolic tangent elasticity. 


given in detail in Tables III and IV. The headings of Columns (1) to (8) 
of Table III are the same as in Table I, Section 1.15. A, is the integral 
under the radical of equation (2.12.2). Column (10) of Table III is the 
integrand of Equation (2.12.2), i.e., it is proportional to the reciprocal 
of the velocity expressed as a function of displacement. The function 
is plotted in Fig. 2.12.2 and its integration is performed in Table IV. In 
columns (11), (12) and (13), intervals of x: are chosen to suit the shape of 
the curve. The values for column (14) are taken from column (10). Col- 
umns (15) and (16) perform the same operations on the integrand 
(Weh — A,) ' that are performed in Columns (5) and (6) of Table III on 
the integrand P. 
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P (POUNDS) 





Oo 1 2 
X> (INCHES) 
Fig. 2.12.1—A load displacement curve for Class F cushioning. 


TABLE III 


1) (2) (3) 4) 5 6 7 8 ) 10 
P I A(x Pp . An 1 I : 1 
n A(x2)n Y2)n n 3 nt Pa | hogy ” Pda hy, = i Gn = ii Wek — A, Via : 
0 0 0.0 | 0 0 0 0 0 540 (0.0431 
1 .20 0.2 |105 10.5 10.5 0.7 7.0 529.5 0.0435 
2. .20 0.4 (155 26.0 36.5 2.4 10.3. 503.5 0.0446 
3 .20 0.6 \192 34.7 ya Be 4.8 12.8 468.8 0.0462 
4 .20 0.8 (217 40.9 112.1 eau 14.5 427.9 0.0483 
5 .20 (1.0 (237 45.4 157.5 10.5 15.8 382.5 0.0511 
6) .20 '1.2 (257 49.4 206.9 13.9 17.1 333.1 0.0547 
7| .20 |1.4 1277 52.9 259.8 17.3 18.5 | 280.2 (0.0597 
8 .20 1.6 |305 58.2 318.0 252 20.3 222.0 0.0671 
9, .20 1.8 (342 64.7 382.7 23.5 22.8 157.3 0.0798 
10, .20 (2.0 (392 73.4 456.1 30.4 26.1 83.9 0.109 
11. .05 (2.05 405 19.9 476.0 31.8 27.0 64.0 0.125 
12 .05 (2.10 422 20.7 496.7 33.2 28.1 43.3 0.152 
13, .05 2.15 |440 21.6 518.3 34.6 29.4 21.7 0.215 
14.01 (2.16 (445 4.42 S22. 34.8 29.7 17.3 0.240 
15) .O1 (2.17 |450 4.48 527.2 35.2 30.0 12.8 ().279 
16. .O1 (2.18 455 4.52 531.7 35.5 30.3 8.3 0.347 
17) .O1 |2.19 457 4.56 536.3 35.8 30.5 ce | 0.521 
2 4.60 540.9 36.1 30.8 0 va 


18, .O1 2.20 462 
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TABLE IV 


11 12) 13) 14 (15) (16) 17) 
fy, = 4 D, = t= 
Pa A(x2)n x2)n 1 X2)n ee ae x2)n dx. JW; 
VWehk — An . [ Wack, OVS 
0 0 0 0.0431 0 0 0 
1 0.4 0.4 0.0446 0.0175 0.0175 0.0024 
2 0.4 0.8 0.0483 0.0185 0.0360 0.0050 
3 0.4 ee 0.0547 0.0206 0.0566 0.0079 
4 0.4 1.6 0.0671 0.0243 0.0809 0.0112 
_ 0.2 1.8 0.0798 0.0149 0.0958 0.0133 
6 0.2 2.0 6.109 0.0189 0.1147 0.0160 
7 0.1 re 0.152 0.0131 0.1278 0.0178 
8 0.05 2.15 0.215 0.0092 0.1370 0.0190 
9 0.03 2.18 0.347 0.0083 0.1453 0.0202 
10 0.01 2.19 0.521 0.0043 0.1496 0.0208 
11 0.01 2.20 a 0.0221 











ie} 2 4 8 1.2 16 2.0 24 


Fig. 2.12.2—Plot of Column (3) vs. Column (10) of Table IIT. 
A difficulty arises because the integrand (Wh — A,)* becomes infinite 
for the maximum displacement (see Column (14)). This is avoided by 
assuming that the acceleration is constant in the last interval? and has the 


* Timoshenko, “Vibration Problems in Engineering,” D. Van Nostrand Co., New York, 
Second Edition (1937) page 123. 
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value given in Column (8), Table III, for the maximum height of drop. 
Then, 


A(x2)n = 3GmgAl (2.12.3) 


DA x» 
At = 4f a (2.12.4) 
Gm & 


or 





1°) 01 02 .03 04 


t 


Fig. 2.12.3—Acceleration-time curve (for the cushioning shown in Fig. 2.12.1) obtained 
by numerical integration. 


In the present instance, 
(Axe), = 0.01 inches 
Gng = 30.8 K 386 = 11900 in/sec.* 


Hence, from (2.12.4), At = 0.0013 sec. and the last entry in Column (17) 
is obtained by adding this value of At to the preceding entry. 

The final curve of acceleration vs. time is obtained by plotting the entries 
of Column (17) against the entries of Column (8), Table III, for correspond- 
ing values of x2. The result is shown in Fig. 2.12.3. 
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PART III 
AMPLIFICATION FACTOR 


3.1 INTRODUCTION 


If the maximum acceleration, of the packaged article as a whole, is 
reached very slowly, the severity of the disturbance experienced by a 
structural element of the packaged article is very nearly proportional to the 
maximum acceleration. Roughly speaking, ‘“‘very slowly” means that the 
time, during which the acceleration undergoes a major change in magnitude, 
is long in comparison with the natural period of vibration of the element 
under consideration. When this is so, no transient vibration is excited in 
the element. The displacement response of an element under very slowly 
varying conditions is called the ‘‘static response”. Under more rapidly 
varying conditions the dynamic response to the same maximum acceleration 
may be greater or less than the static response. The ratio (A) of the maxi- 
mum dynamic response to the static response is called the amplification 
factor. In general, for a given acceleration disturbance, very low-frequency 
elements have amplification factors less than unity, while the amplification 
factors are greater than unity for elements whose natural frequencies are 
near or above the disturbing frequencies. The numerical value of the 
amplification factor depends not only on the manner in which the disturbing 
acceleration varies with time, but also on the “reference acceleration”’, i.e., 
the value of acceleration for which the static response is calculated. Usually 
the reference acceleration chosen for calculating the static response is the 
maximum value (G,,) of the disturbing acceleration. However, when 
special circumstances are being investigated, such as the effect of damping 
or abrupt bottoming, the reference acceleration is taken to be Go , which is 
the acceleration that would be reached if the damping or bottoming were 
absent. In such cases the amplification factor includes both the effect of 
rate of change of acceleration and the effect of the special conditions. 

When the reference acceleration is G,, the amplification factor will be 
denoted by A,, and when the reference acceleration is Go the amplification 
factor will be denoted by Ay. The symbol G, will be used to designate the 
slowly applied acceleration that would produce the same maximum dis- 
placement as the transient acceleration, i.e., G. = AnG» or G, = AgGo. 
The symbol G, will be used to denote the safe value of G, , for an element 
of the packaged article, as determined by a strength test or by calculation. 
In specifying G, some judgement is required to take into account the effects 
of plastic deformation in comparing tests made on greatly different time 
scales. Good judgement is also necessary in deciding whether or not the 
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assumptions listed in Section 0.2 are valid in each application. The general 
procedure for using amplification factors is as follows. We first find the 
value of the reference acceleration (in units of number of times gravity) from 
Part I. From Part II we find the properties of the acceleration-time rela- 
tion which give us the information required for entering one of the curves 
of Part III and finding the amplification factor. Then, the product of the 
reference acceleration and the amplification factor (4,,G,, or ApG») is a 
number (G,) by which the weight of the structure is to be multiplied when 
calculating its deflection or stress by the usual static methods of elementary 
strength of materials. Alternatively, G, must be found not to exceed G, . 
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Fig. 3.2.1—Idealized system used in calculating amplification factors for linear undamped 
cushioning with perfect rebound. (a) initial position, (b) first contact with floor. 


In the following sections the amplification factors for typical transient 
accelerations encountered in package drop tests are calculated. The ampli- 
fication factor curves that are plotted are entirely analogous to the familiar 
“resonance curves” for steady sinusoidal vibration, except that in this 
case the disturbing forces are transients of various shapes. It will be seen 
from the curves that the maximum acceleration, as calculated by the 
methods of Part I or as measured by an accelerometer, is not always a true 
measure of the severity of the disturbance. 


3.2. AMPLIFICATION FACTORS FOR A HALF-SINE-WAVE PULSE 
ACCELERATION 


The first case to be treated is the response of an element of the packaged 
item to the transient acceleration that would occur in a package with linear 
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undamped cushioning and perfect rebound. Figure 3.2.1 illustrates the 
idealized system, and it may be noted that the mass ms is omitted, as is 
required for perfect rebound (Section 2.3). At first we shall consider that 
the mass m is undamped and later we shall consider the effect of damping 
in this element. 

The mass m, is taken to be small in comparison with mz, so that the motion 
of the latter is the same as we found it to be in Section 2.2 where m, was not 
considered. Hence the acceleration of mz is a half-sine wave pulse: 


dt. = —wo V2ghsinwst, (0 = t = x/w). (3.2.1) 
The equation of motion of my is 
mx, + k(x; — Xo) = 0. (3.2.2) 
Let x be the relative displacement of m, with respect to my , i.e., 
X= X— mM. (3.2.3) 


x is proportional to the force in the spring (Aix) and to the acceleration 
of m, and hence is proportional to the deflection, strain and stress in the 
element which the system m;, , k; represents. 

Substituting (3.2.3) in (3.2.2), we find: 


mx + kx = —m ko. (3.2.4) 


This equation holds for the duration w/w, of the pulse #:. The initial 
conditions for x are 


[x].no0 = [Z],0 = 0 (3.2.5) 
so that the solution of (3.2.4) is 


ins wr V 2¢h [2 sin w,f — sin wt], (0 ztiz =), (3.2.6) 


AN 
| 


w; — wi Lor 


It may be seen that x is composed of a forced displacement at the accelera- 
tion frequency w2, on which is superposed a free vibration at the natural 
frequency, w; , of the element. The maximum value of the relative displace- 
ment is 


212 x) 
We : (3.2.7) 


in which n is a positive integer chosen so as to make the sine term as large 
as possible while the argument remains less than 7. 

(3.2.7) gives the maximum dynamic response of the element m, during 
the interval of impact. To find the amplification factor we must compare 
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Xmax With the ‘‘static response” i.e. with the value (x,,) that x would have 
if the acceleration #2 reached the same maximum value (w21/ 2gh) in a very 
long time. The resulting value may be found from (3.2.4) by omitting 
the transient term m,%. Then 


/ m 
Xs = we V 2gh - 
1 


or 
we 
- Py 1] 
Xe = “9 V2gh. (3.2.8) 
@) 


The amplification factor for the interval 0 = ¢ < m/w. is then 


W) 


Se we . inex 
= “= sin — (0 


=f ey 
We 


_— 


1/2). (3.2.9) 


All 
Ail 


«itm 


Xst @1 
We 


It should be observed that A,, depends only on the frequency ratio w;/w» . 
That is, since w;/w: = T2/t,, the amplification factor depends only on the 
ratio of the pulse duration to the half period of vibration of the element. 

Thus far we have studied only the motion in the interval 0 = ¢ = r/ae. 
We must not, however, overlook the possibility of larger displacements of m, 
with respect to m, occurring after rebound. In fact, examination of (3.2.6) 
reveals that x has no maximum in the interval 0 <= ¢ = r/w2 whenw; < we. 
It is very likely, then, that larger values will occur at later times. 

After rebound, m, executes free vibrations with respect to m.. We have 
to compare the magnitude of xmax, in the interval 0 = ¢ = w/we , with the 
amplitude of the free vibration. Calling the relative displacement during 
free vibration x’ and measuring a time coordinate ¢’ from the instant the 
package leaves the floor, we have 


mx’ + kyx’ = 0, (3.2.10) 
with initial conditions 


[x’] p20 tig [x] mofo ’ 
(3.2.11) 


[x] e-mo = [tee we 


The solution of (3.2.10) with initial conditions (3.2.11) is 


vi /sen( 1 + cos“) . 
: igen, Si (wr > 
(3.2.12) 


z= 


We 


wi(w} — w}) 


us 
we 











426 BELL SYSTEM TECHNICAL JOURNAL 


Then 
WI TW, 
2— cos D 
YX me Ww a) a 
i, Se SS oar. ‘: : (3.2.13) 
Xst 1 wy We 


ws 
We find, on comparing (3.2.13) with (3.2.9) that for a; < w. equation 
(3.2.13) gives the larger value of 4,,, while for w; > w2 equation (3.2.9) 
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Fig. 3.2.2—Amplification factors for linear undamped cushioning with perfect rebound. 
See Fig. 3.2.1 and equations (3.2.9) and (3.2.13). 


gives the larger value of A,,. That is, when the duration of impact is 
shorter than the half-period of vibration of the element, the maximum 
displacement (and stress) in the element occurs after the impact is over. 

The curve marked 8, = 0 in Fig. 3.2.2 is a plot of the largest value of A», 
from (3.2.9) and (3.2.13) with the frequency ratiow:/w: as abscissa. (3.2.13) 
was used for wi/we2 S 1 and (3.2.9) for wi/w2 2 1. The maximum value 
of A, is 1.76 and occurs at w;/w2 = 1.6. Hence, at this frequency ratio, the 
deformation of the element is 1.76 times as great as would be expected from 
a calculation using the maximum value of acceleration alone as in Part I. 
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On the other hand, for frequency ratios w:/w: < 0.5 the severity of the 
shock can be very much less than might be expected from the calculations 
of Part I. For very small values of w;/w2 the amplification factor may be seen 
from (3.2.13) to be equal to 2w;/w2. For large values of w;/w: (stiff elements) 
Fig. 3.2.2 shows that the amplification factor is very nearly unity and the 
methods of Part I can be used without additional calculation. 

When damping of the element of the packaged article is considered, the 
amplification factors are less than without damping. The applicable 
equations of motion during and after impact are obtained by inserting 
velocity damping terms in (3.2.4) and (3.2.10): 


me + ox + kx = —m, xo, O2tZ (3.2.14) 
we 
-f a, ae * = v g 
m* +ox% + kx’ = 0, iS —. (3.2.15) 
we 


If we express the damping of the element m, as the fraction of critical 
damping 


1 
= (3.2.16) 
Pm ka 
(as in Section 2.5) equations (3.2.14) and (3.2.15) become 
#+ Bot+wix=—h, OF LS, (3.2.17) 
We 
ef -f | = 7 > 
# + Biot’ + wiv = 0, i alban (3.2.18) 
We 


The amplification factors for equations (3.2.17) and (3.2.18), with boundary 
conditions (3.2.5) and (3.2.11), respectively, were obtained on the Westing- 
house Mechanical Transients Analyzer* for 8, = 0.005, 0.01, 0.05, 0.10, 0.30, 
0.50 and 1.00. The curves are shown in Fig. 3.2.2. 


3.3. APPLICATION OF HALF-SINE-WAVE AMPLIFICATION FACTORS 


As an example of the use of the amplification factor curves of Fig. 3.2.2, 
let us consider the following problem: 


3 Arrangements for performing these calculations were made through the courtesy of 
Mr. A. C. Monteith, Manager of Industry Engineering, and Mr. C. F. Wagner, Manager 
of Central Station Engineering, Westinghouse Electric and Manufacturing Co. Dr. G. D. 
McCann, Transmission Engineer, was in immediate charge of the project. For a descrip- 
tion of the analyzer see ‘““A New Device for the Solution of Transient-Vibration Problems 
by the Method of Electrical-Mechanical Analogy” by H. E. Criner, G. D. McCann and 
C. E. Warren, Journal of Applied Mechanics, Vol. 12, No. 3 (1945) pp. A-135 to A-141. 
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It is required to judge the suitability of a proposed package for a large 
vacuum tube weighing 10 pounds. Strength tests have been made on the 
tube in a shock testing machine which produces a half-sine-wave acceleration 
pulse of 25 milliseconds duration. The weakest element of the tube is 
found to be the cathode structure, for which the safe maximum acceleration 
in the drop testing machine is 200g. The cathode structure has a natural 
vibration frequency of 120 cycles per second and has 1% of critical damping. 
The proposed package has essentially linear, undamped cushioning with a 
spring rate of 3300 pounds per inch and an available displacement of } inch. 
The outer container weighs much less than the tube so that the package may 
be expected to rebound. Is the cushioning suitable for protecting the 
cathode in a drop of 5 feet? 

First find the maximum G that the tube will experience in a 5 ft. drop of 
the package (equation 1.3.3): 


Kate 4/2 y/? LORS. 


The accompanying maximum displacement is, from equation (1.3.4), 


2h _ 2X 60 
dm = me = 0.6 in. 

The available displacement (} inches) is therefore sufficient and the maxi- 
mum acceleration (199g) is slightly less than the safe maximum acceleration 
(200g) found with the shock testing machine. However, before the cush- 
ioning is approved it is necessary to investigate the frequency effects. The 
duration of acceleration in both the shock machine and in the package 
must be considered. 

The amplification factor for the element tested in the shock machine is 
found as follows. First find the frequency corresponding to the 25 milli- 
second pulse: 


1 


Se = 9025 


= 20 c.p.s. 


The ratio of the element frequency to the shock machine frequency is 


fram 


Entering Fig. 3.2.2 with w;/w. = 6, we read, from the curve 8; = 0.01, An = 
1.14. The 200g test in the shock machine is, therefore, equivalent to a 
slowly applied acceleration of G, = 200 XK 1.14 = 228,. 
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To find the corresponding quantity for the package drop, first find the 
cushion frequency: 


ee a 
= OR m, 24 10 ee 


The ratio of the element frequency to the package frequency is therefore 


fo we 57 
Entering Fig. 3.2.2 with wi/w. = 2.1 we read, from the curve 8; = 0.01, 
Am = 1.59, The 199g acceleration pulse in the package drop is therefore 
equivalent to a slowly applied acceleration of G. = 199 & 1.59 = 316g. 
This is almost 40% in excess of the value (228g) found to be safe from the 
shock machine data. The cushioning is therefore judged to be inadequate. 
The procedure for finding the correct spring rate for the cushioning is as 
follows. It is known that we must have 


G. #6.. 


Therefore, take 


Now 
a= of Bony 
Therefore 
Anw, = 409 rad/sec. 
Also 


w, = 2n X 120 = 754 rad/sec. 


Then, with successive trial values of we, we calculate w;/w2, enter Fig. 3.2.2, 
read the corresponding value of A,, from curve 8; = 0.01 and test to see if 
the product A,,w2 = 409. The combination which satisfies the test is found 
to be 


we = 280 rad/sec. 
w/we = 2.69 


Am = 1.47. 
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Then 
‘ ))? 0 ; 
ky = wom, = cia. = 2030 lbs./in. 
Gy, = / > = 155 
? 
jo 
Ge 


Hence the spring rate of the cushioning should be reduced from 3300 Ibs./in. 
to 2030 Ibs./in. and the available space should be increased to accomodat, 
the 0.77 inch maximum displacement before bottoming. 
3.4. SpecraL TREATMENT OF STRONG, Low FREQUENCY ELEMENTS 
The product of the amplification factor (4,,) and the maximum accelera- 
tion (G,,) must be equal to or less than the maximum allowable slowly 
applied acceleration (G,): 
G, = AnGm = G; . (3.4.1) 
For frequency ratios 
1 
<z, (3.4.2) 
Figure 3.2.2 shows that, approximately, 
w 
Am = i (3.4.3) 


We 


for 8, < 0.10. When this is so, we may combine (3.4.1) and (3.4.3) to 
obtain the criterion 


6. &.. (3.4.4) 
We 
Now, 
—_— / ahs. 4/28. (3.4.5) 
Wu 2 g 


Hence the criterion (3.4.4) may be written as 


‘) ae 
2w) 7 = G, (3.4.6) 


; 1.1 G, 
tS Vh (3.4.7) 


or 


where / is in inches and /; is in cycles per second. 
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It may be observed that (3.4.7) is independent of the properties of the 
cushioning. Hence, as long as (3.4.2) is satisfied, any cushioning at all 
may be used for an element that satisfies (3.4.7) regardless of the magnitude 
of the maximum acceleration G,,._ In particular, rigid mounting is suitable 
for such an element. The only precaution to be observed is that the maxi- 
mum acceleration and duration must not be unfavorable for other elements 
of the packaged article. 

Example: A 9-pound vacuum tube has an anode structure for which the 
safe maximum acceleration is 200g as determined in a centrifuge test. 
The natural vibration frequency of the anode is 35 cycles per second and the 
damping is 1% of critical. What cushioning around the tube is required to 
protect the anode from damage in a package drop of 3 feet? 

Calculate 


1G, 11X20 0 
ag 7, il. )./ C.p.s. 


This is greater than f; = 35 c.p.s and hence any cushioning is safe for the 
anode. The results of calculations for cushioning with spring rates of 50, 
500, 5000, 5 X 10° and 5 X 10° pounds per inch are given in the following 


table: 


@) 


ko(Ibs./in.) Gm = Aw AmGm 
50 20 4.74 1.12 22 
500 63 1.47 1.65 106 
5 X 10° 200 474 0.9 180 
5 &K 10° 2,000 .0474 0.09 180 
5 x 107 20,000 .0047 0.009 180 


In each case the product of A,,G,, is less than the allowable 200 and, as long 
as the combination of G,, and the amplification factors for other elements 
does not exceed the allowable A,,G,, for those elements, the cushioning is 
suitable. The precaution to observe is that higher-frequency elements 
shall not have amplification factors such that 4,,G,, may be excessive for 
them. 


3.5 AMPLIFICATION FACTORS FOR DAMPED SINUSOIDAL ACCELERATION 


If the outer container of the package is heavy enough (see Section 2.3) 
there will be no rebound and the packaged item will vibrate in the cushion- 
ing afterimpact. For linear cushioning with velocity damping, the accelera- 
tion produced by the vibration will be a damped sinusoid (equation (2.5.5) 
and Fig. 2.5.2). The system to be considered is shown in Fig. 3.5.1. 
To determine the effect of the damped vibration of m, on the mass m, , 
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Fig. 3.5.1—Idealized system for linear damped cushioning with no rebound. 
(a) initial position, (b) first contact with floor. 
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Fig. 3.5.2—Amplification factors for linear damped cushioning with no rebound. 
B2 = 0.005. See equations (3.5.1) and (3.5.2). 


we note that the equation of motion and initial conditions are identical 
with (3.2.17) except that for the acceleration # we use the damped sinusoid, 
equation (2.5.5), instead of the half-sine pulse (3.2.1). The solution of 


(3.2.17), i.e. the relative displacement (x; — x2) of m, with respect to mz, 
is 
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3 V2gh , —s.0 ag ee 
pe VER mts sin (ult + y — 8) + Bain (wt — 7 — Dd) 
2w5 W} 


w 
na 
— 


— ¢ 14 sin (wt + y — 8)— Bsin (wet + 7 + O)]! (3. 


ki 
gy = 

my, 

he 
@ = 

Me 


wo = oV1 — 8 
w =u V1 — B 


where 


l/A = V (Bowe six Bia)" + (wy = ws)" 
1'B = V (Bow. — Bun)? + (w) + 3)? 
283 — 1 
tan y = 26, V1 — B 
; wr — ws 
an 0 = 
Bow. — Biwi 
tan ¢ = = + 


Bowe — Biwi 


The relative displacement of m, with respect to mz is seen to consist of a 
forced, damped vibration (w: , 82) on which is superposed the free damped 
oscillations (w; , 81) of m,. 

The amplification factor 


9 
x W) Ymax a 
{j= = : (3.5.2) 


Xe 2 WV Agh 

is plotted in Figs. 3.5.2 to 3.5.7 for six values of 8; and six values of 82 . 
These curves were obtained by direct solution of the differential equation 
on the Westinghouse Mechanical Transients Analyzer.’ The amplifica- 
tion factor in this case includes the effect of damping; i.e., the reference 
acceleration is not the maximum acceleration of m2, but is the maximum 
acceleration that m2 would reach if the damping £2 were zero. Conse- 
quently, the amplification factors for large values of w;/w2 do not approach 
4 See footnote, Section 3.2. Only enough data were obtained with the analyzer to find 
the general shapes of the curves, so that the fine structure is not revealed. Checks on 


the analyzer results were made by computing Ao from equations (3.5.1) and (3.5.2) for 
w1/w2. = 1, B, = Bo; w/o. = 0; w/w. ©. 
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Fig. 3.5.3—Amplification factors for linear damped cushioning with no rebound. 
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Fig. 3.5.4—Amplification factors for linear damped cushioning with no rebound. 


B: = 0.05. See equations (3.5.1) and (3.5.2). 
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See equations (3.5.1) and (3.5.2). 
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unity. For example, the curve for 8; = 0.005, 8; = 1 (Fig. 3.5.7) approaches 


a value of nearly four as w;/w,;— ~. 


The factor four is composed of two 
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Fig. 3.5.5—Amplification factors for linear damped cushioning with no rebound. 
B2 = 0.1. See equations (3.5.1) and (3.5.2). 
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Fig. 3.5.6—Amplification factors for linear damped cushioning with no rebound. 
Bz = 0.5. See equations (3.5.1) and (3.5.2). 
factors of two. The first arises from the fact that the maximum value of 
acceleration, for 82 = 1, is twice the value that would be reached if 8» 
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were equal to zero (see Fig. 2.5.3). The second factor (of nearly two) js 
due to the fact that the maximum acceleration is reached at time ¢ = () 
when 82 = 1 (see Fig. 2.5.2) and the response of an almost undamped sys- 
tem (8, = 0.005) to a suddenly applied and subsequently maintained 
acceleration is double the response to a slowly applied acceleration (see 
curve (a) Fig. 3.8.1). For 8; > 0 and 6. < 1 the amplification factor is 
less than four, as w;/w. — *, in accordance with the curves plotted in Fig 
3.5.8. 

Example: A 1.5-pound vacuum tube is to be packed in a container whose 
estimated weight will be at least 50 pounds. The cathode structure of the 
tube has a natural frequency of 25 c.p.s. with damping 0.5% of critical 
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Fig. 3.5.7—Amplification factors for linear damped cushioning with no rebound. 
B2 = 1.0. See equations (3.5.1) and (3.5.2). 


and its safe acceleration, as determined in a centrifuge, is 90g. What 
spring rate of cushioning is suitable for protecting the cathode in a drop of 
five feet? It is specified that the cushioning shall have damping 50% of 
critical. 

Assuming linear cushioning, the spring rate that would be prescribed, by 
considering maximum acceleration alone, is 


_ WrGm _ 15 X (90)" 


s Festi = 10 ./in. 
ko ah 2x 60 101 Ibs. /in 


Considering damping, Fig. 2.5.3 shows that 50% of critical damping does 
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not change G,,. To find the amplification factor we must first decide if the 
package will rebound. With 50°% of critical damping, the maximum ac- 

40 
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Fig. 3.5.8—Limiting values of amplification factors for linear damped cushioning with no 
rebound. w;/w2— ©. See equations (3.5.1) and (3.5.2). 


celeration on the first upstroke is 0.164 G,, (see Section 2.6 and Fig. 2.6.1). 
Then, 0.164 X 90 X 1.5 = 22 lbs. which is less than the estimated weight 
of the outer container. The package will not rebound and Fig. 3.5.6 
should be used for the amplification factor. 
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The frequency of vibration of the tube in its cushion will be 


w. = /* = ” de x 2 159 rad. /sec. 
Me 1.5 


Hence w;/w, = 2m XK 25/159 = 0.99 and, from Fig. 3.5.6, Ao = 1.4. Hence 
G. = 90 XK 1.4 = 126, which is greater than the allowable G, = 90, so 
that the 101 lb./in. cushion is unsatisfactory. 

To obtain satisfactory cushioning, set 


ApGo = 90, 
that is 


90 


1, = = 5¢ ‘d ve 
Agwe NE: 159 rad./sec 
g 


Noting that w; = 2x X 25 = 157, we find from Fig. 3.5.6 that there are two 
values of w. (90 and 600 rad/sec.) that satisfy the criterion Ayws = 159 
rad/sec. The first gives 


we = 90 rad/sec. 
Ao = 1.8 . 


ky = 31.5 lbs./in. 


Go = 50 
G, = 90 
d,, = 2.4 in. 


The second gives 
we = 600 rad/sec. 
Ay = 0.27 
ky = 1400 Ibs./in. 
Go = 335 
G,. = 90 
d», = 0.36 in. 


The second solution requires less space for cushioning than the first but 
should be used only if the remainder of the tube can endure the high ac- 
celeration of 335g. Otherwise the 50g package should be used . 
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3.6 AMPLIFICATION FACTORS FOR THE PULSE ACCELERATION OF 
CuBIC CUSHIONING 
In a rebounding package with undamped Class B cushioning, the pack- 
aged article (mz) will undergo a pulse acceleration of duration m/w as given 
by equation (2.8.11). The shape of the pulse is illustrated in Fig. 2.8.2 
and its functional form is 


2 2 , ; 
, « Ss | 28 ot CA a K) ts 1 ]or ea “ K). (2.8.14) 
™ Tv vs 


To determine the influence of the shape and duration of this pulse on the 
amplification factor, we proceed as before by substituting (2.8.14) in the 
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Fig. 3.6.1—Idealized system used in calculating amplification factors for non-linear, 
undamped cushioning with perfect rebound. 


differential equation governing the relative displacement (x = x; — x2) 
between m, and mz (see Fig. 3.6.1): 
+ wix = —d#e. (3.6.1) 


With boundary conditions x(0) = #(0) = 0, the solution of (3.6.1) may be 
written as 


t 
t= i #2 (A) sin aw, (A — ¢t) dd (3.6.2) 


@1 “0 


and the maximum value of « may be expressed by 


tm 
Sines: = [ ¥2(A) sin wi(A — tm) dr, (3.6.3) 
wr Jo 


where /,, is the time at which the largest value of x occurs. 
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The amplification factor, in this case, will be taken as the ratio of x 
to the relative displacement (x,,) resulting from a slow application of th: 
maximum value of #:. From (3.6.1), 


“ ’ 
2 | max Gm g 
u 


9 ? 


Ww} W] 
where G,, is given by equation (1.5.6). Then 


Vinax _ Vmax W] 
fu a= — 


Nee Gm 8 


Amplification factors for undamped cushioning with cubic elasticity. Perfect 
rebound. See equation (3.6.6). 


@) 


Am 


tm 
[ ¥o(A) sin wi(\ — ty) dr. (3.6.6) 
0 


iil gs 


A, was evaluated, mostly by graphical methods, for four values of B 
(0, 2, 20 and ~) and the results are plotted in Fig. 3.6.2. Observing that 
B = 0 corresponds to linear cushioning, it may be noted that cubic non- 
linearity in the cushioning does not change the amplification factor by 
more than 35% even in the most extreme case (B — ~). The severity 
of the shock, however, may be much greater for the cubic cushioning than 
for linear cushioning with a spring rate equal to the initial spring rate (ho) 
of the cubic cushioning. This is because 4,, is multiplied by G,, to obtain 
G, and, for large values of B, G,, may be much larger than the maximum 
acceleration for the linear case. In other words, in comparing Class B 
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with Class A cushioning the difference in maximum acceleration, rather than 
the difference in amplification factors, is usually more important. 

Example: Consider the example given in Section 1.6 and let it be required 
to determine the effect of pulse duration on a cathode structure with a 200 
c.p.s. natural frequency of vibration. In Section 1.6 we found that 


B = 5.4 Ry = 255 
Go = 28.6 y = 108. 
Gm = 55 


With B = 5.4, enter Fig. 2.8.2 and find 


oo = 0.88. 
we 
Now 
/p ‘255 X 386 
wy = 0§ _ ———— = 66.1 rad./sec. 
W. 2Ei9 
Hence 
_ 66.1 


5 rad./sec. 


~I 


ag = _ 

0.88 
Then, with w;/w2. = 27 X 200/75 = 16.7, enter Fig. 3.6.2 and find A,, = 
approximately 1.0. Hence G, is abou’. the same as G, and the conclusions 
reached for this problem in Section 1.6 are not altered. 


3.7 AMPLIFICATION FACTORS FOR ABRUPT BOTTOMING 


The amplification factors for bilinear elasticity have not been computed 
in complete detail. They can be obtained approximately by using the dura- 
tion curves (Fig. 2.10.2) and the amplification curves for the linear case 
(Figs. 3.2.2 and 3.5.2 to 3.5.7). It is useful, however, to calculate the am- 
plification factors for extremely abrupt bottoming (4, — ~) to obtain a 
general understanding of the accompanying phenomena. 

The system to be considered is illustrated in Fig. 3.7.1. It is assumed 
that the impact between m, and the base (occurring at ¢ = ¢,, %2 = d,) 
has a coefficient of restitution of unity. Hence mp will strike the base with 


velocity 
[te]imt, = p/2x(1 _ ) 
d? 


(see equation (2.10.8)) and leave it at a velocity of the same magnitude but 
opposite sign. Perfect rebound of the whole package is also assumed. 
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The acceleration pulse will then look like the curve marked k,/ky — » 
in Fig. 2.10.1. 


There will be three regions in which to consider the relative displacement .: 
Region 1 0<t<t, 
Region 2 t,<t< 2, 


Region 3 t> 2t, 


The relative displacement (x = x; — x2) of m, with respect to me will have 























f 2 


Fig. 3.7.1—Idealized system representing abrupt bottoming. 


the same functional form for Region 1 as in the linear case (see equation 
(3.2.6)), and the amplification factor is, by analogy with (3.2.9), 


O<t<h, 


For Region 2, we use the differential equation 


+ wrx = wo V 2gh sin wo(t —_ 2t,) (3.7.2) 


and, as initial conditions at ¢ = ¢,, we use the terminal conditions for 
Region 1 with the sign of [%2],.., reversed. The amplification factor for 
this region is found to be (by the same’ method as in Section 3.2): 


where 
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w? 


ge V A + ‘B? sin (witm + 7 — wile) 


1 in( - (3.7.3) 
“ye SIN | Wo , — wot, — sin : 
rs “*) dy 


<4 < 2h 


where 


_ os 
woV 2h 


n 


and t, is the root of 


3 


SAE FB 008 (itm + ! 
we V/ 2gh A? + B? cos (w) n — wile) 
1 
_ FX? Cos (coin — wot, — sin’ *) = 0 
1 ~ ES " 
(=) 
that yields the largest value of Ao in equation (3.7.3). Region 3 is gov- 
erned by 
# + wix = 0 (3.7.4) 

and the initial conditions are the terminal conditions of Region 2. By 
the same method as was used in Section 3.2, we find 


7 ae 
(2) " — @ 
t> 2t,. 


The largest value of Ay from equations (3.7.1), (3.7.3) and (3.7.5) is 
plotted against w:/w in Fig. 3.7.2 for several values of d,/do . 
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Notice that the amplification factor is Ao rather than An. That is, 
the reference acceleration is Go rather than G,,. This is necessary because 
G, is infinite in the present instance. Hence Fig. 3.7.2 cannot be com- 
pared directly with Figs. 3.2.2 and 3.6.2. However, it is interesting to 
observe that, for w:/w2 < 0.5, (low frequency elements) abrupt bottoming 
has no harmful effect. For high-frequency elements, the severity of bottom- 
ing is very great even when very nearly all of the required space (d,) is 
available. For example, if 90°% of the required space is available (d,/d) = 
0.9) and the frequency of the element is ten times the package frequency, 

ni bee GAR ae | | 
w——p— tf 


| 
| 
} 

} 


| 
| eee 


20,-—___-—__-+— 


paaee 


























LA 
Op 
Fig. 3.7.2—Amplification factors for abrupt bottoming. See equations (3.7.1), (3.7.3 


ae: 


and (3.7.5). 


the severity of the shock is almost ten times as great as it would be if the 
additional 10°% of space were available. 


3.8 GENERAL INFLUENCE OF SHAPE OF ACCELERATION-TIME CURVE 
ON AMPLIFICATION FACTOR 


When amplification factor curves are not available for a special shape of 
acceleration-time curve, an approximate value of A,, may be obtained by 
interpolation between or extrapolation from the curves of the preceding 
sections. The shape of the acceleration-time curve and its duration (r2) 
or frequency (w2) should be found, first, by the methods described in Part II. 
The shape found should then be compared with the standard shapes shown 
in Part II, for which amplification factors are given in Part III. 

The amplification factor found in this way will generally be within 25% 
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of the true value because amplification curves for pulse accelerations do not 
differ greatly even for very different acceleration-time curves as long as the 


5 
or 
We 


Fig. 3.8.1—Dependence of amplification factor on shape of symmetrical acceleration pulse. 


2. 


Fig. 3.8.2—Effect of asymmetry of an acceleration pulse on amplification factor. 


amplitudes and frequencies are adjusted to the same scales. This is illus- 
trated in Fig. 3.8.1 where the amplification factor curves are drawn for 
square wave, half-sine wave, triangular and cubic pulses. 
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Amplification factors for small values of w:/w2 may be calculated very 
accurately if it is observed that the initial slope of the amplification factor 
curve for a pulse acceleration is proportional to the area under the accelera- 
tion-time curve. For example, noting that the initial slope of the amplifica- 
tion factor curve for the half-sine wave pulse is 2, we assign the value 2 to 
the area under the half-sine wave. On the same scale, the area under a | 
square wave pulse is x and under a triangular pulse is 7/2. Accordingly, 
the initial slopes of the amplification factor curves for the latter two pulses 
are x and 2/2 respectively. 

As an additional aid in finding amplification factors for unusual cases, 
Fig. 3.8.2 is given to show the effect of asymmetry of an acceleration pulse. 
The pulse is triangular in shape but the time (rp) taken to reach the peak 
value of acceleration may have any value from zero to the total duration 
(r2) of the pulse. 


PART IV 


DISTRIBUTED MASS AND ELASTICITY 


4.1 INTRODUCTION 


It is important to be aware of the conditions under which the assumption 


of lumped parameters is permissible. In Parts I and II the cushioning 
medium was assumed to be massless, so that wave propagation (or surges) 
through it was ignored. Such surges will contribute to the acceleration 
imposed on the packaged article and we should be able to predict both the 
magnitudes and frequencies of the additional disturbances. If this is done, 
the information in Part III may be used to obtain at least a rough estimate 
of the resulting effects. In Part III itself the effects of accelerations were 
determined by studying the response of a system having only one degree oi 
freedom; that is, an element of the packaged article was assumed to be a 
single mass supported by a massless spring. Every real element, of course, 
has an infinite number of degrees of freedom, so that it is important to 
discover the contribution, of the higher modes of vibration of an element, 
to the overall response. 

Both of these problems (distributed parameters of mass and elasticity in 
the cushioning medium and in an element of the packaged article) are 
studied in this part. One example of each type is considered, and the choice 
of the example in each case was influenced by considerations of expediency, 
namely that the mathematical derivations be relatively simple and lead to 
solutions for which not too lengthy computations are necessary to yield 
results that can be applied practically. At the same time, the examples 
chosen are believed to give some insight into several of the physical phe- 
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nomena involved. The treatment is by no means complete, but a more 
detailed investigation is beyond the scope of this paper. 


4.2 Errect oF DISTRIBUTED MASS AND ELASTICITY OF CUSHIONING 
ON ACCELERATION OF PACKAGED ARTICLE 


Referring to Fig. 4.2.1, we consider the packaged article, of mass m2, 
to be supported by distributed cushioning of mass m, and depth ¢. The 
cushioning may be a pad, say of rubber, in which case ¢ is the pad thickness, 
or it may be a helical metal spring, in which case ¢ is the coil length. The 
package is dropped vertically from a height / and has attained a velocity v 
at the instant of contact (¢ = 0) of the outer container and the floor. The 
outer container is assumed to be heavy enough so that there is no rebound. 
A horizontal plane in the cushioning is located by a coordinate x measured 
from the end of the cushioning attached to the outer container. The vertical 


} M2 | 
i 4 re 


MMe: | i | £ (cushion) 


= x 





| 


| 
| Floor 
VID) | PITT 7 LER E ELEM E A 
Fig. 4.2.1—Packaged article of mass m2, supported on distributed cushioning of depth ¢ 
and mass m., depicted at the instant of first contact of the outer container (m:) and 
the floor. 








displacement of the plane x is designated by u. The undamped motion of 
the cushioning after contact is governed by the one-dimensional wave 


equation: 
au 20 u 
— = @ —— 4.2.1) 
ar Ox?’ ( 
in which a is the velocity of propagation of longitudinal waves in the cushion- 
ing. If the cushioning is continuous, 


¢=-, (4.2.2) 
p 


where E is the modulus of elasticity and p is the density of the cushioning. 
If the cushioning is a helical spring, 
P= — (4.2.3) 


where & is the spring rate. 
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The initial conditions of the system are 


[u} <0 anal 0, 


I = 
ot t=0 


The boundary conditions are 


(tt]en0 = 


he ou » ie Ou 
Ox |rmf “1 Ot? |rmt 


Equation (4.2.7) expresses the requirement that the force on the upper end 
of the cushioning must balance the inertia force of the packaged article. 
For continuous cushioning kf should be replaced by EA, where A is the 
cross-sectional area of the cushioning. 

A solution of (4.2.1) satisfying conditions (4.2.4) and (4.2.6) is 


x 
 n® 
u = >> A, sin = sin w, t, (4.2.8) 


n=l 


‘ root of a transcendental equation to be obtained from 


where w,, is the 7 

(4.2.7) and A, is a constant to be determined by (4.2.5). Substituting 

(4.2.8) in (4.2.7) and equating coefficients of like terms of the series, we 
obtain the transcendental equation 

worl wr mM. 

tan = — 


: (4.2.9) 
a a My 


Substituting (4.2.8) in (4.2.5) we obtain, by the usual methods of expansion 
into trigonometric series, 


2v 


A,=- (-2! we 3) (4.2.10) 
Wn 2 sin —— 
a a 


Hence the complete solution of the problem is 


dn cin, nk .- 
™ 2v sin - SsIn w, lt 


at (~ +i dnt 
a a 


Our chief interest is in the acceleration of mz. Making use of (4.2.7) 
and (4.2.9) we find, from (4.2.11), that this acceleration is 


k “| _ ates Dd Ba sin wnt, (4.2.12) 


(4.2.11) 


EYE 


n=l 





yer end 
article, 


is the 


(4.2.8) 
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) /™ (™: 4 en =) 
ms \m5 a” 
Bn = ie Se 
M. , Me , Wt 
-+— + — 
Mo ‘ a* 

The acceleration of my is, therefore, a sum of sinusoids of frequency w, 
and amplitude tw B,. Now, two is the maximum acceleration that my» 
would attain if the mass of the cushioning were negligible. Calling G, 
the maximum acceleration in the m‘® mode and Gp the maximum accelera- 
tion neglecting the mass of the cushioning, as in Part I, we have 


Gr - 
— = B,. wh 
G (4.2.15) 


But B, depends only on the ratio m./mz, as may be seen from equations 
(4.2.9) and (4.2.14). Similarly the ratio of the frequency (w,) of any mode 
to the frequency (wo) with massless cushioning depends only on m,/mz , 
as may be seen from equations (4.2.3), (4.2.9) and (4.2.13). Hence, both 
the amplitude and frequency ratios for the acceleration in any mode depend 
only on the ratio of the mass of the cushioning to the mass of the packaged articie. 
The ratios G,,/Go and w,,/w) are plotted against m./mz in Figs. 4.2.2 and 4.2.3 
for the first five modes. Jt may be seen from these figures that the accelerations 
in the higher modes can be very important. For example, if the cushioning 
weighs half as much as the packaged article the maximum acceleration in 
the second mode is about 40% of the acceleration in the first mode and the 
latter is about the same as found by the elementary method of Part I. 
This could have a disastrous effect on an element of the packaged article if 
the latter had a fundamental frequency near that of the second mode of the 
cushioning, the latter being found, from Fig. 4.2.3, to be about five times 
the fundamental frequency of the package. 

It must be remembered that damping has been neglected in the above 
investigation and damping in the cushioning will serve to mitigate the se- 
verity of the higher mode accelerations to a great extent. However, the 
danger is always present at the start of a design and the possibilities of un- 
favorable combinations should be studied in every case. 
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Fig. 4.2.2—Influence of ratio of mass of cushioning (m,.) to mass of packaged article 
(m2) on acceleration ratio. The numerator of the acceleration ratio is the maximum 
acceleration (G,) in the nt mode of vibration transmitted through the cushioning. The 
denominator of the acceleration ratio is the maximum acceleration (Gp = V 2hk2/mag) that 


the mass mz would experience if the mass of the cushioning were negligible. See equations 
(4.2.15), (4.2.14), (4.2.9). 


20 festesst: 
18% 


16> 


© O02 04 06 08 10 12 14 16 18 20 


Me 


m 
2 
Fig. 4.2.3—Influence of ratio of mass of cushioning (m.) to mass of packaged article 
on frequency ratio. The numerator of the frequency ratio is the frequency (wn) of the 
nt» mode of vibration transmitted through the cushioning. The denominator of the 
frequency ratio is the frequency (w, = k2/m:) of vibration of the mass ms neglecting 
the effect of the mass of the cushioning. See equations (4.2.9), (4.2.13), and (4.2.3). 











article 
ximum 
. The 
g) that 
lations 


ticle 
the 
the 
ting 
3). 








DYNAMICS OF PACKAGE CUSHIONING 451 


4.3 EFFECT OF DISTRIBUTED MAss AND ELAstTIcitTy, OF AN ELEMENT 
OF THE PACKAGED ARTICLE, ON THE AMPLIFICATION FACTOR FOR 
A HA.r-StnE-WAVE PULSE ACCELERATION 


In this section we shall determine the contribution of the higher modes of 
vibration of a structural element to its total response to a half-sine-wave 
pulse acceleration. For the shape of the element, we choose a prismatic 
bar because this leads to the simplest mathematical formulation of the 
problem and such a bar is also a common structural element. Other con- 
siderations influence the choice of direction of acceleration with respect to 
the axis of the bar. The transverse direction (cantilever) is the most 
practical from a physical standpoint, but, for purposes of comparison with 
the one-degree-of-freedom system, the parallel (axial) direction of accelera- 
tion is the more logical. Both problems lead to solutions in the form of 
infinite series, but, in the latter case, the expression for the strain at a fixed 





x L 
-V 
2 
= ko 
= 
a as 
Fig. 4.3.1—The system studied in Section 4.3 depicted at the instant of contact with 
the floor. 


end can be summed in terms of elementary functions without difficulty. 
Since it is necessary to determine maximum values of strain over a wide 
range of frequency ratios for the plotting of an amplification factor curve, an 
enormous reduction in the time required for accurate computations is 
obtained by choosing the axial case. Furthermore, the axial case appears to 
contain the essential features which might result in differences between the 
response of a one-degree-of-freedom system and a continuous one. 

The complete system to be studied is illustrated in Fig. 4.3.1. To the 
mass m;, supported on massless cushioning of constant spring rate ke, 
is attached one end of an elastic prismatic bar, of length ¢, cross sectional 
area A, modulus of elasticity E, and density p, with its axis oriented verti- 
cally. The system is dropped from a height / so that its velocity is 2 
at the instant of contact of the cushioning with the floor. The mass of the 
bar is supposed to be small in comparison with m2 and perfect rebound is 
assumed, so that the motion of m, during contact is a half-sine wave of 
frequency 
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ws = 4/®. (4.3.1 
Mo 


The maximum acceleration of ms is thus tw:. If this magnitude of ac 
celeration were reached very slowly, so as not to excite transient longitudinal 
waves in the bar, the maximum force between the bar and mz would be the 
product of the acceleration and the mass of the bar: 

F = wopAl. (4.3.2) 


Hence the strain at the end of the bar attached to ms would be 


We f 
@ = a (4.3.3) 
Our problem is to find the ratio of the maximum transient strain to ¢. 


18 ; 


et ae Oe rh 





O/ . Ma 
t/2 2Wol 


Fig. 4.3.2—Amplification factors for an element of the packaged article having dis- 
tributed mass and elasticity. The package has linear undamped cushioning and perfect 
rebound. See equations (4.3.15), and (4.3.14). 

Let u be the displacement of a transverse plane section of the bar distant x 
from the end attached to m,.. Then, the equation of motion of the bar is 


Ou OU 


oe ws rs 
_ (4.3.4) 


ax?’ 


where a is the velocity of propagation of longitudinal waves in the bar: 


aq =, (4.3.5) 


* 


Taking the instant of first contact of the cushioning with the floor to be 
t = 0, we know, from Part II, that the system will leave the floor when ¢ = 





and 








ac 
nal 
the 


io 





DYNAMICS OF PACKAGE CUSHIONING 453 


t/wo. We shall therefore treat separately, as in Part III, the motion 
during contact 


= = 7T 
Q<=t<— 
Wwe 
and after rebound 
= 
ie ; 
we 


During the first interval, the initial and boundary conditions are 


[1] 4.0 = Q, (4.3.6) 
Ou i 
E = —0, (4.3.7) 
Ol Jr=0 : 
[ulzn0 = —— sin wel, (4.3.8) 
we 


, 
Hi = 0. (4.3.9) 
OX _|r-1 


The first and second conditions state that, at the instant of contact, all 
points in the bar are moving with the approach velocity v, without relative 
displacement. The third condition prescribes the half-sine wave motion 
of the end of the bar that is attached to m,. The fourth condition states 
that the strain at the free end of the bar is always zero. 
By the usual methods, a solution of (4.3.4) satisfying conditions (4.3.6) 
to (4.3.9) is found to be 
v cos (€ — x) sin wol sin ~~ sin = 


8e = ) 2 
sc lp calla 4. z. a (43.10) 


wet Td a=l,3b ° nra \* 
@2 COS —— n ~}—1 
a 2uol 


(0 212 >, 1). 

a. url 
The displacement is seen to be a forced vibration at the frequency (we) 
of the applied acceleration, on which are superposed the free vibrations of 
the bar given by the series expression. The frequency of the fundamental 
mode of vibration of the bar is ra/2¢ and the frequencies of the higher modes 
are the odd integral multiples of the fundamental. 
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The strain at the attached end of the bar is 


. nat 
sin => 


2¢ 


wel .. 4 o 
@ tan = a mt T n= te. ted - J (4.3.11) 


2w2 £ 


= Ww oe 


It may be verified that the sum of the series in (4.3.11) is given by 





; nrat 
ra sin 
ne ot = tan = SiN wet + COS wot — 1, (4.3.12) 


(0<:<%). 


It should be observed that the summation is valid only in the interval 0 < 
t < 2¢/a. However, the series is periodic with half period 2¢/a and includes 
only the odd terms, so that the function repeats itself with reversed sign 
after each interval 2¢/a. Hence the summation, valid for all ¢, can be 
written 

nnat 


sin 2 


4 wD 

x Pi... P (ara ] 
wot 6) 

= (-1)* [tan oe in we (: - = is + co son ( - 2H) _ 1] (4.3.13) 


9 
k =m _ when = fie _—— 


m = OQ, 1, 2, 3, - 
We may, therefore, rewrite (4.3.11) in the form 


= = ~tan “sin ud + (~i1P tan on in 2 ( _ a) 


a 


+ cos o2(1 _ 2m) _ 1] (4.3.14) 
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k = m when ‘a “t< =e 


m = 0,1, 2,3--- . 


The expression (4.3.14) is simple enough so that the maximum value (€m) 
of the strain at the attached end can be obtained without difficulty for any 
ratio of the fundamental frequency (w; = wa/2¢) of the bar to the frequency 
(we) of the disturbing acceleration. The amplification factor 


ém | E 2w; 4 = 
: - = a Em | 5 (4.3.15) 
€ vwopl Tw. V 


Em | 


An = 


may then be calculated. The results of these calculations are plotted in 
Fig. 4.3.2. The important feature of this curve is that the amplification 
factor is everywhere less than the corresponding amplification factor for the 
one-degree-of-freedom system (Fig. 3.2.2, 8; = 0). Hence the assumption 
of lumped parameters is on the side of safety as regards amplification factor. 

It is interesting to observe that the curve of A, vs. w:/we2 , for this case, 
is a straight line between w:/w2 = 0 and w:/w2 = 1. This arises from the 
fact that, for w:/w, < 1, equation (4.3.14) reduces to 


< = cos wt — 1, 51). (4.3.16) 


Hence, when the duration of shock is less than the half period of the funda- 
mental mode of vibration, the maximum value of strain occurs at the end of 
impact and is equal to twice the ratio of the approach velocity to the velocity 
of wave propagation in the bar. 

The whole solution of the problem is not yet completed; for, although it is 
fairly evident from the fact that there is at least one maximum in the inter- 
val 0 = ¢ = w/w for all values of w:/we , it must be verified that the maxi- 
mum strain (and, therefore, the amplification factor) is never greater after 
t = m/we than before. Defining a new time coordinate 


(4.3.17) 


we have, for the initial and boundary conditions of equation (4.3.4) for 
tS r/o, 
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_ nar | nnx 
oan. . Sa ~ 


¢ = 2, 
lek = Sof i Qwol t 


Ta n=1,3,5 9 na : 1 
nN Sa - 
2we { 


nar | mrx 


We ‘ 
au ‘ won aa 
= - (4.3.19) 
t’ =0 my n »3,5° 


ot nia \* 1 
"1 \2et) — 


[u]e.o = vt’, (4.3.20) 


Ou 
= Qj; 3.2 
Bia (@ 


The first and second conditions state that the displacement and velocity of 
every point in the bar must be the same at the beginning of the second inter- 
val as at the end of the first interval; the expressions in (4.3.18) and (4.3.19) 
are obtained from (4.3.10). The third condition prescribes the constant 
velocity of departure from the floor of the mass ms and, therefore, of the 
end of the bar attached to it. The last condition states, again, that the 
strain at the free end of the bar is zero. 

It may be verified that a solution of (4.3.4) satisfying conditions (4.3.18) 
to (4.3.21) is 


_ , mex , {nrat’ 
C, sin sin - + 


f 2¢ 


where 


Sef sin Deal 


C, sin yn = a 
3T an i& :) _ | 

gor (1 nan 

8v to. -— 

+ 2w { 


3 2.8 nna ¥ 
ed 1 OY) li (4.3.24) 


(4.3.23) 


Cr COS Yn 


una 


2we sf 
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Hence, the strain at the attached end of the bar is 


_ nrat 
Ou 40 = O28 
nie. .2 ¥ 
r=) 


Ox TA n=1,3.5 nna . 1 
nN maa 
2w»2 f 
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The two series may be summed, as before, with the result 
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Once more, the expression for the strain at the attached end of the bar is 
in a form suitable for rapid calculation and it can be shown the ¢€ in equation 
(4.3.26) for ¢ > m/w is never greater than the ¢ in equation (4.3.14) for 0 = 
t S m/we for the same w,/w2. Hence, Fig. 4.3.2 and the conclusions follow- 
ing equations (4.3.15) and (4.3.16) need not be modified. 


NOTATIONS 


Cross sectional area of a bar element of the packaged article. Also, a 
constant of integration. 


Amplification factor when the reference acceleration is Go. Ratio of 
maximum dynamic response to the response to a slowly applied ac- 
celeration of magnitude G,g. 


Amplification factor when the reference acceleration is Gn. Ratio of 
maximum dynamic response to the response to a slowly applied accelera- 
tion of magnitude Gng. 


In Section 1.15, the sum of all the trapezoidal areas from x. = 0 to x2 = 
(x2)n. Also, in Section 4.2, the coefficient of the nth term of a series. 


The area of a trapezoid with altitude A(x2), and sides P,_; and Pp. 


xo/l in the tension spring package. Also, in Part IV, the velocity of 
propagation of longitudinal waves. 
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A parameter of cushioning with cubic elasticity defined in equatio; 
(1.5.3). Also, a constant of integration. 


Coefficient in the n‘* term of a series. 

f/l in the tension spring package. 

A constant of integration. 

Coefficient of the nt» term of a series. 

A constant defined in equation (1.7.11) 

Damping coefficient of an element of a packaged article. 
Damping coefficient of linear cushioning. 

The elliptic cosine function. 


Hypothetical displacement that would result if initial spring rate were 
maintained. 


Maximum possible displacement of packaged article in cushioning with 
tangent elasticity. 


Maximum displacement of packaged article. 
Value of dm when ko = ko. 


Displacement of bi-linear cushioning at which the spring rate changes 
from ko to kp. 


Modulus of elasticity. 


In the tension spring package the stretch of a spring when the displace- 
ment is dm. 


e ‘’, where e is the Naperian base 2.718: -- 
In section 2.7, a frictional force. 
In the tension spring package, the maximum force on a spring. 


In the tension spring package, the difference between / and the distance 
between hooks of an unstretched spring. 


Frequency of vibration of an element of the packaged article. 


Frequency of vibration of the packaged article on its cushioning. 


Hypothetical maximum acceleration (in number of times g) that would 
result if initial spring rate were maintained. 


AmGm or AoGo, i.e. the slowly applied acceleration (in number of times g) 
that will produce the same maximum response as a transient acceleration 
of maximum value G,», or Go. 


Maximum acceleration (in number of times g) in cushioning with fric- 
tion and spring rate kp. 


Absolute value of maximum acceleration of packaged article in units of 
“number of times gravitational acceleration.” 


Value of G, when ko = ki. 





K,, Ko, Kz 


DYNAMICS OF PACKAGE CUSHIONING 459 


In section 1.15, the maximum acceleration (in number of times g) ex- 
perienced by the suspended mass when dropped from a height A». In 
Part IV, the maximum acceleration (in number of times g) of the nt 
mode of vibration. 

Maximum acceleration (in number of times ¢) after rebound. 

Safe value of G,. 

Gravitational acceleration. 


Height of drop. 


In Section 1.15, the height of fall that will cause the cushioning to dis- 
place an amount (x2) n. 


In the tension spring package, the initial spring rate of the suspension. 
In Section 2.8, the complete elliptic integral of the first kind. 


The initial spring rates in the three mutually perpendicular directions 
normal to the faces of the package frame. 


In the tension spring package, the spring rate of a spring. In Section 
2.8, the modulus of an elliptic integral. 


In Section 4.3, 0, 1, 2,3, ---. 
Initial spring rate of non-linear cushioning. 
Optimum value of initial spring rate ko. 


Spring rate of lumped elasticity of element of packaged article. 


Spring rate of linear cushioning. 
Spring rate of bilinear cushioning after bottoming. 
Spring rate defined in equation (2.7.7). 


Constant defined in equation (1.8.2). 


In the tension spring, the projection of /; on a horizontal plane. In 
Section 4.2, length of cushioning. In Section 4.3, length of element of 
packaged article. 


In the tension spring package, the distance between the two support 
points of a spring when the suspended article is in the equilibrium 
position. 

Constant defined in equation (1.8.4), equal to Gm /Go. 

Reduced mass defined in equation (2.4.5). 

In Section 4.3, 0, 1, 2,3, --:. 

Lumped mass of element of packaged article. 

Lumped mass of packaged article. 


Lumped mass of outer container. 


Mass of cushioning. 


M?. 





BELL SYSTEM TECHNICAL JOURNAL 
0, 1, 2,3, .--° 
Force transmitted through cushioning. 


Asymptotic value of force transmissible through cushioning with hyper - 
bolic tangent elasticity. 


Maximum force exerted on packaged article by cushioning. 
In Section 1.15, the load that produces displacement (x2) n. 
Force between package and floor. 


Coefficient of cubic term in load-displacement function for cushioning 
with cubic elasticity. 


The direction cosines of the acceleration direction with respect to the 
normals to the faces of the package frame. 


The elliptic sine function. 
The period of vibration of the packaged article on its cushioning. 


Time coordinate. 
T 
we 


Time of first contact of package with floor. 

Time at which maximum displacement or acceleration occurs. 
Time at which package leaves floor on rebound. 

Time at which the displacement reaches d,. 

Displacement in x direction. 

Approach velocity. 

Weight of packaged article. 

Weight of outer container. 


Xx, — X2; relative displacement of m, with respect to mo. 


Relative displacement of m, with respect to m2 at time ?’. 


In the tension spring package, the perpendicular distance from an inner 
spring support point to the nearest plane, perpendicular to the displace- 
ment direction and containing four outer spring support points. 


Displacement of my. 
Velocity of m. 
Acceleration of m). 
Displacement of mz. 


Velocity of mz. 





Xe 
Xmax 


(x2) n 
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Acceleration of m3. 


Maximum value of x. 


In Section 1.15, the displacement associated with the n*® point. 


The value x would have if the acceleration reached its maximum value 
in a very long time. 


In Section 1.15, equals (x2), — (2) n-1. 
X2_X3. 
x2/l (tension spring package). 

Phase angles. 

Fraction of critical damping of an element of the packaged article. 
Fraction of critical damping of package cushioning. 

Phase angle of nt» term of series (equation (4.3.22) 

Strain at attached end of element under transient conditions. 

Strain at attached end of element under non-transient conditions. 
Maximum strain at attached end of element under transient conditions. 


Angle between the displacement direction and the acceleration di- 
rection. 


3.14159--+. 

Density (mass per unit of volume) 

Pulse duration of a half-sine-wave acceleration. 

Pulse duration associated with non-linear cushioning. 

Duration of bottoming of cushioning with bi-linear elasticity. 

Time required to reach peak value of a triangular acceleration pulse. 
Radian frequency defined in equation (2.4.6). 

Radian frequency of vibration of an element of the packaged article. 
Radian frequency of vibration of damped element of packaged article. 
Radian frequency of vibration of the packaged article on its cushioning. 


Radian frequency of vibration of the packaged article on damped 
cushioning. 


A frequency defined in equation (2.10.10). 
A frequency defined in equation (2.8.8). 


Radian frequency of nth mode. 





Abstracts of Technical Articles by Bell System Authors 


Dimensional Stability of Plastics. RoBert Burns. Because of inherent 
insulating properties, rigid plastics play an important part in the design and 
manufacture of precision electrical apparatus. Almost invariably, practical 
design considerations require that the plastics have reasonable structural 
possibilities since it is rarely practicable to disassociate completely electrical 
and structural functions. 

This paper discusses one of the important factors in the successful use of 
plastics in precision devices, namely, dimensional stability. Since plastics 
are organic compounds, one must be prepared to accept a degree of insta- 
bility not usually encountered in metals. The measurement of this property 
is therefore of prime importance to the user of plastics since the data provide 
a basis for design adjustment which frequently is the difference between 
failure and success. 

The various types of dimensional change are reviewed. Data illustrating 
the separate effects of humidity, drying, and cycling procedures are sub- 
mitted. The influence of fabricating processes such as compression or 
injection molding, and sheeting, is included. 

Some Numerical Methods for Locating Roots of Polynomials? THuorNton 
C. Fry. It is the purpose of this paper to discuss the location of the roots 
of polynomials of high degree, with particular reference to the case of com- 
plex roots. This is a problem with which the Laboratories has been much 
concerned in recent years because of the fact that the problem arises rather 
frequently in the design of electrical networks. Attention is not given to 
strictly theoretical methods, such as the exact solution by elliptic or auto- 
morphic functions: nor to the development of roots in series or in continued 
fractions, though such methods exist and one at least—development of the 
coefficients of a quadratic factor—is of great value in improving the accuracy 
of roots once they are known with reasonable approximation. 

Instead, the paper deals with just two categories of solutions: one, the 
solution of the equations by a succession of rational operations, having for 
their purpose the dispersion of the roots; the other, a method depending on 
Cauchy’s theorem regarding the number of roots within a closed contour. 

Thermistor Technics2 J.C. JouNsON. This paper is confined to a study 
of how the three basic types of thermistors, namely, externally-heated or 
ambient temperature type, the directly-heated type, and also the indirectly- 

1A.S.T.M. Bulletin, May 1945. 


2 Quarterly Applied Mathematics, July 1945. 
® Electronic Industries, August 1945. 
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heated type, are used in simple feedback amplifiers as regulation and control 
devices to effect the economies inherent in an entirely electrical system by 
eliminating such mechanical devices as motor-driven condensers, sliding 
contacts and rotary switches. 

Dynamic Measurements on Electromagnetic Devices* E.L. Norton. A 
method is presented by which measurements of flux may be made at any 
desired time during the operate cycle of an electromagnet. Apparatus is 
described which operates the magnet cyclically at an accurately held rate, 
and provides a means for measuring flux either by the use of a search coil 
or by the operating winding of the magnet itself. When using a search coil, 
it is connected to a direct-current milliammeter at the time in the cycle at 
which the value of the flux is desired and disconnected at the end of the 
cycle or just before the magnet is energized for the next pulse. If proper 
precautions are taken, the steady reading of the instrument is an accurate 
measure of the difference in the flux in the coil between the time it is con- 
nected to the meter and the time it is removed, or, since the latter is zero 
except for residual flux, the reading is a direct measure of flux. 

The same apparatus may be used for the measurement of instantaneous 
current by the addition of an air core mutual inductance, and its use is 
extended to the measurement of armature position and velocity by the 
addition of a photoelectric cell! and the proper amplifiers. 


A form of vacuum tube filter is described which effectively filters the 
pulses from the indicating instrument without affecting the accuracy of the 
measurements. 


Coaxial Cables and Television Transmission.6 HarotD S. OSBORNE. 
Communication techniques and facilities useful to the entertainment 
industry have evolved naturally from the Telephone Companies’ main 
objective—the transmission of speech. The development of carrier sys- 
tems for long-distance transmission and technical features involved in the 
latest carrier medium—the coaxial cable—are reviewed. The television 
transmission capabilities of this medium, both now and what may be 
expected shortly after the war, are mentioned. The extensive system of 
such cables planned for the next five years, supplemented by radio relay 
systems to the extent that these prove themselves as a part of a communica- 
tions network, will provide an excellent beginning for a nation-wiue tele- 
vision transmission network. Planned primarily to meet telephone 
requirements, this network of cables will be suitable to meet the transmission 
needs of the television industry. 

The Performance and Measurement of Mixers in Terms of Linear-N etwork 
Theory® LL. C. PETERSON AND F. B. LLEWELLYN. This paper discusses 

4 Elec. Engg., Transactions Section, April 1945. 


5 Jour. S.M.P.E., June 1945. 
6 Proc. I.R.E., July 1945. 
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the properties of mixers in terms of linear-network theory. In Part I the 
network equations are derived from the fundamental properties of nonlinea: 
resistive elements. Part II contains a résumé of the appropriate formulas 
of linear-network theory. In Part III the network theory is applied, first 
to the case of simple nonlinear resistances, and next to the more genera! 
case where the nonlinear resistance is embedded in a network of parasiti: 
resistive and reactive passive-impedance elements. In Part IV application 
of the previous results is made to the measurement of performance properties 
The “impedance” and the “incremental” methods of measuring loss are 
contrasted, and it is shown that the actual loss is given by the incremental 
method when certain special precautions are taken, while the impedance 
method is in itself incomplete. 

A Figure of Merit for Electron-Concentrating Systems.’ J. R. Pierce. 
Electron-concentrating systems are subject to certain limitations because 
of the thermal velocities of electrons leaving the cathode. A figure of merit 
is proposed for measuring the goodness of a device in this respect. This 
figure of merit is the ratio of the area of the aperture which, in an ideal 
system with the same important parameters as the actual system, would 
pass a given fraction of the cathode current to the area of the aperture which 
in the actual system does pass this fraction of the cathode current. Ex- 
pressions are given for evaluating this figure of merit. 


A 60-Kilowatt High-Frequency Transoceanic-Radiotelephone Ampblifier.® 
C.F. P. Rose. Here is described a high-frequency radio amplifier recently 
developed for the transoceanic-telephone facilities of the Bell System at 
Lawrenceville, New Jersey. In general, the amplifier is capable of delivering 
60 kilowatts of peak envelope power when excited from a 2-kilowatt radio- 


“ec 


frequency source. It is designed to operate as a “class B’”’ amplifier for 
transmitting either single-channel double-sideband or twin-channel single- 
sideband types of transmission. Features are described which permit 
rapid frequency-changing technique from any preassigned frequency to 
another lying anywhere within the spectrum of 4.5 to 22 megacycles. 

Some Notes on the Design of ElectronGuns.2 A. L.SAmMuEL. A method is 
outlined for the design of electron guns based on the simple theory first 
published by J. R. Pierce. This method assumes that the electrons are 
moving in a beam according to a known solution of the space-charge equa- 
tion, and requires that electrodes exterior to the region of space charge be 
shaped so as to match the boundary conditions at the edge of the beam. 
An electrolytic tank method is used to obtain solutions for cases which are 
not amenable to direct calculation. Attention is given to some of the 

7 Proc. I.R.E., July 1945. 


8 Proc. 1.R.E., October 1945. 
’ Proc. I.R.E., April 1945. 
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complications ignored by the simple theory and to some of the practical 
difficulties which are encountered in constructing guns according to these 
principles. An experimental check on the theory is described, together 
with some information as to the actual current distribution in abeam 
produced by a gun based on this design procedure. 

Microwave Radiation from the Sun!° G. C. Sournworrn. During the 
summer months of 1942 and 1943, a small but measurable amount of micro- 
wave radiation was observed coming from the sun. This appeared as ran- 
dom noise in the outputs of sensitive receivers designed to work at wave- 
lengths between one and ten centimeters. Over a considerable portion of 
the range, the energy was of the same order of magnitude as that predicted 
by black-body radiation theory. 

Attempts were made to determine the effect of the earth’s atmosphere on 
this radiation. Measurements made near sunrise or sunset, when the 
path through the earth’s atmosphere was relatively long, differed only 
slightly from those made at noon. This suggested that any absorption 
that may have been present was small. In this connection it is of interest 
that small temperature differences could be noted between points below the 
horizon and the sky immediately above. This also suggested that the 
earth’s atmosphere was relatively transparent. 

In another kind of measurement the parabolic receiver was centered on 
the sun and its output was observed as the sun’s disc moved out of the aper- 
ture of the receiver. The directional pattern so obtained indicated that 
at the shorter wave-lengths the sun’s apparent diameter was considerably 
larger than that measured by ordinary optical means. This suggested that 
there may have been some refraction or perhaps scattering by the earth’s 
atmosphere. 

Resistive Attenuators, Pads and Networks—An Analysis of their A pplica- 
tions in Mixer and Fader Systems (Part Eight of a Series)’ Pavur B. 
Wricut. In last month’s discussion, the series-connected fader and the 
parallel-connected fader systems were considered, together with an analysis 
of their performance expressed both algebraically and in terms of the 
hyperbolic functions of a real variable. In this installment, the series- 
parallel-connected fader system discussion is continued and equations 
describing the complete behavior of this type network system are developed. 
This is followed by further analytical work dealing with the parallel-series- 
connected fader and mixer system and several lesser known systems which 
are quite useful to use. These are the multiple bridge and the lattice network 
systems which may be utilized to advantage for some applications. All of 

10 Jour. Franklin Institute, April 1945. 


4 Communications, September 1945. (Preceding parts of this Series appeared in earlier 
issues of Communications.) 
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the equations which are derived are shown in the algebraical, hyperbolica! 
and symbolical forms. The key chart which was presented earlier in this 
series may be used to great advantage when checking the definitions of th: 
symbols used which are not specifically defined in the text. This procedure 
also may be directly applied to the hyperbolic equations shown. It is of 
course necessary to take into account that, in general, subscripts are used 
in most of the equations in the text while the key chart does not have any 
subscripts. This does not, however, alter the fundamental forms nor their 
definitions in terms of the propagation function, theta. To avoid the neces- 
sity for extensive interpolation of the hyperbolic function tables to find the 
correct numerical values for the various functions used throughout the text, 
a series of tables providing all of the functions required is presented. 


RAY 
C.E., 
Profes 
ing, ' 
Resea 
1940- 
mathe 
tories 
a 
exper 


J. 
Tec 
Eng 


A 


in F 
Add 


tor 

Tecl 
bee! 
prin 
ultr 





Contributors to this Issue 


Raymonp D. Minpiin, B.A., Columbia University, 1928; B.S., 1931; 
C.E., 1932; Ph.D., 1936. Assistant 1932-38, Instructor 1938-40, Assistant 
Professor 1940-45, Associate Professor 1945-, Department of Civil Engineer- 
ing, Columbia University. Consultant, Section T, National Defense 
Research Committee (later Office of Scientific Research and Development), 
1940-45, on the development of the rugged radio proximity fuze and on 
mathematical problems in fire control. Consultant, Bell Telephone Labora- 
tories, 1943-44, Member of the Technical Staff 1944-45, Consultant 1945-. 
Dr. Mindlin has been concerned with the fields of mathematical and 
experimental mechanics. 


J. R. Preece, B.S. in Electrical Engineering, California Institute of 
Technology, 1933; Ph.D., 1936. Bell Telephone Laboratories, 1936-. 


Engaged in study of vacuum tubes. 


A. L. Samuet, A.B., College of Emporia (Kansas), 1923; S.B. and S.M. 


in Electrical Engineering, Massachusetts Institute of Technology, 1926. 
Additional graduate work at M.I.T. and at Columbia University. Instruc- 
tor in Electrical Engineering, M.I.T., 1926-28. Mr. Samuel joined the 
Technical Staff of the Bell Telephone Laboratories in 1928, where he has 
been engaged in electronic research and development. Since 1931, his 
principal interest has been in the development of vacuum tubes for use at 
ultra-high frequencies. 





